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1  Summary 

This  document  is  the  final  report  for  the  ONR  project  N00014-90-J-1678,  the  third  phase  of 
a  study  entitled  “Modelling  bathymetric  control  of  near  coastal  wave  climate”  The  effort  in 
this  phase  of  the  study  was  concentrated  on  the  extension  of  the  angluar  spectrum  modelling 
technique  to  nonlinear  shallow  water  waves,  and  to  further  studies  of  the  Bragg  reflection 
mechanism  for  surface  waves  being  reflected  by  artificial  or  natural  bars. 

Data  from  a  laboratory  experiment  on  the  mean  surfzone  flows  associated  with  a  passively- 
driven  syphon  system  is  presently  being  reduced  and  will  be  described  in  a  separate  data 
report. 

2  Bragg  Reflection 

The  major  effort  in  this  area  was  related  to  describing  the  reflection  process  for  bars  with 
shapes  that  would  correspond  to  artificial  bars  to  be  installed  in  the  field.  Appendix  A 
contains  a  paper  which  describes  the  analysis  of  reflection  from  a  bar  field  using  the  Fourier 
decomposition  of  the  bar  shapes  in  the  analysis.  Appendix  B  describes  an  attempted  field 
installation  (see  Kirby,  1988)  of  bars  of  the  type  described  in  Appendix  A.  Finally,  Appendix 
C  assesses  the  usefullness  of  artificial  bars  as  coastal  protection  devices.  The  analysis  in 
this  paper  utilizes  the  bar  relection  program  developed  in  Kirby  (1987).  Finally,  Appendix 
D  presents  a  manuscript  which  addresses  some  theoretical  points  arising  in  the  analysis  of 
resonant  reflection  of  waves  by  sinusoidal  bars. 

3  Angular  Spectrum  Modelling 

The  development  of  the  angular  spectrum  methodology  for  modelling  nearshore  wave  prop¬ 
agation  was  extended  to  the  case  of  nonlinear  shallow  water  waves;  this  case  is  described  in 
Appendix  E.  The  work  leading  up  to  this  point  has  culminated  in  an  invited  review  article 
covering  work  on  linear  waves,  nonlinear  Stokes  waves,  and  shallow  water  Boussinesq  waves, 
which  is  presented  here  as  Appendix  F. 
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CHAPTER  58 

BRAGG  REFLECTION  OF  WAVES  BY  ARTIFICIAL  BARS 
James  T.  Kirby^,  Jeffrey  P.  Anton^ 


Abstract 

We  consider  the  extension  of  previous  theories  for  Bragg  reflection  of  sur¬ 
face  waves  by  parallel  bars  to  the  case  of  artificial  bars  placed  discretely  on  the 
seabed.  The  case  of  non-resonant,  weak  reflection  is  considered  first,  followed 
by  a  consideration  of  the  application  of  resonant  interaction  theory  to  the  dom¬ 
inant  Fourier  mode  of  the  bar  field.  Both  theories  are  compared  to  numerical 
results,  and  discrepancies  are  seen  in  both  cases.  Finally,  experimental  results 
are  compared  to  theory. 

Introduction 

The  discovery  that  the  Bragg  reflection  mechanism  leads  to  strong  reflection 
of  incident  surface  waves  by  periodic  bottom  undulations  has  led  to  speculation 
that  artificial  bars  could  be  constructed  which  would  partially  shelter  shores  or 
localized  structures  from  wave  attack.  Possible  bar  configurations  of  this  sort 
have  been  discussed  previously  by  Mei  et  al.  (1988)  and  Naciri  and  Mei  (1988). 
The  paper  by  Bailard  et  al.  (1990)  in  this  conference  describes  an  effort  which 
was  made  to  install  and  test  an  artificial  bar  field  offshore  of  a  natural  beach. 

The  purpose  of  the  present  study  was  to  extend  the  scope  of  available  theory 
and  techniques  which  were  available  for  predicting  wave  reflection  from  bars,  in 
support  of  the  proposed  field  study.  Here,  we  discuss  the  application  of  ana¬ 
lytic  perturbation  methods  for  both  non-resonant  and  resonant  cases.  We  also 
discuss  numerical  results,  which  point  out  limitations  present  in  both  analytic 
approaches.  Finally,  experimental  results  largely  provide  a  qualitative  verifica¬ 
tion  but  in  turn  show  some  limitation  of  the  small  amplitude  bar  theory. 

Theory  for  Small  Amplitude  Bars 

The  theory  which  provides  the  framework  for  analysis  here  is  given  by  an 
extended  mild-slope  equation  derived  by  Kirby  (1986). 

We  treat  the  w’ater  depth  A'(x,  y)  as  the  superposition  of  a  mildly-sloping 
bottom  A(x,y)  and  a  rapidly- varying  but  small-  amplitude  undulation  6(x,y): 

h\x,y)  =  b{x,y)-S{x,y)  (1) 

‘Assoc.  Prof.,  Center  for  Applied  Coastal  Research,  Dept,  of  Civil  Engrg.,  Univ.  of 
Delaware,  Newark,  DE  19716 

’Formerly,  Grad.  Stud.,  Coastal  and  Oceanographic  Engineering  Department,  University 
of  Florida,  Gainesville,  FL  32611 
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X 

Figure  1:  Bar  field  with  four  discrete  bars. 


Using  h{x,y)  as  the  reference  depth  in  the  mild-slope  sense,  the  model  equation 

v»  •  (CC.Vk#)  +  =  0  (2) 

is  obtained,  where  ^  is  the  value  of  the  linear  wave  potential  at  the  still  water 
surface.  The  model  coefficients  are  obtained  from 

w’  =  ^fctanh  5  ^  J  5  ^  (3) 

and  are  determined  by  the  value  of  h{x,y)  in  all  cases. 


In  the  absence  of  appropriate  field  data,  we  have  restricted  our  attention  here 
to  the  study  of  periodically  spaced  bars  {S  =  ^(x))  and  an  otherwise  uniform 
depth  h  =  constant.  In  principal,  6  is  arbitrary  aside  from  the  small  amplitude 
restriction.  In  the  present  study,  we  have  chosen  a  bottom  consisting  of  rectified 
sine  waves,  given  by 

D  cos  ^{x  -  N  L)  NL-^-<x<NL  +  ^ 

0  otherwise  ^  ^ 

TV  =  0,  . . . .,  iVk  -  1 

where  is  the  number  of  bars,  L  is  the  periodic  bar  spacing,  bi  is  the  footprint 
of  the  bar  on  the  bottom,  and  D  is  the  bar  height.  The  rectified  cosine  form 
is  chosen  mainly  for  its  convenience  in  later  analysis.  An  example  bar  field  is 
shown  in  Figure  1.  The  bax  field  is  periodic  over  intervals  of  width  L,  and  can 
be  conveniently  represented  by  the  cosine  series. 


=  51  ^nCOs(nAx) 
n=0 


where 


Do  = 


»(1  -  7.’) 


(1  +  cos  htt) 


(6) 
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Non-rssgpant  Rcflwtion 

For  the  case  of  h  constant,  the  model  equation  (2)  may  be  written  as 


4>  =  ■  iSVh<i>) 

(7) 

where 

Ak 

^  2kh  +  s\nh2kh 

(8) 

With  6{x)  representing  bars  varying  in  the  x-direction, 
waves  according  to 

we  represent  oblique 

<f>{x,y)  =  ^x)e*’"’' ;  m  =  ksm9 

(9) 

and  obtain 

^xx  +  =  Qf{^^x)x  ~  vn? ciS^  , 

(10) 

e  =  =  k^  cos^  e 

This  equation  has  been  obtained  by  Miles  (1981)  who  used  it  to  study  reflection 
from  a  single  isolated  obstacle. 

For  6(i)  confined  in  a  finite  region  of  space,  we  may  write 

^(x-foo)  =  Te"' 

^(x-*-oo)  =  +  (11) 

where  T  and  R  are  complex  transmission  and  reflection  coefficients.  With  6 
small,  we  expand  o,  T  and  R  as  series  in  the  small  parameter  e  =  D/h.  and 
obtain 

do=c’'^  Ro  =  0  To  =  l  (12) 

at  leading  order.  At  second  order,  the  reflection  coefficient  Ri  for  an  arbitrary 
topography  ^(x)  is 


(13) 


as  found  by  Miles  (1981 ).  Note  that  Ri  is  singular  in  the  limit  ^ls  0  — ♦  ir/2.  This 
effect  has  not  been  previously  noted  and  its  practical  implications  are  unclear. 
For  the  case  of  a  simple  sinusoidal  bottom 


6{x) 

we  obtain  the  expression 


D  sin  (Ax)  ;  0  <  X  <  Nt,L 


1  .  U-\ 

\  J  p  2  ’  X  ~  * 


(14) 


(15) 


This  result  extends  the  non-resonant  theory  of  Davies  and  Heathershaw  (1984) 
to  include  obliquely  incident  waves. 
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For  the  periodic  bar  field  described  by  (5)  and  (6),  we  substitute  (5)  in  (13) 
and  obtain  the  expression 


(16) 


where  the  integrals  /„  are  given  by 


~  J-k 


cos(nAx)e*’^*<ix- 


(17) 


As  in  the  case  of  a  single  sinusoidal  bar,  the  integral  In  takes  on  special  values 
when  2ifnX  =  1  for  the  corresponding  value  of  n.  We  further  simplify  the 
notation  by  setting 

2/ 

7  =  y  (18) 

Then,  for  7  ^  n,  we  obtain  the  expression 


(19) 


For  n  =  7,  we  obtain  the  expression 


/n(7  =  n)  = 


N^L 


(20) 


We  thus  obtain  the  general  solution  for  obliquely  incident  waves 
^  lo  (C'-rn}\  f  ^  7^£>„  j  p.  r, 

(21) 

where  6(n  —  7)  is  the  delta  function,  and  there  n  =  7  can  only  occur  for  one 
wavenumber  component  for  a  fixed  value  of  1.  For  the  case  of  normally  incident 
waves  (studied  further  below),  we  let  f  — ♦  fc,  m  —»  0  and  obtain 


kNi,L\ 


«.  =  -f  {£  / 


(22) 


From  the  form  of  the  solution,  it  is  apparent  that  each  harmonic  of  the  bar  field 
contributes  to  the  reflection  process,  with  the  dominant  contribution  of  the  n*** 
harmonic  coming  from  the  neighborhood  \  fan.  An  example  plot  of  reflection 
coefficient  |  i?i  (  is  shown  in  Figure  2  for  a  case  of  4  bars  with  crest- to-crest 
spacing  equal  to  the  unrectified  wavelength  (6i,  =  f«/2).  Waves  are  normally 
incident  on  the  bar  field.  The  peak  in  |  |  at  2A:/A  =  1  corresponds  to  the'u^uaT 

Bragg  interaction  between  the  surface  wave  and  the  fundamental  harmonic  of 
the  oar  field,  when  the  surface  wave  length  is  twice  the  bar  spacing.  A  second 
prominent  peak  is  located  at  24:/A  =  2,  corresponding  to  a  surface  wavelength 
equal  to  twice  the  length  of  the  second  harmonic  of  the  bar  field  (and  thus  equal 
in  length  to  the  bar  spacing).  This  strong  second  peak  is  absent  when  the  bar 
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Figure  2:  Reflection  for  four  discretely  placed  bars.  Non-resonant  theory,  equa¬ 
tion  (22). 


field  being  considered  is  a  simple  sinusoid,  as  in  Davies  and  Heathershaw  (19S4) 

and  Mei  (1985).  in-  m  ■ 

In  general,  the  relative  amplitude  of  the  peaks  in  the  reflection  coefficient 
may  be  adjusted  by  changing  the  spacing  of  artificial  bars,  assuming  the  cross- 
section  of  each  bar  to  stay  the  same.  Pushing  bars  closer  together  makes  the  bar 
field  more  sinusoidal  and  reduces  the  importance  of  higher  harnronics.  Placing 
the  bars  further  apart  makes  them  into  relatively  more  “solitary”  features,  and 
thus  emphasizes  the  relative  importance  of  higher  harmonics.  Two  cases  illus¬ 
trating  these  extremes  were  investigated  in  the  experiments  described  below. 
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Resonant  Reflection 


The  reflection  coefficient  described  in  (21)  is  defective  for  the  cases  2t.  w  nA, 
where  the  coefficient  can  become  arbitrarily  large  as  — »  oo.  While  this 

limit  would  never  be  reached  in  practice,  the  result  shows  that  the  theory  is 
not  strictly  valid  in  the  neighborhood  of  the  resonances.  The  problem  lies 
with  assuming  that  is  '  '^(1)  in  the  perturbation  series  used  above.  Mei 
(1985)  hais  developed  a  resoni..'.ce  theory  which  allows  for  0(1)  reflection  in  the 
neighborhood  of  each  resonance.  Mei  et  al.  (1988)  further  suggested  that,  for 
the  case  of  a  bar  field  with  multiple  Fourier  components,  the  reflection  could  be 
estimated  using  the  resonance  theory  applied  to  the  Fourier  mode  corresponding 
to  the  bar  wavelength.  This  approach  would  not  account  for  the  occurrence  of 
multiple  strong  peaks.  In  the  present  study,  we  define  a  neighborhood  of  each 
resonance  21  fX  =  n  to  be  the  range  n  —  1/2  <  2£/A  <  n  +  1/2.  Then,  in  each 
range,  Mei’s  theory  is  used  with  2^/A*  replacing  2^/A,  with  A*  =  nA.  We  refer 
the  reader  to  Mei  (1985)  for  the  expressions  defining  the  reflection  coefficients. 
The  only  necessary  modification  to  the  theory  account  for  oblique  incidence  and 
the  presence  of  multiple  resonant  peaks.  The  frequency  of  the  n‘^  resonant 
peak  is  given  by 

gn\  .  (  nXh 


Uil  = 


2cosd 


tanh 


/  nXh  N 
\  2  cos  0  ) 


The  cutoff  condition  fl©  defined  by  Mei  is  replaced  by 


Hon  = 


0  —  m} 

P  2sinh2jth 


k  = 


cosB 


(23) 


(24) 


where  the  are  the  amplitudes  of  the  bar  Fourier  coefficients,  and  there  flon 
refers  to  the  resonant  peak. 

An  example  of  the  reflection  calculated  for  the  case  of  normal  incidence  is 
given  in  the  following  section,  in  comparison  with  numerical  results  and  results 
of  the  non-resonant  theory. 

Nuinsrical  Solutions 

In  order  to  study  the  validity  of  each  of  the  perturbation  solutions,  direct 
numerical  solutions  of  equation  (10)  were  also  performed.  For  a  bar  field  in  the 
region  0  <  i  <  an  incident  wave  boundary  is  established  at  i  =  A  <  0, 
and  a  downwave,  transmitting  boundary  is  established  hi  x  —  B  >  N^L.  For 
an  incident  wave  the  appropriate  boundary  conditions  are 


=  I 


i({2^i  -  ^)  ;  X  =  A 
;  X  =  B 


(25) 


Equations  (10)  and  (25)  are  finite-differenced  using  central  differences,  leading 
to  a  tridiagonal  system  which  is  solved  using  the  Thomas  algorithm. 

Figure  3  shows  a  sample  of  calculated  reflection  coefficients  obtained  with 
the  numerical  solution  and  the  two  analytic  solutions,  for  the  bar  field  described 
in  Figure  2.  As  expected,  the  non-resonant  solution  over-predicts  reflection  at 
2k fX  =  1,  in  comparison  with  the  resonant  theory  of  Mei  (1985).  The  discrep¬ 
ancy  is  relatively  minor  at  the  second  peak  2k/ X  —  2,  where  the  resonance  is 
relatively  weaker.  In  contrast  to  both  analytic  theories,  the  numerical  results 
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Figure  3:  Comparison  of  non-resonant,  resonant  and  numerical  solutions  for  a 

bar  field  with  N(,  =  4.  . .  ,  non-resonant  theory; - ,  resonant  theory: 

—  -  — ,  numerical  results. 

show  a  strong  downshift  of  the  reflection  peaks  to  lower  values  of  2k/ X.  This 
downshift  is  related  to  the  simultaneous  interaction  between  the  wave  field  and 
several  bottom  modes;  a  similar  effect  occurs  for  the  case  of  a  sinusoidal  bar 
field,  but  it  is  much  more  subtle.  The  numerical  scheme  also  predicts  a  higher 
reflection  coefficient  at  each  peak.  The  large  downshifts  and  higher  peaks  are 
largely  validated  by  data  described  below.  These  results  indicate  that  either 
of  the  two  analytic  solutions  are  at  best  qualitatively  accurate  when  used  to 
describe  reflection  from  the  type  of  bars  that  could  be  built  in  an  actual  con¬ 
struction  project. 

Experimental  Results 

Experiments  were  conducted  in  the  60  cm  wide  wave  flume  in  the  Coastal 
and  Oceanographic  Engineering  Laboratory,  University  of  Florida,  in  order  to 
verify  the  basic  aspects  of  the  theory  for  normally  incident  waves.  For  the 
experiments,  a  water  depth  h  =  lb  cm  was  used.  Bar  height  D  was  5  cm, 
giving  D/h  =  0.33,  which  is  relatively  large  and  could  contribute  to  some  of 
the  discrepancies  between  theory  and  data  noted  below.  The  bar  footprint 
bi  =  50cm.  Two  bar  spacings,  L  =  80cm  and  120cm,  were  tested,  corresponding 
to  cases  with  bar  field  higher  harmonics  of  low  importance  and  great  importance. 
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INCIDENT  ENERGY  SPECTRUM 

1  REFLECTED  ENERGY  SPECTRUM  TRANSMITTED  ENBIGY  SPECTRUM 
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Figure  4:  Bar  field  and  wave  gage  placement 


respectively.  The  bar  fields  tested  contained  four  discrete  bars. 

In  order  to  maintain  a  close  correspondence  between  the  assumed  linearity 
of  the  wave  theory  and  the  experiments,  incident  waves  on  the  order  of  1  cm  in 
height  were  generated.  (Actual  height  varied  with  wave  frequency,  as  would  be 
expected  from  wavemaker  theory.)  Wave  heights  were  measured  using  capac¬ 
itance  wave  gages  mounted  with  6  cm-long  wires,  which  were  calibrated  over 
their  full  length.  Data  was  sampled  using  12-bit  digitization,  giving  a  resolution 
of  0.012  mm/division.  The  wavefield  was  sampled  at  a  lO/fz  frequency,  with 
experimental  waves  being  generated  in  the  range  0.4 /fz  <  /  <  l.6Hz. 

Due  to  the  small  wave  heights  being  used,  there  was  an  additional  source  of 
noise  in  the  data  associated  with  mechanical  vibration  in  the  wavemaker  and 
other  high-frequency  effects.  In  retrospect,  it  would  be  better  to  use  slightly 
higher  waves  in  future  experiments  unless  great  care  were  taken  to  isolate  me¬ 
chanical  vibration.  (For  a  particularly  spectacular  example  of  clean  data  in  a 
related  low-amplitude  wave  experiment,  see  Benjamin  et  al.  (1989)). 

The  three-gage,  least  squares  method  developed  by  Funke  and  Mansard 
(1980)  was  used  to  separate  incident  and  reflected  waves.  The  gage  layout 
relative  to  the  bar  field  is  shown  in  Figure  4.  The  incident-reflected  separation 
was  performed  both  upwave  and  downwave  of  the  experimental  bar  field.  The 
downwave  separation  indicated  a  reflection  from  the  absorbing  beach  on  the 
order  of  5-6%.  This  reflected  energy  was  neglected  in  subsequent  processing 
and  the  downwave  region  was  assumed  to  be  perfectly  transmitting. 

Figure  5  shows  the  measured  reflection  coefficient  for  the  case  of  L  =  SOcm, 
when  bar  field  harmonics  are  relatively  unimportant.  Also  included  in  the  fi^re 
is  the  prediction  of  the  numerical  model  described  in  the  previous  section.  The 
data  largely  validates  the  theory,  although  there  is  a  great  deal  of  scatter.  (It 
also  appears  that  shifting  the  data  to  higher  values  of  2k /X  would  bring  the  data 
into  fairly  close  agreement  with  theory.  No  systematic  error  was  ever  detected  in 
the  experimental  procedure  whidi  could  account  for  such  a  shift,  unfortunately.) 

Figure  6  shows  corresponding  data  for  the  case  of  L  —  120cm,  where  the 
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Figure  6;  Experimentally  measured  reflection  coefficient,  L  =  120cm. 

second  harmonic  component  of  the  bar  field  is  comparable  in  height  to  the  fun¬ 
damental  harmonic.  Due  to  time  constraints,  the  data  here  is  relatively  sparse, 
and  tests  were  grouped  in  order  to  show  the  relative  heights  of  the  two  reflection 
peaks.  Both  the  presence  and  the  relative  impoitance  of  the  two  numerically- 
predicted  peaks  are  substantiated  by  the  data. 

Response  of  a  Closed-end  Channel 

One  question  that  arises  in  response  to  the  realization  that  bars  can  reflect 
significant  amounts  of  incident  energy  is  whether  or  not  a  region  downwave  of  a 
bar  field  experiences  a  less  severe  wave  condition  than  the  region  on  the  incident 
site.  The  answer  to  this  question  can  be  positive  or  negative,  depending  on  the 
geometry  of  the  downwave  region  and  the  reflectivity  of  the  end  boundary.  For 
cases  where  reflection  from  the  end  wall  is  nearly  complete,  waves  travelling 
back  towards  the  bar  field  are  partially  re-reflected  into  the  sheltered  region. 
The  possibility  of  resonating  the  sheltered  region  exists,  as  does  the  possibility 
of  reducing  the  wave  activity,  and  depends  primarily  on  whether  the  sheltered 
region  contains  an  integer  multiple  of  one-half  the  surface  wavelen^h. 

Figure  7  shows  the  numerically  predicted  amplitude  at  a  vertical  wall  sit¬ 
uated  four  barfield  wavelengths  downwave  of  a  bar  field  with  four  bars,  as  in 
the  previous  examples.  The  incident  wave  has  an  amplitude  of  unity,  and  so  an 
amplitude  of  2  represents  simple  reflection.  The  figure  shows  that  the  amplitude 
at  the  wall  can  reach  as  high  as  3.6  and  as  low  as  1,  representing  a  range  of 
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Figure  7:  Wave  Amplitude  at  End  Wall  of  a  Channel  Sheltered  by  an  Artificial 
Bar  Field 


resonance  and  sheltering  conditions. 

For  the  case  where  the  channel  end  is  primarily  absorbing  (due  to  wave 
breaking  or  frictional  effects),  the  possibility  of  resonating  the  sheltered  region 
is  greatly  reduced.  This  virtually  guarantees  that  a  bar  field  designed  to  shelter 
a  beach  from  wind- wave  band  waves  would  not  resonate  the  shoreline.  How¬ 
ever,  as  the  wave  frequency  becomes  low,  even  a  mildly-sloped  shoreline  can 
become  essentially  reflective.  It  is  thus  possible  that  broad,  low  bars  contribute 
significantly  to  amplifying  long  wave  energy  on  the  beach  face.  The  long  w’aves 
that  could  be  amplifi^  or  resonated  by  this  mechanism  may  be  locally  gener¬ 
ated  by  nonlinear  processes  in  the  surfzone,  or  they  may  be  arriving  as  part  of 
the  forced  or  free  long  wave  climate  incident  from  offshore.  Problems  of  this 
sort  need  further  investigation  to  determine  the  relative  importance  of  bottom 
interaction  in  influencing  the  nearshore  wave  climate. 
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ABSTRACT 

The  feasibility  of  a  new  type  of  breakwater  is  explored 
through  a  combined  prograim  of  theorteical  analysis,  lab* 
oratory  experiments  and  a  prototype  field  test.  The 
breakwater  consists  of  a  series  of  low  height,  shore- 
parallel  bars  which  are  placed  just  outside  the  surf- 
zone.  The  incident  wave  field  is  strongly  reflected 
when  the  incident  wave  length  is  equal  to  twice  the  dis* 
tance  between  adjacent  bars.  The  breakwater  acts  to 
shelter  the  beach  against  storm  wave  attack  and  to  build 
a  tombolo  sand  deposit  behind  the  breakwater.  Theoreti¬ 
cal  and  laboratory  studies  indicated  that  the  Bragg  re¬ 
flection  breakwater  concept  may  have  considerable  merit, 
however,  a  small  scale  field  test  served  to  demonstrate 
the  many  practical  difficulties  in  implementing  the 
concept. 

INTRODUCTION 

The  Naval  Civil  Engineering  Laboratory  (NCEL)  has  been 
involved  in  a  multi-institutional  program  to  explore  low 
cost,  rapidly  deployable  techniques  for  reducing  beach 
erosion  during  storms.  One  concept  under  consideration 
is  the  Bragg  reflection  breakwater.  Theoretical  and 
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laboratory  studies  have  shown  that  a  series  of  periodi¬ 
cally  spaced,  low  height  bars,  oriented  parallel  to  the 
beach,  can  produce  a  strong  reflection  of  the  incident 
wave  field.  A  reasonant  condition  (termed  Bragg  reflec¬ 
tion)  occurs  when  the  incident  wave  length  is  equal  to 
twice  the  separation  distance  between  bars. 

The  Bragg  reflection  breaJcwater  concept  involves  placing 
a  series  of  artificial  bars  just  outside  the  surfzone. 
Referring  to  Figure  1,  the  effect  of  the  bar  field  is 
two- fold:  to  shelter  the  beach  from  storm  wave  attack; 
and,  to  create  a  large  sand  volume  behind  the  breakwater 
by  way  of  the  induced  nearshore  circulation.  Both  ef¬ 
fects  serve  to  protect  the  beach  against  storm  wave  at¬ 
tack:  the  first  by  reducing  the  rate  of  offshore  sand 
transport,  the  second  by  increasing  the  erodable  sand 
volume . 

The  objective  of  the  research  program  was  to  explore  the 
feasibility  of  the  Bragg  reflection  breakwater  concept 
through  a  combined  program  of  numerical  and  laboratory 
modeling  and  a  prototype  field  test.  The  present  paper 
presents  an  overview  of  the  program  with  an  emphasis  on 
the  planning  and  results  of  a  small  scale  field 
experiment. 


BACKGROUND 

Theoretical  and  Experimental  Basis 

Davies  and  Heathershaw  (1984),  Mei  (1985)  and  others 
have  shown  that  a  series  of  low  amplitude  sinusoidal  un¬ 
dulations  on  the  seabed  can  be  an  effective  reflector  of 
surface  gravity  waves.  The  reflected  wave  energy  varies 
as  a  function  of  the  ratio  of  the  incident  wave  length 
to  the  spacing  between  bottom  undulations.  A  reasonant 
condition,  termed  Bragg  reflection,  occurs  when  the  in¬ 
cident  wave  length  is  equal  to  twice  the  distance  be¬ 
tween  undulations.  When  resonance  occurs,  a  small  num¬ 
ber  of  bars  can  reflect  a  substantial  portion  of  the 
incident  wave  energy. 

Technical  Issues 

At  the  onset  of  the  research  program,  there  were  a  num¬ 
ber  of  unresolved  technical  issues  relating  to  the  fea¬ 
sibility  of  the  Bragg  reflection  breakwater  concept. 
These  included:  the  response  characteristics  of  non-si- 
nusoidal  bars  on  a  sloping  beach;  the  effects  of  finite 
bar  length  and  the  resulting  wave-induced  circulation; 
the  morphological  response  of  the  beach  in  the  presence 
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Figure  1.  Schematic  drawing  of  Bragg  reflection 
breakwater  concept. 
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Figure  2.  Beach  profile  response  model  results. 
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of  the  breakwater;  the  practical  design,  installation 
and  recovery  of  a  brealcwater;  and,  the  effects  of  waves 
on  bar  stability  and  scour.  An  additional  issue  was  how 
to  neasure  spectral  wave  reflection  in  the  field. 

PRELIMINARY  ASSESSMENT 

As  a  means  of  exploring  the  feasibility  of  the  Bragg  re¬ 
flection  brealcwater  concept,  a  series  of  preliminary 
studies  were  conducted.  These  studies  were  to  culminate 
in  a  small  scale  field  test  of  the  concept.  The  results 
of  the  preliminary  studies  can  be  summarized  as  follows. 

Wave  and  Current  Modeling 

Numerical  and  laboratory  modeling  studies  were  conducted 
at  the  University  of  Florida  to  examine  the  effects  of 
bar  shape,  bar  placement,  beach  slope  and  bar  length  on 
the  reflection  characteristics  of  a  Bragg  reflection 
breakwater.  Details  of  these  studies  may  be  found  in 
Anton  et  al.  (1990).  Specific  results  included: 

o  The  primary  effect  of  beach  slope  is  to  require  an 
adjustment  in  the  spacing  between  each  bars  to 
accommodate  the  change  in  incident  wave  length. 

o  The  reflection  from  a  non-sinusoidal  bar  field  can 
be  calculated  by  expanding  the  bottom  shape  func¬ 
tion  as  a  fourier  series  and  linearly  summing  the 
response  characteristics  for  each  component  of  the 
series. 

o  Finite  bar  length  causes  a  longshore  variation  in 
the  wave-induced  setup  leading  to  the  formation  of 
a  rip  current  flowing  outward  over  the  bar  field. 

The  rip  current  acts  to  broaden  the  reasonant  peak. 

These  studies  indicated  that  a  Bragg  reflection  breakwa¬ 
ter,  constructed  from  a  small  number  of  practical -shaped 
(i.e.  non-sinusoidal)  bars,  could  be  expected  to  gener¬ 
ate  a  reflection  coefficient  in  the  range  of  0.2  to  0.4. 

Beach  Response  Modeling 

A  Bragg  reflection  breakwater  protects  a  beach  from 
storm  wave  attack  by  reducing  the  rate  of  offshore  sand 
transport  via  decreased  wave  height  and  increasing  the 
erodable  sand  volvuae  via  creation  of  a  tombolo  behind 
the  brealcwater.  The  former  is  a  two-dimensional  effect 
which  is  independent  of  breakwater  length.  The  latter 
is  a  three-dimensional  effect  which  results  from  the 
formation  of  a  nearshore  circulation  cell  with  a  rip 
current  flowing  outward  over  the  top  of  the  breakwater. 
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NCEL  conducted  a  numerical  model  study  to  detexmine  the 
impact  of  a  Bragg  reflection  brealcvater  on  the  two-di¬ 
mensional  response  of  a  beach  profile  during  a  simulated 
storm.  The  model  study  utilized  an  energetics-based 
beach  profile  response  model  developed  by  DeVries  and 
Bailard  (1988).  The  model  was  nin  for  a  simulated  5  day 
storm  assuming  varying  wave  reflection  coefficients. 

The  effectiveness  of  the  breedcvater  was  expressed  in 
terms  of  the  normalized  erosion  volume.  This  was  de¬ 
fined  as  the  ratio  of  the  eroded  sand  volume  in  the 
presence  of  the  brea]cwater  divided  by  the  eroded  sand 
volume  without  the  breakwater. 

Model  inputs  were  as  follows: 

o  Initial  equilibrium  beach  slope  ■  0.03 
o  Grain  size  «  0.4  mm 
o  Storm  duration  «  72  hr 
o  Peak  wave  height  ^  2  m 
o  Peak  wave  period  »  9  sec 

Without  the  breakwater,  the  beach  eroded  approximately 
205  cubic  meters  per  meter  of  beach.  Referring  to  Fig¬ 
ure  2,  the  breakwater  reduced  the  storm  erosion  volume 
by  an  amount  which  was  inversely  proportional  to  the 
wave  reflection  coefficient,  K^.  The  degree  of  reduc¬ 
tion  is  quite  significant,  amounting  to  35%  for  a  re¬ 
flection  coefficient  of  0.2. 

A  movable  bed  physical  model  study  was  conducted  at  the 
University  of  Florida  to  qualitatively  examine  the 
three-dimensional  beach  changes  induced  by  a  Bragg  re¬ 
flection  breakwater.  The  model  confirmed  the  presence 
of  a  nearshore  circulation  cell  with  a  rip  current  flow¬ 
ing  outward  over  the  bar  field.  The  circulation  cell 
caused  a  tombolo  to  form  behind  the  breakwater,  with  the 
outer  edge  of  the  tombolo  perched  on  the  shoreward-most 
bar. 


Bar  Module  Tests 

It  was  anticipated  that  developing  a  practical  breakwa¬ 
ter  would  be  a  significant  design  challenge.  Our  lim¬ 
ited  budget  necessitated  coming  up  with  a  breakwater 
which  could  be  rapidly  assembled,  installed  and  recov¬ 
ered  using  a  minimum  of  people  and  equipment.  In  order 
to  explore  various  design  and  installation  concepts, 

NCEL  conducted  field  tests  of  two  bar  designs  at  a  Port 
Hueneme  beach.  Both  designs  featured  modular  construc¬ 
tions  consisting  of  skid-mounted  bar  modules  which  could 
be  assembled  on  the  beach  and  dragged  offshore. 
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Referring  to  Figure  3,  the  first  design  consisted  of  a 
geotextile  bag  attached  to  a  skid-mounted  steel  frame. 
The  concept  was  to  drag  the  modules  into  place  with  the 
bags  empty,  and  then  fill  the  bags  with  sand  using  a 
small  dredge  pump.  When  the  breakwater  was  no  longer 
needed,  the  bags  could  be  slit,  allowing  the  sand  to 
disperse  and  the  bars  dragged  back  ashore.  The  advan¬ 
tage  of  this  design  was  that  the  bar  modules  would  be 
light  and  easy  to  move  when  empty,  but  heavy  and  diffi¬ 
cult  to  move  when  full. 

An  attempted  field  test  of  the  geotextile  bar  design 
ended  in  failure.  Although  assembling  and  positioning 
the  bar  nodule  proved  simple,  filling  the  geotextile  bag 
with  a  small  dredge  pump  proved  difficult.  The  pump 
was  mounted  on  an  amphibious  LARC  vehicle,  parked  in  the 
surfzone.  The  principal  difficulties  were  maintaining 
an  adequate  supply  of  sand  to  the  pump  and  handling  the 
intake  and  discharge  hoses.  VThen,  after  a  few  hours  of 
filling,  little  sand  had  been  pumped  into  the  bag,  the 
test  was  abandoned.  In  retrospect,  the  test  might  have 
been  more  successful  had  a  larger,  perhaps  land-based 
dredge  pump  had  been  used. 

Nevertheless,  the  installation  procedure  was  judged  to 
be  too  cumbersome  for  the  planned  field  experiment. 

Referring  to  Figure  4,  the  second  bar  module  design  con¬ 
sisted  of  a  corrugated  steel  arch  attached  to  a  skid- 
mounted  steel  frame.  The  concept  was  to  assemble  the 
modules  on  the  beach  and  drag  then  into  position.  When 
the  breakwater  was  no  longer  needed,  the  modules  would 
be  dragged  back  onto  the  beach  and  removed.  The  advan¬ 
tage  of  this  design  was  that  the  modules  were  ready  to 
go  once  they  were  moved  into  position.  The  main  draw¬ 
back  was  that  the  modules  were  relatively  light  weight 
and  could  be  moved  about  by  large  waves.  Although  pin¬ 
ing  the  modules  with  uplift  resisting  anchors  was  con¬ 
sidered,  it  was  judged  unnecessary  for  the  anticipated 
wave  climate. 

The  field  test  demonstrated  that  the  arch  module  design 
could  be  rapidly  assembled  and  installed.  Provided  wave 
heights  were  less  than  1  m  (rms) ,  the  module  tended  to 
remain  in  place.  The  primary  problem  was  a  scouring  of 
the  seabed  underneath  the  bar  module.  This  %ras  caused 
by  wave- induced  flow  passing  through  the  narrow  gap  be¬ 
tween  the  leading  and  trailing  edges  of  the  bar  and  the 
sand  bottom.  The  scour  depression  acted  to  further  en¬ 
hance  the  venting  flow,  reducing  the  reflectivity  of  the 
bar  module. 
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Figure  3.  Geotextile  bar  module  design. 
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Figure  4.  Steel  arch  bar  module  design  i 
field  test. 
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The  venting/scour  problem  was  solved  by  adding  steel 
plates  to  the  leading  and  trailing  edges  of  the  bar. 

The  plates  were  attached  to  the  bar  by  hinges,  effec¬ 
tively  sealing  off  the  gap  vindemeath  the  bar.  Although 
sooe  scour  continued  to  occur  around  the  bar  module,  the 
plates  continued  to  function  by  dropping  down  into  the 
scour  depressions.  After  a  period  of  two  days,  the 
scour  depressions  appeared  to  stabilize  at  a  depth  of 
about  0.3  m.  A  few  days  later,  a  significant  storm 
passed  through  the  area  destroying  the  bar  modules. 
However,  the  successful  performance  of  the  module  prior 
to  the  storm  convinced  us  to  that  the  steel  arch  design 
would  be  suitable  for  the  planned  field  test. 

Detection  Method 

An  inverse  method  was  developed  to  estimate  the  direc¬ 
tional  spectrum  of  the  reflected  wave  field  (see  Herbers 
and  Guza,  1990) .  The  method  requires  simultaneous  mea¬ 
surement  of  the  incident  wave  field  at  a  point  offshore 
and  the  reflected  wave  field  at  a  point  immediately  in 
front  of  the  Bragg  reflection  brealcwater.  The  procedure 
is  as  follows: 

0  The  incident  wave  field  is  estimated  using  the 
offshore  wave  measurements. 

o  The  incident  wave  field  is  transformed  to  shallow 
water  using  linear  refraction. 

o  The  reflected  wave  field  is  estimated  using  the 
inshore  wave  measurements  and  subject  to  maximum 
compatibility  with  the  refracted  incident  wave 
field. 

PROTOTYPE  FIELD  TEST 
Site  Selection 

In  order  to  maximize  the  detectable  reflected  wave  en¬ 
ergy,  we  wanted  to  construct  a  breakwater  that  was  sig¬ 
nificantly  longer  than  the  incident  wave  length  and 
which  produced  a  reflection  coefficient  of  about  0.4. 

To  minimize  the  size  of  the  brealcwater  and  therefor  its 
cost,  we  decided  to  carry  out  our  experiment  in  "super 
laboratory"  conditions.  These  can  be  defined  as  signif¬ 
icant  wave  heights  less  than  0.25  m,  a  wave  period  of 
about  5  seconds  and  an  installation  depth  of  about  1.5 
meters. 

A  number  of  sites  were  considered  for  the  field  experi¬ 
ment.  The  desired  site  required  small  amplitude,  short 
period  waves,  a  long  planar  beach,  limited  public  access 
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and  good  logistical  support.  Based  on  these  criteria, 
Cape  Canaveral  Beach,  Florida,  was  selected  as  the  pre¬ 
ferred  site  for  the  field  experiment.  The  site  was  lo¬ 
cated  within  the  confines  of  Cape  Canaveral  Air  Force 
Base,  limiting  public  access  and  providing  the  necessary 
l^istical  support.  Two  years  of  wave  data  from  the 
site  indicated  that  during  the  month  of  July,  we  could 
expect  significant  wave  heights  of  about  0.2  m  with  a  5 
second  period. 

Breakwater  Design 

Referring  to  Figure  5,  the  prototype  Bragg  reflection 
breakwater  was  composed  of  three  bars,  90  meters  long. 
The  offshore  and  middle  bars  were  0.7  meters  high  and 
2.2  meters  wide.  The  inner  bar  was  0.6  meters  high  and 
2.4  meters  wide.  The  spacing  between  the  outer  and  mid¬ 
dle  bar  was  9  meters  versus  8  meters  between  the  middle 
and  inner  bars.  The  breakwater  was  designed  to  be  in¬ 
stalled  in  a  water  depth  of  about  0.76  meters  MLW. 

Referring  to  Figure  5,  each  bar  module  was  7.3  meters 
long.  The  modules  were  fabricated  from  a  pair  of  corru¬ 
gated  steel  arches  attached  to  a  steel  skid-mounted 
frame.  The  frame  was  composed  of  a  pair  of  steel  angle 
members  attached  to  two  steel  skid  beams.  The  angle 
members  served  as  attachment  points  for  the  edges  of  the 
corrugated  steel  arches  and  the  a  series  of  steel  scour 
plates  hinged  to  the  front  and  rear  edges  of  the  bar. 

The  bar  modules  were  designed  to  be  assembled  by  a  small 
crew  using  a  crane  and  hand  tools.  Assembly  was  facili¬ 
tated  by  prefabricating  the  components  off-site  and  us¬ 
ing  nuts  and  bolts  for  fasteners. 

The  total  weight  of  the  bar  modules  was  approximately 
3300  lb.  An  stability  analysis  indicated  that  the  mod¬ 
ules  would  remain  in  place  without  anchoring  if  rms  wave 
heights  were  less  than  1  m  (assuming  a  wave  period  of  5 
to  6  sec) . 

Anticipated  Performance 

Figure  6  shows  a  plot  of  the  estimated  reflection  char¬ 
acteristics  for  the  prototype  breakwater.  These  charac¬ 
teristics  were  estimated  using  a  computer  program  devel¬ 
oped  by  Kirby  (1987).  The  response  function  shows  a 
broad  peak  centered  at  about  6  seconds  with  a  reflection 
coefficient  value  of  0.4.  The  adjacent  narrower  peak  is 
the  first  harmonic  resulting  from  the  non-sinusoidal 
shape  of  the  bar  field  (i.e.  discrete  bars  resting  on  a 
planar  bottom) . 
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Installation  and  Recovery 


The  bar  modules  were  assembled  in  a  staging  area  located 
about  200  meters  from  the  beach.  After  assembly,  the 
modules  were  carried  to  the  beach  with  a  large  front-end 
loader.  The  loader  placed  the  modules  on  the  upper 
beach  face  in  groups  of  three,  with  their  skids  oriented 
perpendicular  to  the  shore.  Scour  plates  were  attached 
to  each  module  and  tied  up  out  of  the  way  to  facilitate 
nodule  dragging. 

Referring  to  Figure  7,  the  original  installation  plan 
called  for  dragging  the  modules  into  position  in  groups 
of  three  using  a  shore-based  winch  and  a  movable  sheave 
and  anchor  assembly.  In  practice,  we  found  that  the 
waves  at  the  site  were  too  large  to  allow  this  method  of 
installation.  In  particular,  the  LARC  vehicle  found  it 
difficult  to  negotiate  the  surfzone  while  pulling  a 
heavy  anchor  and  cable.  As  a  result,  only  one  group  of 
nodules  was  installed  using  this  technique.  The  rest  of 
the  modules  were  installed  individually  at  extreme  low 
tide  using  the  front-end  loader.  Although  rapid  (all  of 
the  modules  were  installed  in  a  single  low  tide) ,  this 
alternative  technique  required  placing  the  breakwater 
approximately  15  meters  closer  to  shore.  As  a  result, 
the  breakwater  was  situated  inside  the  surfzone  for  a 
substantial  portion  of  the  tidal  cycle. 

Monitoring  Plan 

A  survey  grid  was  established  in  the  area  of  the  break¬ 
water  and  in  an  adjacent  control  area.  Daily  wading 
profiles  were  to  be  conducted  in  the  test  and  control 
areas.  These  were  to  be  supplemented  by  combined  wading 
and  fathometer  surveys  at  the  start  and  end  of  the  ex¬ 
periment. 

The  incident  directional  wave  spectrum  was  measured  at  a 
depth  of  6  meters  using  a  linear  array  of  pressure  sen¬ 
sors.  A  second  array  of  pressure  sensors  and  current 
meters  was  to  measure  the  reflected  wave  field  and  the 
induced  nearshore  circulation  field  around  the  breakwa¬ 
ter.  The  offshore  array  was  constructed  prior  to  the 
installation  of  the  breakwater.  The  inshore  array  was 
to  be  constructed  once  the  breakwater  was  in  place.  Un¬ 
fortunately,  this  never  occurred  due  to  problems  de¬ 
scribed  below. 

RESULTS 


The  Bragg  reflection  breakwater  was  deployed  on  11  July 
1988.  During  the  week  prior  to  installation,  the  sig¬ 
nificant  wave  height  at  the  site  averaged  about  0.4  m. 
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Figure  7.  Planned  installation  method. 


Figure  8.  Observed  scour  pattern  in  field  test. 
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This  was  approximately  twice  the  anticipated  wave  height 
based  on  two  years  of  measured  wave  data  and  essentially 
equal  to  the  maximum  design  wave  height  for  bar  stabil¬ 
ity. 

Immediately  following  brealcwater  installation,  the  sig¬ 
nificant  wave  height  increased  to  2d>out  0.6  m  and 
remained  at  this  level  for  the  next  36  hours.  These 
waves  were  significantly  greater  than  the  maximum  design 
wave,  causing  the  bar  modules  to  slowly  shift  their 
positions.  The  modules  rotated  towards  the  direction  of 
wave  attack  and  moved  shoreward  about  1  to  2  m.  The 
movement  began  with  the  outer-most  bar  and  progressed 
inward  through  the  middle  and  inner  bars. 

Another  problem  which  developed  almost  immediately  after 
installation  of  the  brea)cwater  was  an  intense  scouring 
around  the  bar  modules.  As  the  tide  rose,  the  bars  were 
observed  to  be  generating  a  significant  degree  of  tur¬ 
bidity.  On  the  following  day  we  discovered  large  scour 
holes  had  formed  in  front  of  and  behind  the  bars  (see 
Figure  8) .  The  scour  holes  slowly  grew,  becoming  deeper 
and  wider.  Eventually  they  reached  more  than  2  meters 
deep,  threatening  the  stability  of  the  bar  modules. 

The  observed  scour  was  approximately  an  order  of  magni¬ 
tude  greater  than  the  scour  observed  in  the  Port  Hueneme 
tests.  Since  wave  conditions  were  similar,  we  believe 
that  the  increased  scour  was  attributable  to  the  differ¬ 
ent  properties  of  the  two  beach  sands.  At  Port  Hueneme, 
the  sand  was  well  sorted,  having  a  median  diameter  of 
0.25  mm  and  quartz  and  feldspar  composition.  At  Cape 
Canaveral,  the  sand  was  poorly  sorted,  having  a  medium 
diameter  of  0.14  mm  and  calcium  carbonate  composition. 
Apparently  the  lighter  specific  gravity  of  the  Cape 
Canaveral  sand,  coupled  with  its  smaller  size  and  higher 
percentage  of  fines,  resulted  in  increased  sand  suspen¬ 
sion  and  enhanced  scour. 

After  two  days,  the  bar  modules  had  become  badly  scat¬ 
tered.  Some  of  the  modules  had  shifted  into  their  scour 
holes  and  begun  to  become  buried.  With  bar  burial  be¬ 
coming  a  growing  problem,  the  decision  was  made  to  ter¬ 
minate  the  experiment  and  recover  the  bar  modules  while 
it  was  still  possible.  A  bulldozer  was  used  to  pull  the 
individual  modules  out  of  the  surf  zone  and  back  up  the 
beach  face.  From  there,  a  front-end  loader  transported 
the  modules  back  to  the  staging  area  where  they  were 
disassembled. 

Because  of  the  short  duration  of  the  breakwater  deploy¬ 
ment,  the  nearshore  sensor  array  was  never  installed. 

As  a  result,  no  quantitative  measurements  were  obtained. 
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visual  observations  during  the  first  few  hours  of  de~ 
ployaent  indicated  that  some  wave  reflection  was  occur** 
ring  as  evidenced  by  a  standing  wave  pattern  in  front  of 
the  breakwater. 

CONCLUSIONS 

Theoretical  and  laboratory  studies  have  indicated  that 
the  Bragg  reflection  breakwater  concept  may  have  merit 
as  an  expedient  shore  protection  method.  The  wave  shel¬ 
tering  produced  by  a  small  number  of  bars  was  found  to 
significantly  reduce  the  estimated  erosion  volume  due  to 
a  model  storm.  The  small  scale  field  test,  however, 
served  to  demonstrate  many  of  the  difficulties  that  will 
need  to  be  overcome  before  the  Bragg  reflection  breakwa¬ 
ter  concept  becomes  practical.  Bar  stability  and  sea 
bed  scour  appear  to  be  the  most  troubling  problems. 
Further  tests  are  needed  to  determine  the  overall  merit 
of  the  concept. 
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Considerations  in  Using  Bragg  Reflection  for 
Storm  Erosion  Protection 


By  Junes  A.  BaOard,'  Member,  ASCE,  Jack  W.  DeVries,*  and 
James  T.  Kirby,*  Membu,  ASCE 

Aasnucr:  It  has  been  suggested  that  Bragg  reflection,  the  combined  coherent 
wave  reflection  from  a  few  low-lying  shore-^rallel  bars,  might  be  used  to  protect 
a  beach  against  storm-wave  attack.  Numerical  models  are  used  to  examine  two 
issues  relating  to  the  feasibility  of  this  concept:  the  degree  of  erosion  protection 
provided  by  reflecting  a  poitioo  of  the  incident  wave  energy,  and  the  degree  of 
wave  reflection  that  can  be  generated  by  a  bar  6eld  of  va^ng  geometry.  The 
results  show  that  a  Bragg-reflectioo  bar  6eld  must  reflect  about  one-quarter  of  the 
incident  wave  energy  to  provide  a  signiScant  measure  of  storm  erosion  protection. 
Bar  fields  with  uniform  spacing  are  capable  of  producing  the  required  magnitude 
of  wave  reflection,  but  lack  suifkieot  bandwidth.  Bandwidth  can  be  increased  by 
staggering  the  spacings  between  bars,  however,  this  produces  a  concurrent  decrease 
in  wave -reflection  magnitude.  Bragg-reflection  bar  fields  appear  capable  of  pro¬ 
viding  a  limited  measure  of  storm  erosion  protection  along  U.S.  Gulf  Coast  and 
Atlantic  Coast  beaches,  but  their  bandwidth  may  be  too  limited  for  use  along  Pacific 
Coast  beaches. 


Introduction 

Davies  ( 1980)  suggested  that  the  combined  action  of  a  few  low-lying  shore- 
parallel  bars  could  reflect  a  significant  portion  of  the  incident  wave  energy, 
thus  protecting  the  beach  behind  the  bars.  This  phenomenon,  dubbed  Bragg 
reflection  after  a  similar  phenomenon  in  optics,  has  received  considerable 
attention  in  the  recent  literature.  Theoretical  and  laboratory  studies  have 
shown  that  wave  energy  can  indeed  be  reflected  by  a  Bragg-reflection  bar 
field;  however,  the  following  question  remains:  is  it  practical  to  use  Bragg 
reflection  to  protect  a  beach  against  storm  wave  attack?  There  are  many 
issues  beyond  mere  reflection  that  have  a  bearing  on  this  question. 

For  Bragg  reflection  to  be  a  useful  shore-protection  method,  a  bar  field 
must  be  relatively  simple  to  construct,  it  must  produce  a  significant  degree 
of  erosion  protection,  and  it  must  be  adaptable  to  a  wide  range  of  wave 
and  beach  conditions.  Bailard  et  al.  (1990)  addressed  elements  of  the  first 
issue.  This  paper  addresses  elements  of  the  second  and  third  issues  via 
application  of  two  numerical  beach-profile  response  models  and  a  numerical 
wave-reflection  model. 

Background 

In  a  series  of  three  papers,  Davies  (1980,  1982a,  1982b)  considered  the 
two-dimensional  interaction  between  surface  waves  and  a  fixed  pattern  of 
undulations  on  the  seabed.  Each  undulation  reflects  a  small hmount  of  wave 
energy,  with  the  combined  reflection  from  a  series  of  undulations  being 
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either  coherent  or  destructive,  depending  on  the  phase  relationship.  The 
latter  is  a  function  of  the  ratio  of  the  surface  wavelength  to  the  spacing 
between  bed  undulations.  Davies  found  that  the  coefficient  of  wave  reflec¬ 
tion  was  oscillatory  with  respect  to  the  wavelength-to-spacing  ratio.  He  also 
found  that  a  resonant  condition  occurred  when  the  surface  wavelength  was 
exactly  twice  the  spacing  of  the  bed  undulations.  Thus,  it  appeared  that  for 
certain  wave  conditions,  a  few  shore-parallel  bars  might  reflect  a  substantial 
portion  of  the  incident  wave  energy. 

In  support  of  these  findings,  a  series  of  laboratory  experiments  (Heath- 
ershaw  1982;  Davies  and  Heathershaw  1984)  were  undertaken  to  measure 
the  reflection  of  monochromatic  waves  from  a  series  of  sinusoidal  bed  un¬ 
dulations.  The  results  confirmed  Davies's  theoretical  predictions  for  all  but 
the  resonant  condition. 

Mei  (1985)  developed  an  alternate  theory  that  applied  to  the  resonant 
reflection  condition.  Later.  Hara  and  Mei  (1987)  extended  this  work  to 
include  the  nonresonant  case  as  well.  In  the  meantime,  Kirby  (1986)  de¬ 
veloped  a  general  wave-equation  solution  that  was  applicable  to  both  res¬ 
onant  and  nonresonant  conditions. 

Directional  Waves  and  Sloping  Bottom 

Most  of  the  work  on  Bragg  reflection  has  been  done  in  two  dimensions 
using  normally  incident,  monochromatic  waves.  The  issue  of  nonnormal 
incidence  was  briefly  addressed,  in  theory,  by  Mei  (1985).  His  results  in¬ 
dicated  that,  relative  to  the  normally  incident  case,  the  effective  bar  spacing 
was  increased  by  an  amount  proportional  to  the  divergence  of  the  incident- 
wave  angle. 

The  effect  of  beach  slope  on  the  Bragg-reflection  process  was  considered 
by  both  Mei  (1985)  and  Kirby  (1986).  Eioth  found  that  a  sloping  bed  had 
the  effect  of  altering  the  ratio  of  the  local  wavelength  to  the  spacing  between 
bed  undulations.  By  varying  the  spacings  between  undulations  to  account 
for  the  local  change  in  wavelength,  results  were  obtained  that  were  identical 
to  the  constant  depth  case. 

Nonsinusoidal  Undulations 

A  sinusoid  is  not  a  practical  shape  for  an  artificial  Bragg-reflection  bar 
held.  Kirby  (1987)  and  Mei  et  al.  (1988)  suggested  that  reflection  from 
nonsinusoidal  undulations  could  be  calculated  by  expanding  the  shape  func¬ 
tion  of  the  bar  field  as  a  Fourier  series.  Kirby  and  Anton  (1990)  developed 
an  extension  to  the  nonresonant  interaction  theory  of  Davies  and  Heath¬ 
ershaw  (1984)  that  allowed  calculation  of  wave  reflection  from  the  individual 
Fourier  components  of  the  bottom  undulations.  A  series  of  laboratory  ex¬ 
periments  were  conducted  to  verify  their  theoretical  results.  They  found 
that  for  a  series  of  smewth  bumps  on  an  otherwise  flat  bottom,  the  peak 
reflection  was  reduced  slightly  relative  to  an  equivalent  sinusoidal  bottom. 

These  studies  provide  the  groundwork  for  consideration  of  Bragg  reflec¬ 
tion  as  a  shore-protection  alternative.  Bragg  reflection  appears  to  a 
workable  concept  in  theory  and  in  the  laboratory,  but  is  it  a  reasonable 
means  to  protect  a  natural'beach? 

PERFOnMANCe  CmTEMA 

Laboratory  studies  have  shown  that  a  Bragg-reflection  bar  field  will  mod¬ 
ify  the  incident-wave  field,  triggering  a  process  leading  to  the  formation  of 


a  tombolo  behind  the  bar  field  (Kirt^  and  Anton  1990).  The  process  begns 
with  the  bar  field  reflecting  a  portion  of  the  incident-wave  energy.  The 
resulting  wave  sheltering  pr^uces  a  long  shore  nadient  in  the  wave-driven 
radiation  stress,  causing  lateral  inflows  behind  the  bar  field  and  a  seaward¬ 
flowing  rip  current  at  its  center.  Laboratory  experiments  have  suggested 
qualitatively  that  this  circulation  cell  will  cause  sand  to  accumulate  behind 
the  bar  field,  leading  to  the  formation  of  a  tombolo. 

It  can  be  hypothesized  that  a  Bragg-reflection  bar  field  will  protect  a 
beach  against  storm-wave  attack  in  three  ways:  by  reducing  the  volume  of 
sand  er^ed  during  the  storm  as  a  result  of  wave  sheltering;  by  providing 
additional  sand  volume  available  for  erosion  via  formation  of  a  tombolo; 
and  by  restricting  the  loss  of  sand  offshore  via  formation  of  a  perched  beach. 
This  paper  deals  with  only  the  first  protective  mechanism;  however,  it  should 
be  noted  that  the  other  two  mecb^sms  may  have  a  significant  bearing  on 
the  ultimate  feasibility  of  the  Bragg-reflection  bar  concept. 

Beach  Erosion  Protection 

The  effect  of  a  Bragg-reflection  bar  field  in  reducing  the  volume  of  sand 
eroded  from  a  beach  by  storm  waves  was  investigated  using  two  numerical 
beach-profile  response  models.  The  first  model  was  developed  by  Kriebel 
(1982)  and  subsequently  included  in  the  Automated  G>astal  Engineering 
System  (1989),  published  by  the  U.S.  Army  Engineers  Coastal  Engineering 
Research  Center.  The  model  is  based  on  the  assumption  that  the  beach 
profile  evolves  towards  a  shape  having  a  uniform  distribution  of  wave-energy 
dissipation.  In  general,  the  model  has  been  shown  to  do  a  good  job  of 
predicting  erosion,  but  a  poorer  job  of  predicting  accretion  (Kriebel  1986; 
Birkemeier  et  al.  1987). 

The  second  beach-profile  response  model  was  developed  by  DeVries  and 
Bailard  (1988).  The  model  is  based  on  Bagnold's  (1%3)  energetics  sediment- 
transport  concepts,  as  generalized  by  Bailard  (1981).  Although  its  principal 
use  has  been  as  a  research  tool  for  understanding  cross-shore  sediment- 
transport  processes,  the  model  has  been  found  to  exhibit  realistic  erosion 
and  accretion  behavior  (Bailard  1985). 

These  models  were  used  to  predict  the  volume  of  sand  eroded  from  an 
initially  planar  beach  during  an  idealized  two-day  storm  (see  Fig.  1).  For 
simplicity,  it  was  assumed  that  the  sole  effect  of  the  Bragg-reflection  bar 
field  was  to  reduce  the  incident-wave  height.  The  transmitted-wave  height 
behind  the  bar  field  was  related  to  the  incident-wave  height  by  a  reflection 
coefficient,  K„  defined  as 

//?  =  (!-  KIW  . (1) 

where  H,  =  transmitted-wave  height;  and  H,  =  incident-wave  height.  Base¬ 
line  conditions  (i.e.,  no  bar  field  present)  were  computed  using  a  wave- 
reflection  coefficient  equal  to  zero. 

Storm  erosion  volumes  were  computed  for  a  range  of  sediment  sizes, 
initial  beach  slopes,  incident-wave  heights,  wave  periods,  and  wave- 
reflection  coefficients  (see  Tables  1  and  2).  The  response  of  the  beach  was 
computed  in  terms  of  a  normalized  erosion  volume,  defined  as 


where  V,  =  erosion  volume  with  wave  reflection  (i.e.,  with  a  bar  field 
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FIQ.  1.  Schdmatic  Otegram  of  Baach>Proflia  Raaponaa  to  Storm-Wava  Attack 


TABLE  1.  Paramatar  Valuaa  for  Baach*Proflla  Modaia 


Vartabla 

(1) 

(2) 

(3) 

(4) 

(5) 

Wave  height*  (ml 

2.00 

3.00 

imiH 

Wave  period*  (sec) 

6.00 

8.00 

— 

Grain  sue*  (mm) 

0.20 

0.40 

— 

Foreshore  slope* 

0.02 

0.03 

0.04 

0.05 

ReOection  coefficient* 

020 

050 

lEliSB 

095 

'Used  in  all  combinations. 


TABLE  2.  Paramatar  Valuaa  for  Baach*Proflla  Modaia 


Cortstanl 

0)  1 

(2) 

Storm  duratkw  (hr) 

48.0 

Berm  height  (m.  SWL) 

8.0 

Slope  of  berm  face 

0.08 

O&hore  slope 

0.0^ 

Initial  depth  of  slope  intercept  (m,  SWL) 

-2.0* 

Drag  ooefRcieni 

0.005* 

Be^oad  effideacy 

0.13* 

Suspended  load  efficieiicy 

0.01* 

'PcVfks  aad  Baflard  aiodd  ooly. 


I 


I 


I 
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present);  end  Vo  erosion  volume  without  wave  reflection  (i.e.,  baseline 
conditions). 

For  both  beach>proflle  response  models,  the  normalized  erosion  volume 
was  found  to  be  a  strong  function  of  the  reflection  coefficient  and  a  weaker 
function  of  the  sediment-^in  size,  the  near>shore  beach  slope  and  the 
incident-wave  height.  The  incident-wave  period  was  found  to  nave  a  neg- 
tigible  effect  on  the  normalized  erosion  volume. 

,  Fig.  2  shows  a  plot  of  the  predicted  normalized  erosion  volume  as  a 
function  of  the  wave-reflection  coefficient.  Due  to  the  weak  dependencv  of 
the  erosion  volume  on  the  sediment  grain  size,  beach  slope,  and  incident 
wave  height,  there  was  considerable  scatter  in  the  predicted  erosion  volumes 
for  each  model.  To  simplify  interpretation,  only  the  average  and  extreme 
values  from  each  model  were  plotted  in  Fig.  2.  For  both  m^els,  avera^ng 
was  done  over  all  cases  producing  erosion  of  the  beach. 

Reflection  Bandwidth 

A  Bragg-reflection  bar  field  will  typically  reflect  wave  energy  over  a 
relatively  nanow  range  of  wave  periods.  To  be  an  effective  shore-protection 
device,  a  bar  field  must  reflect  wave  energy  over  the  full  ranee  of  storm- 
wave  periods  anticipated  at  a  given  beach.  A  survey  was  made  of  typical 
wave  climates  along  three  U.S.  coastlines  to  determine  the  required  band- 
widths  for  potentialBragg-reflection  bar  fields.  Utilizing  data  from  Thomp¬ 
son  (1977)  and  Jensen  (1983),  wave  periods  were  noted  for  past  storms 
having  a  recurrence  interval  of  2  years  or  less.  Table  3  contains  a  list  of  the 
range  of  wave  periods  for  Atlantic  Coast,  Gulf  Coast,  and  Pacific  Coast 
storms.  Gulf  Coast  storms  have  the  narrowest  range  of  wave  periods  (S-9 
s).  with  Atlantic  Coast  storms  having  a  somewhat  wider  ranee  (6-11  s), 
and  Pacific  Coast  storms  having  the  widest  range  (6-17  s).  Normalizing 
these  wave-period  ranges  by  their  central  wave  period,  the  normalized  band- 
widths  of  storms  are  ^proximately  0.70  for  Gulf  Coast  and  Atlantic  Coast 
beaches,  and  1.0  for  Pacific  Coast  beaches. 

Reflection  Performance 

In  assessing  the  anticipated  performance  of  a  Bragg-reflection  bar  field, 
a  few  practical  limits  were  placed  upon  its  configuration.  First,  it  was  as- 
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TABLE  9.  Wavt  Ptrlod  Rangat  for  Commonly  Occurring  Storma  along  U.S. 
Coaatllnaa 


FIQ.  3.  Definition  Diagram  tor  Uniformly  Spaced  Bar  Geometry 


sumed  that  the  bar-height-to-water-depth  ratio.  BIh.  could  not  exceed  0.5. 
Beyond  this  vajue,  the  bars  begin  to  look  more  like  submerged  breakwaters 
than  undulations  on  the  seafloor.  Second,  it  was  assumed  that  a  bar  field 
could  include  no  more  than  9  bars.  Even  with  this  limitation,  seaward  bars 
tend  to  become  large  on  beaches  having  moderate  slopes  (for  constant 
I  BIh  ratio).  Finally,  it  was  assumed  that  the  bar  field  resided  in  shallow 

water.  This  last  assumption  ensures  that  the  wave-reflection  characteristics 
of  the  bar  field  are  solely  functions  of  the  bar  height-to-depth  ratio.  B//i. 
the  bar  width-to-spacing  ratio,  If/E.  and  the  total  number  of  bars. 

I  Uniform  Bar  Spacing 

The  reflection  characteristics  of  Bragg-reflection  bar  fields  having  a  uni¬ 
form  spacing  between  bar  elements  were  investigated  using  Kirby's  (1987) 
numerical  wave-reflection  model.  For  simplicity,  constant  depth  was  as¬ 
sumed.  however,  the  same  results  could  have  been  obtained  by  adjusting 
the  bar  spacing  to  account  for  wave  shoaling  on  a  sloping  bottom. 

I  Referring  to  Fig.  3,  the  bar  field  consisted  of  cosine-shaped  bumps  on  an 

otherwise  flat  bed  such  that  the  shape  of  the  bottom,  z{x),  can  be  described 
as: 

^  z{x)  =  Bsin  ^^(x  -  nZ.  +  (nL  >  j  <=  x  <=  nL  + 

. (3fl) 

z(x)  ■  0;  otherwise  . (3b) 

where  B  »  bar  height;  IF  -  bar  width;  L  -  bar  spacing;  and  n  -  number 
^  of  the  bar  element. 


Wive<icflection  properties  were  estimated  for  three  different  Bragg* 
reflectioo  bar  fields  over  a  ranae  of  bar  height-to-deplh  ratios  (BIh)  and 
bar  widtlhto-bar  spacing  ratios  (H7L).  Results  were  computed  for  bar  fields 
consistin|  of  three,  five,  and  nine  elements. 

Fig.  4ls  a  plot  of  the  normalized  reflection  bandwidth  as  a  function  of 
the  ratios  BIh  and  WIL  for  each  bar  field.  The  normalized  bandwidth  was 
defined  as  the  widui  of  the  principal  wave-reflection  peak  (in  s)  at  the  half¬ 
value  of  the  peak  reflection  coefficient,  divided  by  the  central  wave  period 
(see  Fig.  5). 

Fig.  6  is  a  plot  of  the  effective  wave-reflection  coefficient  as  a  function 
of  the  ratios  BIh  and  WH  for  each  bar  field.  The  effective  wave-reflection 
coefficient  was  defined  as  the  average  value  of  the  reflection  coefficient 
over  the  bandwidth  of  the  principal  reflection  peak  (see  Fig.  7). 

Staggered  Bar  Spacing 

Bragg-reflection  bar  fields  with  uniform  bar  spacings  have  relatively  nar¬ 
row  bandwidths  (0.1  <  BWITq  <  0.5).  These  bandwidths  are  less  than  the 
normalized  bandwidths  of  frequently  occurring  storms  on  most  U.S.  beaches 
(0.7- 1.0),  thus  limiting  the  utility  of  a  uniformly  spaced  bar  field.  It  is 
possible,  however,  to  alter  the  bandwidth  of  a  Bragg-reflection  bar  field  by 
staggering  the  spacings  between  bar  elements. 

Kirby  (1987)  briefly  considered  the  idea  of  staggered  bar  spacings  and 
found  that  the  bandwidth  of  a  Bragg-reflection  bar  field  could  be  increased, 
but  only  at  the  expense  of  reducing  the  wave -reflection  coefficient.  In  the 
present  study  this  concept  was  further  explored  by  systematically  varying 
the  spacing  geometry  of  three  Bragg-reflection  bar  fields. 

Effective  wave-reflection  coefficients  and  reflection  bandwidths  were  cal¬ 
culated  for  three-,  five-,  and  nine-element  bar  fields  having  constant  BIh 


no.  4.  Olm«nsionl«M  Reflactlon  Bandwidth  {BWITJ  vartut  Bar-Width-to-Bar- 
Spacing  Rrito  {WiL)  «id  Bar-Halght-to-Wator-Oapth  Ratio  (Bih)  for  Uniformly  Spaead 
Bar  Flalda  Compel  of  3,  S,  and  •  Etamanta 


FIG.  5.  Definition  Diagram  for  Effective  Wave-Reflection  Bandwidth 


3  Bars  s  bars  9  bars 


FIG.  C.  Effective  Coefficient  of  Wave  Reflection  veraus  Bar-Wkfth-to- 

Bar-Speclng  Ratio  {WiL)  and  Bar-Helght-to-Water-Oegth  Ratio  {B/h)  for  Uniformly 
Spaced  Bar  FleWa  Compoaed  of  3,  S.  and  t  Elements 

and  WIL  ratios  of  0.3  and  0.4.  respectively.  The  geometries  of  the  bar  fields 
were  varied,  as  shown  in  Fig.  8.  by  applying  a  stretching  parameter.  A, 
defined  as  follows: 

L'  •  Ul  ±  nA);  n  «  1.  2,  3.  4  . (4) 


99 


I 


> 


3  BARS 
- ?— 


(mt 


L - X - J 

m-w  m»*) 


5  BARS 


1 - L - - - - - X - - - - X - J 

UI-4AI  UI-34)  Ul-2«l  Ul-M  Ul*«  UI*2A)  U>*3>) 
t*l 


RQ.  1.  Definition  Dlegrem  for  Staggered  Spaced  Bar  Geometriea 


where  L'  -  the  local  staggered  bar  spacing;  and  L  =  the  equivalent  uniform 
bar  spacing  defined  as  the  overall  width  of  the  bar  field  divided  by  the 
numMr  of  bar  spaces. 

Fig.  9  shows  a  plot  of  the  predicted  bandwidths  for  each  bar  field  as  a 
function  of  the  stretching  parameter,  A.  Fig.  10  shows  a  plot  of  the  predicted 
)  wave-reflection  coefficient  for  each  bar  field  as  a  function  of  the  stretching 

parameter.  In  both  figures,  the  (vedicted  bandwidths  and  wave-reflection 
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FKa.  1  Rtlitive  Bandwidth  (BW/BlVa)  versus  Bar  Geometry  Stretching  Paramettr 
(A)  for  Staggered  Bar  FieMs 


FIG.  10.  Reiattve  Wave  Reflection  Coefficient  {KJKr^)  versus  Bar  Geometry 
Stretching  Parameter  (A)  for  Staggered  Bar  Fields 


coefficients  have  been  normalized  by  their  corresponding  values  for  uniform 
bar  spacings. 

^  Alternative  staggered  bar  geometries  were  also  evaluated,  however,  the 

*  best  results  were  quantitatively  similar  to  those  in  Figs.  9  and  10.  Similarly. 

the  ratio  BIh  was  found  to  have  a  relatively  small  influence  on  these  results. 

Comment 

^  The  aforementioned  results  suggest  that  a  Bragg-reflection  bar  field  re¬ 

duces  storm  erosion  a^roximately  in  proportion  to  the  percentage  of  wave 
energy  that  it  reflects.  Referring  to  Fig.  2.  a  bar  field  having  a  modest  wave- 
reflection  coefficient  of  0.2  will  reduce  the  predicted  storm  erosion  volume 
by  about  4%.  Similarly,  a  bar  field  having  a  wave-reflection  coefficient  of 
B  O.S  win  reduce  the  predicted  storm  erosion  volume  by  about  25%.  Con¬ 

sequently,  as  a  guideline,  it  appears  that  a  Bragg-reflection  bar  field  must 
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produce  «  wtve*reflectk>n  coefficiefit  of  about  O.S  in  order  to  be  an  effective 
shore-protection  device. 

Referring  to  Figs.  4  and  6.  Bragg-reflection  bar  fields  having  uniformly 
spaced  bar  elements  are  capable  of  producing  the  desired  level  of  wave 
reflection  {K,  «  O.S),  but  lack  sufficient  bandwidth  for  practical  application. 
Bar  fields  having  a  small  number  of  elements  have  a  wide  (althou^  insuf¬ 
ficient)  bandwidth,  but  a  modest  level  of  wave  reflection.  Bar  fields  having 
a  large  number  of  elements  have  a  high  level  of  wave  reflection,  but  very 
nanow  bandwidths. 

To  be  of  practical  use  along  U.S.  coastlines,  a  Bragg-reflection  bar  field 
must  have  a  dimensionless  bandwidth  of  at  least  0.7.  Fig.  4  indicates  that 
the  maximum  bandwidths  produced  by  the  uniformly  spaced  Bragg- 
reflection  bar  fields  range  from  a  high  of  0.45  (3  bars)  to  a  low  of  0.20 
(9  bars).  These  bandwid^s  can  be  significantly  increased  by  staggering  the 
spacing  between  bar  elements  (see  Fig.  9),  but  this  also  decreases  the  wave- 
reflection  magnitudes  (see  Fig.  10). 

Assuming  that  data  in  Figs.  9  and  10  can  be  applied  to  B/h  ratios  between 
0.1  and  O.S,  and  WiL  ratios  between  0.3  and  0.9  (spot  checks  indicate  that 
this  is  approximately  true),  maximum  bandwidths  and  wave-reflection  coef¬ 
ficients  were  estimated  for  staggered  bar  fields  consisting  of  three,  five,  and 
nine  elements.  Referring  to  Table  4,  it  was  found  that  minimally  acceptable 
performance  (i.e.,  a  wave-reflection  coefficient  of  0.5  and  a  normalized 
bandwidth  of  0.7)  was  just  possible  using  a  nine-element  staggered  bar  field 
having  a  BIh  ratio  of  O.S,  a  W/L  ratio  of  0.6,  and  a  A-value  of  0.07. 

These  findings  suggest  that  a  Bragg-reflection  bar  field  may  have  some 
potential  as  a  shore-protection  device  for  Gulf  Coast  and  Atlantic  Coast 
beaches,  but  little  potential  for  Pacific  Coast  beaches.  Considering  just  the 
effects  of  two-dimensional  wave  sheltering,  the  number  and  size  of  the  bar 
elements  required  to  produce  a  significant  level  of  storm  erosion  protection 
will  make  construction  of  a  Bragg-reflection  bar  field  a  significant  under¬ 
taking.  The  difficulty  and  expense  of  constructing  an  effective  Bragg- 
reflection  breakwater  could  easily  approach  that  of  a  conventional  detached 
breakwater,  however,  the  latter  has  greater  bandwidth  and  wave-reflection 
capabilities. 

The  foregoing  findings  are  based  on  a  consideration  of  simple  two- 
dimensional  wave  sheltering.  Other  mechanisms  that  may  provide  additional 
erosion  protection  include:  the  formation  of  a  tombolo  ^hind  the  bar  field, 
and  the  formation  of  a  perched  beach.  The  tombolo  would  provide  addi¬ 
tional  sand  to  the  beach  profile,  thus  reducing  the  magnitude  of  the  induced 
shoreline  retreat.  The  perched  l^ach  may  result  fi-om  the  retention  of  eroded 


TABLE  4.  Maximum  Availabla  Bandwidths  for  Staggared  Bar  Fields 


Site 

(1) 

Required 

bandwidth 

(s) 

(2) 

Availabis  Bandwidth* 

Three  bars 
(9) 

(3) 

Rve  bars 
(s) 

(4) 

Nine  bars 
(s) 

(5) 

Atlantic  Coast 

3.2 

3.4 

5.2 

Gulf  Coast 

HSH 

3.8 

4.2 

6.4 

Pacific  Coast 

5.2 

5.8 

8.6 

*Cslculsied  sssuming  K,  -  O.S. 


sand  against  the  shoreward  bar  element,  limiting  the  induced  shoreline 
retreat. 

Conclusions 

Bragg-reflection  bar  fields  appear  capable  of  providing  a  limited  measure 
of  storm-erosion  protection  on  Gulf  Coast  and  Atlantic  Coast  beaches. 
Sufficient  bandwidth  can  be  obtained  by  staggering  the  spacing  between  bar 
elements;  however,  adequate  levels  of  wave  reflection  will  require  more 
than  just  a  few  low-lying  bars. 

Numerical  model  results  indicate  that  a  staggered  nine-element  bar  field 
is  necessary  to  provide  a  25%  reduction  in  storm  erosion  volume  for  Gulf 
Coast  and  Atlantic  Coast  beaches.  To  produce  this  measure  of  protection, 
the  bar  elements  must  have  a  hei^t  of  approximately  half  the  local  water 
depth.  The  difficulty  in  constructing  a  bar  field  with  this  number  of  large 
bar  elements  suggests  that  the  Bragg-reflection  bar  concept  may  not  be 
practical  for  most  beaches. 

Laboratory  studies  indicate  that  a  Bragg-reflection  bar  field  may  produce 
a  tombolo  or  a  perched  beach  condition.  Although  outside  the  scope  of  this 
study,  these  effects  may  significantly  enhance  the  ultimate  feasibility  of  the 
Bragg-reflection  bar  concept. 
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Appendix  II.  Notations 

The  following  symbols  are  used  in  this  paper: 

B  =  bar  height; 

BW  =  bandwidth; 

BWq  =  bandwidth  with  uniformly  spaced  bars; 

Cf  =  drag  coefficient  for  bed; 
h  =  water  depth; 

Hi  =  incident-wave  height; 

H,  =  transmitted-wave  height; 

K,  -  wave-reflection  coefficient; 

=  reflection  coefficient  with  uniformly  spaced  bars; 

L  =  spacing  between  bar  crests; 
n  =  number  of  bar  element; 

To  =  wave  period  at  central  reflection  peak; 

=  normalized  beach-erosion  volume,  V^Vq; 

Vo  =  beach-erosion  volume  without  wave  reflection; 

V,  =  beach-erosion  volume  with  wave  reflection; 

W  =  bar  width; 

X  =  horizontal  distance; 
z  =  vertical  elevation; 

A  =  staggered  bar  stretching  parameter; 
e*  =  bed-load  efficiency; 
e,  =  suspended  load  efficiency;  and 
w  »  3.14. 
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Appendix  D:  A  Note  on  Bragg  Scattering  of  Surface  Waves 
by  Sinusoidal  Bars 

I  to  be  submitted  to  Physics  of  Fluids. 
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A  note  on  Bragg  scattering  of  surface  waves  by  sinusoidal  bars 
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Abstract 

The  reflection  of  linear  surface  waves  by  sinusoidal  bottom  undulations  is  considered 
in  the  case  where  the  incident  wave  is  not  necessarily  close  to  the  resonant  frequency. 
For  finite  detuning  away  from  the  resonant  frequency,  two  previous  solutions  are  shown 
to  give  results  which  are  inconsistant  with  direct  numerical  solutions,  especially  when 
the  results  are  extended  to  oblique  incidence.  The  correction  to  the  methods  is  given, 
and  various  consequences  of  the  new  results  are  examined. 


1  Introduction 

In  recent  years,  the  problem  of  the  reflection  of  surface  waves  by  undular  bottom  features 
has  drawn  considerable  attention,  owing  primarily  to  the  mechanism's  possible  relation  to 
coastal  morphology.  Kirby  (1986)  has  derived  an  extension  to  the  usual  mild-slope  equation 
(Berkhoff,  1972)  which  provides  a  theoretical  umbrella  for  most  of  the  existing  analytical 
results.  Davies  and  Heathershaw  (1984)  solved  for  the  reflection  coeflicient  for  the  case 
of  weak  reflection  (far  from  resonance)  and  sinusoidal  bars  in  otherwise  constant  depth. 
Their  solution  was  recovered  from  Kirby’s  equation  by  Kirby  and  Anton  (1990),  who  also 
provided  the  extension  to  the  case  of  oblique  incidence  over  the  bar  field.  Close  to  the 
resonant  condition,  where  the  bar  wavenumber  is  twice  the  wavenumber  component  of  the 
surface  wave  in  the  direction  normal  to  bar  crests,  Mei  (1985)  provided  a  solution  using 
a  multiple  scales  approach.  His  solution  was  extended  to  oblique  incidence  by  Dalrymple 
and  Kirby  (1986)  and  Mei  et  al  (1988,  hereafter  referred  to  as  MHN).  Kirby  (1986)  showed 
that  Mei’s  analysis  could  be  obtained  from  the  extended  mild-slope  formulation.  Further 
explication  of  the  problem  (along  with  some  additional  solution  techniques)  may  be  found 
in  Davies  et  al  (1989),  Kirby  (1989),  and  Benjamin  et  al  (1987).  As  in  these  papers,  we 
restrict  our  attention  here  to  the  case  of  perfectly  sinusoidal  bars. 

The  resonant  reflection  solution  of  Mei  proceeds  by  assuming  that  the  waves  are  de¬ 
scribed  by  carrier  waves  which  are  exactly  in  resonance  with  the  bar  field,  together  with 
amplitudes  which  vary  slowly  in  both  space  and  time,  to  account  for  both  frequency- 
wavenumber  detuning  and  slow  evdution  over  the  bar  field.  The  analysis  results  in  a 
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set  of  perfectly  tuned,  coupled  evolution  equations  for  the  wave  amplitudes;  these  may  be 
solved  for  the  reflection  coefficient  and  for  the  wave  heights  over  the  bar  in  the  case  of 
time-periodic  motions. 

For  the  case  of  finite,  arbitrary  detuning,  it  is  natural  to  extend  the  approach  given 
by  Mei  to  consider  the  frequency  and  wavenumber  detuning  as  part  of  the  carrier  wave; 
as  a  result,  the  coupled  amplitude  equations  become  explicitely  detuned,  but  may  still  be 
easily  solved  for  the  case  of  periodic  waves.  Liu  (1987)  and  Yoon  and  Liu  (1987,  referred 
to  as  YL)  took  this  approach  in  two  related  studies.  Liu  considered  the  resonant  reflection 
of  linear  monochromatic  waves  in  a  channel  with  corrugated  sides  and  bottom,  while  YL 
considered  the  scattering  of  incident  cnoidal  waves  by  a  sinusoidal  bar  field  in  the  context 
of  Boussinesq  theory.  (The  latter  problem  has  also  been  addressed  by  Kirby  and  Vengayil, 
1988).  For  the  case  of  normal  incidence  on  the  bar  field,  Lin’s  solution  gave  “cut-ofi^ 
conditions  (demarcating  the  boundary  between  exponential  and  sinusoidal  behavior  of  the 
incident  and  reflected  wave  envelopes)  which  differed  markedly  from  the  values  given  by 
Mei  (1985).  These  differences  were  attributed  to  the  effect  of  finite  detuning. 

In  this  note,  we  formulate  the  problem  of  detuned  scattering  by  a  field  of  sinusoidal  bars, 
and  obtain  a  solution  with  different  properties  than  would  be  found  using  the  methods  in 
Liu  and  YL.  In  section  2,  a  new  formulation  of  the  detuned-interaction  theory  is  developed. 
This  is  followed  in  section  3  by  a  discussion  of  previous  results  and  of  their  breakdown  in  the 
oblique  incidence  case.  The  methods  of  Liu  and  YL  are  extended  to  the  oblique  incidence 
case  in  the  Appendix.  In  section  4,  numerical  results  are  used  to  establish  the  validity  of 
the  present  formulation.  We  also  show  that  the  resonant  interaction  theory  of  Mei  (1985) 
is  reasonably  robust  for  the  entire  range  of  physically  relevant  cases. 

2  Solution  of  the  generalized  mild-slope  equation  for  de¬ 
tuned  resonance 

In  this  section,  we  present  a  perturbation  solution  of  the  equation  of  Kirby  (1986)  which 
differs  from  the  previous  detuned  interaction  results  given  by  Liu  and  YL.  The  deriv^ation 
here  follows  the  notation  of  MHN  closely,  although  the  governing  modulation  equations  are 
arrived  at  by  a  different  method.  We  consider  the  case  of  a  patch  of  sinusoidal  bars  ^(x) 
given  by 

«(x)  =  ZJsin(Ax);  A  =  ^  (1) 

id 

where  L  is  the  bar  wavelength  and  D  is  the  bar  amplitude.  A  finite  patch  of  n  bars  rests  in 
the  interval  0  <  x  <  nL,  superimposed  on  a  uniform  depth  h.  We  generalize  the  previous 
results  somewhat  by  considering  the  case  of  oblique  incidence  on  the  bar  field,  and  take  9 
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to  represent  the  angle  between  z  and  the  incident  wavenumber  vector.  The  general  two- 
dimensional  form  of  the  governing  equation  is  then  written  as 


where 


V(pV^)  +  k^CCg^  =  0 

p  =  cc,- 


(2) 


(3) 


cosh^  kh 

and  where  the  angular  frequency  w,  wavenumber  ifc,  phase  speed  C,  and  group  velocity  Cg 
are  determined  based  on  the  mean  depth  h  and  the  actual  (tuned  or  detuned)  wave  period. 
For  obliquely-incident  wives,  we  take 


^  =  ^z‘)e*”**';  Tn  =  k  sind. 

The  governing  equation  is  reduced  to  the  second-order  ODE  in  z; 


where  the  factor  7^  is  given  by 


4>xx  +  —<l>x  +  7^^  =  0 

P 


7^  =  k^p  ^CCg  - 


which,  for  the  case  of  small  bottom  perturbations,  may  be  approximated  by 


gkH 


CCg  cosb^  kh' 


I  =  \/k'^  -  =  k  cos  0 


As  in  Kirby  (1986),  we  seek  coupled  first  order  equations  of  the  form 

Oj  =  i70''’ -I- F'{c!>'*',o~) 

<Px  =  -*7^"  - 

where  the  total  wave  field  is  divided  into  incident  and  reflected  components, 

<t>  =  <i>+  +  (p- 


(4) 

(5) 

(6) 

(7) 

(8) 
(9) 

(10) 


and  -l-(-)  denotes  the  incident  (reflected)  wave.  Substituting  (8)-(10)  in  (5)  leads  to  the 
result 

(11) 

To  the  leading  order  of  approximation  in  powers  of  S,  the  coefficients  in  (8)  and  (9)  are 
given  by 


^  ^  2CCg  cosh*  kh  cos*  d) 


(12) 


(7P)x 

27P 


2qQq6x 

~kC^ 


(13) 


where 


(14) 


_  cos 29 
^  ~  €08^9 
and 

Note  that  Hq  here  differs  from  the  value  given  in  Mei  (1985)  in  that  it  is  evaluated  at 
the  detuned  wavenumber  rather  than  the  resonant  wavenumber.  We  also  note  that  no 
approximations  have  been  made  in  obtaining  the  basic  form  of  the  coupled  equations  (8) 
and  (9);  the  splitting  is  reversible  and  no  information  has  been  suppressed.  We  introduce 
the  explicit  form  of  6x: 

^x=^(e'**+e-^)  (16) 

and  seek  a  solution  for  slowly- varying  incident  and  reflected  waves  of  the  form 


—ik  cot  Sx 


=  >l(i)e’^  =  >i(x)e‘*“'**;  4>-  =  B(i)e-'*  =  5(i)e 


Following  Liu  and  YL,  we  construct  am  approximate  detuned-interaction  model  by  sub¬ 
stituting  (17)  in  (8)-(9)  and  retaining  terms  which  come  closest  to  satisfying  resonance 
conditions.  This  leads  to  the  coupled  evolution  equations 


Xx  =  - 


Bx  =  - 


flflOCOS  9  n 


iloocosB  (k, 


where 


is  the  resonant  wavenumber,  and 


kres  — 


2  cos  9 


=  X^2k  cos  9 


is  the  detuning  parameter,  as  introduced  in  Liu  and  YL.  For  I  =  fccos9  =  A/2,  we  re¬ 
cover  the  exactly-resonant  interaction  equations  following  from  MHN,  equation  (2.5).  The 
detuned  equations  differ  from  the  results  in  MHN  by  the  inclusion  of  the  factor  {kre$/k), 
which  introduces  am  aisymmetry  in  the  coupling  coefficient.  We  see  that  waves  with  lower 
wavenumbers  (indicating  relatively  shaUower  water)  become  more  strongly  coupled,  as  may 
be  expected  on  physical  grounds. 

We  can  increase  the  correspondence  between  the  results  of  MHN  and  the  present  deriva¬ 
tion  by  introducing  a  detuning  wavenumber 

K  =  k-kr^  (22) 
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&Dd  a  corresponding  frequency  parameter 

ft  =  KCg  (23) 

Note  that  ft  again  differs  from  Mei’s  ft  since  detuning  is  finite  and  the  relation  between  the 
detuned  wavenumber  and  the  detuned  frequency  here  is  not  linear,  except  in  the  limit  of 
extremely  shallow  water.  The  detuning  parameter  &  may  now  be  written  as 


=  —2K  COS0. 


We  introduce  the  transformation 


(;)■(* 


ie-‘^co«s 
^^ifC  cot  9 


and  obtain  the  coupled  equations 


A  ■L'  a  A  ftoOCOS^  /*re,\  5 

A.  =  ,Kco.0A 

(26) 

iA-cos.B 

(27) 

We  may  reduce  these  to  two  second-order  equations  for  A  and  B: 

ii) 

(28) 

1 

\°  /  XX  \°  / 

The  parameter  P  is  determined  from 

P  ftocos^f/fty 

c,  [UJ  u  jj 

(29) 

and  may  be  either  real  or  imaginary  depending  on  the  term  in  square  brackets.  The  ex¬ 
pression  for  P  here  is  also  similar  to  the  expression  given  in  MHN,  with  the  exception  of 
the  appearance  of  the  krt,/k  ratio.  The  cutoff  condition  in  the  solution  corresponds  to  the 
value  P  =  0. 

Equations  (28)  are  solved  with  the  boundary  conditions 

/4(0)  =  Aq-,  B{nL)  =  0  (30) 

which  indicates  an  incident  wave  amplitude  of  Aq  at  the  start  of  the  bar  field  and  a  reflection 
of  zero  at  the  downwave  end.  The  solution  is  given  by 

.  _.iKc^t,PCgCwP{nL-x)-iUcosBs\nP{nL-x) 

A[x}  -  >lo€  _  in  cos  sin  nPL 

and 

-  A  -iKco,9.^{lfres/k)cos$siaP{nL-x) 

B(x)  -  Aoe  pCgCOinPL  -  ift costf  sin nPL' 
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The  reflection  coefficient  upwave  of  the  bar  field  is  given  by 
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► 


\ 


1 
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> 


5(0) 

(loQikretlk)  cos^sin  nPL 

Ao 

PCf  cos  nPL  -  I'fi  cos  ^  sin  nPL 

We  illustrate  the  present  solution  in  comparison  to  the  resonant  solution  of  MHN  by 
plotting  contours  of  reflection  coefficient  as  a  function  of  0  and  2fccos0/A,  where  resonance 
occurs  at  2A:cos^/A  =  1.  The  solution  of  MHN  is  implemented  exactly  as  given  in  their 
paper,  with  K  specified  based  on  the  chosen  wavenumber  parameter  and  fl  computed  from 
Q  =  KCg^^rea)-  The  choice  of  21:cosd/A  as  a  plotting  parsimeter  is  inconvenient  from  a 
design  point  of  view,  where  a  fixed  range  of  wave  periods  would  lead  to  consideration  of  a 
fixed  range  of  2k  jX.  This  alternate  choice  will  be  followed  below;  it  is  avoided  in  the  initial 
comparisons  here  since  the  resonance  curve  goes  to  large  values  of  21;/A  as  0  increases. 

To  maintain  correspondence  with  experimental  results  given  in  Davies  and  Heathershaw 
(1984),  we  choose  n  =  A,D/h  =  0.32,  Z)  =  5cm,  L  =  Im,  and  h  =  15.625cm.  Solutions  are 
computed  for  the  parameter  ramge  .25  <  2A:cosd/A  <  1.75,  0“  <  ^  <  85®.  Solutions  for 
the  resonant  theory  of  MHN  are  shown  in  Figure  1.  The  symmetry  of  the  theory  about 
the  resonance  condition  2kcos0fX  =  1  is  apparent,  as  is  a  tendency  for  the  solution  to  die 
out  for  large  values  of  0.  Figure  2  shows  reflection  coefficient  contours  and  the  locus  of  the 
cutoff"  condition  for  the  present  analytic  theory.  The  results  here  are  markedly  asymmetric 
about  the  resonance  condition.  Agreement  between  the  two  analytic  solutions  is  good,  both 
for  reflection  coefficients  and  cutoff"  conditions,  up  to  9  =  55®.  For  larger  O's  and  for  values 
of  2kcos0/\  <  0.5  and  >  1.5.  the  asymmetry  of  the  solution  becomes  marked  and  the  two 
approximate  solutions  deviate. 

The  apparent  agreement  between  the  cutoff" conditions  predicted  here  and  the  conditions 
predicted  by  the  Mei  (1985)  and  MHN  theory  is  in  marked  contrast  to  the  large  deviation 
between  cutoffs  seen  in  comparing  the  theory  of  Liu  and  Mei  (1985)  (see  section  4  in  Liu). 
Direct  numerical  determination  of  the  detuning  0  corresponding  to  the  cutoff"  condition  has 
been  made  for  the  bar  geometry  considered  here,  but  with  a  range  of  D/h  values.  The 
cutoff"  condition  on  the  high  frequency  side  of  resonance  is  compared  to  Mei's  (1985)  theory 
in  Figure  3,  where  we  plot  the  %  error  in  Mei’s  result  relative  to  the  present  theory.  The 
results  typically  agree  to  within  1%  or  so.  in  contrast,  Liu  shows  deviations  on  the  order 
of  50%  between  his  and  Mei’s  predicted  cutoffs. 

For  angles  less  than  about  60®,  we  will  evaluate  the  disagreement  between  the  approxi¬ 
mate  solutions  here  by  comparison  with  numerical  solutions  for  the  fuU  problem.  We  close 
this  section  by  remarking  that  both  of  the  approximate  solutions  studied  here  would  be 
expected  to  break  down  far  from  resonance,  since  the  effect  of  neglected  non-resonant  terms 
in  the  equations  would  become  as  important  as  the  effects  of  the  retained  terms. 
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3  Comparison  with  previous  results 


In  this  section,  we  compare  the  present  results  to  the  results  obtained  by  extending  the 
derivation  techniques  used  in  Liu  and  YL.  Those  derivations  are  described  in  Appendix  A. 
The  resulting  derivations  produce  coupled  equations  with  the  same  form  as  equations  (18)- 
(19)  above,  and  hence  we  can  consider  the  correspondence  between  the  form  of  the  coupling 
coefficients,  the  values  of  the  factor  P  and  loci  of  the  cutoff  conditions,  and  numerical 
values  of  predicted  reflection  coefficients.  We  limit  the  comparison  here  to  an  examination 
of  the  coupling  coefficients  and,  in  particular,  to  the  determination  of  the  zeroes  of  those 
coefficients,  where  we  see  a  decoupling  between  the  wave  and  bar  fields.  Both  the  Liu 
and  YL  theories  have  the  interesting  feature  of  predicting  decoupling  for  normally  incident 
waves  at  a  finite  value  of  the  parameter  2k/ X.  When  extended  to  oblique  incidence,  both 
theories  produce  patterns  of  decoupling  which  are  at  odds  with  the  present  theory,  with  the 
MHN  resonant  theory,  with  the  nonresonant  theory  of  Kirby  and  Anton  (1990),  and  with 
numerical  results. 

The  interaction  coefficient  for  the  present  theory  may  be  written  as 


QoCOs2^ 
CgCosO  \  k  ) 


(34) 


which  has  a  zero  at  ^  =  45*  for  all  incident  wave  frequencies.  This  result  also  was  found 
by  MHN  for  the  resonant  case,  and  is  a  feature  of  a  wide  range  of  non-resonant  scattering 
theories  for  small  obstacles  (see  Kirby  and  Anton,  1990,  for  a  review).  The  locus  of  the  zero 
in  the  present  coupling  coefficient  is  shown  in  Figure  4  along  with  the  locus  of  the  resonance 
curve  2k /X  =  l/cos0. 

The  coupling  coefficient  Oj  obtained  from  the  YL  derivation  has  zeroes  on  the  curve 
2k /X  =  2  cos  9,  shown  in  Figure  4.  The  intersection  of  this  curve  with  the  resonance  curve 
occurs  at  ^  =  45“  and  2k/X  =  \/2,  as  does  the  curve  for  the  zeroes  of  the  present  theory. 
Away  from  this  intersection,  the  predicted  zeroes  deviate  markedly. 

For  the  second  theory  of  Appendix  A,  the  interaction  coefficient  02  leads  to  decoupling 
on  a  curve  2k/ X  =  v^2)/cos^,  which  lies  to  the  right  of  the  resonant  interaction  curve  in 
Figure  4  and  doesn’t  intersect  it  anywhere. 

The  two  theories  described  in  Appendix  A  clearly  predict  patterns  of  wave-bottom 
interaction  which  are  at  odds  with  the  remaining  body  of  theory  on  this  problem.  This 
result  points  out  the  sensitivity  of  the  results  here  to  the  method  of  derivation;  in  particular, 
theories  which  neglect  terms  of  the  order  of  terms  that  appear  eventually  in  the  equations 
to  be  solved  (Az,  and  Bxx  here)  lead  to  incorrect  results.  The  splitting  method  employed 
by  Kirby  (1986)  and  Kirby  and  Vengayil  (1988)  appears  to  be  a  robust  method. 
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4  Numerical  results 


Numerical  solutions  are  obtained  for  the  parameter  range  0.5  <  2A;/A  <  2.5,  0°  <  0  <  85°, 
using  second-order  accurate  centered  finite  differences  applied  to  (5).  Boundary  conditions 
for  the  numerical  problem  are  discussed  in  Kirby  (1989).  Results  for  the  n  =  4  Davies  and 
Heathershaw  barfield  are  shown  in  Figure  5.  We  see  that  the  numerically  predicted  solution 
shows  nearly  total  reflection  at  large  angles  of  incidence,  as  well  as  very  small  reflection 
at  ^  =  45°,  as  predicted  by  the  analytic  theories  developed  here.  The  reflection  coefficient 
does  not  drop  identically  to  zero  for  6  =  45“  owing  to  the  effect  of  nonresonant  modes  in 
the  solution  which  are  not  accounted  for  in  the  analytic  approaches. 

In  Figure  6,  reflection  coefiicents  for  normally  incident  waves  obtained  numerically, 
from  Mei's  solution,  and  from  the  present  analytical  solution  are  compared.  The  numerical 
solution  (solid  line)  is  seen  to  obtain  a  somewhat  higher  maximum  than  in  the  analytic 
solutions,  and  the  maximum  is  shifted  to  a  slightly  lower  wavenumber.  The  two  analjrtic 
solutions  agree  at  resonance  but  deviate  sbghtly  away  from  resonance.  The  numerical 
solution  also  indicates  that  reflection  is  slightly  stronger  for  the  first  lobe  on  the  short  wave 
side  of  the  main  resonance,  in  contrast  to  the  prediction  of  the  present  anal\'tic  theory. 
The  solution  of  Mei  is  therefore  a  somewhat  more  accurate  representation  of  the  numerical 
solution,  although  the  two  analytic  results  agree  more  closely  with  each  other  than  they 
do  with  the  numerical  result.  Overall,  the  deviation  between  the  numerical  and  analytical 
results  would  not  be  resolvable  within  the  accuracy  of  the  Davies  and  Heathershaw  (1984) 
data  set. 

5  Effects  of  frictional  damping 

The  propagation  of  waves  over  a  bar  field  is  affected  by  a  number  of  mechanisms  which 
are  not  included  in  the  analysis  presented  above.  YL  and  Vengayil  and  Kirby  (1988)  have 
discussed  some  of  the  consequences  of  nonlinear  interaction  on  the  scattering  process.  In 
real  applications,  it  is  expected  that  a  number  of  dissipative  mechanisms  could  come  into 
play.  Bottom  boundary  layer  damping  would  affect  the  waves  moderately,  causing  a  small 
decrease  in  the  transmitted  and  reflected  waves  predicted  by  the  theory  described  above. 
For  bars  fronting  a  beach  face  in  relatively  shallow  water,  the  protruding  bar  crests  could 
also  induce  breaking  of  steep  incident  waves,  leading  to  significant  reduction  of  incident 
energy  over  the  bar  Add. 

Here,  we  consider  the  effect  of  laminar  frictional  damping  on  the  predicted  reflection 
coefficient.  Following  Booij  (1981)  and  Kirby  (1986),  we  may  extend  the  generalized  mild- 


8 


slope  equation  to  include  damping  effects  according  to 


-  V  •  (pV^)  +  (u>*  -  k^CCg)4  =  0 

or,  for  time-harmonic  waves, 

V  •  (pV<A)  -H  k^CCg^  +  iuW^  =  0 


(35) 


(36) 


Here,  VF  is  a  generalized  complex  damping  term.  Following  the  derivation  in  section  2  leads 
to  the  coupled  evolution  equations 


A,  =  - 


W 


2Cg  COS  0 


A  +  iK  cob6A  — 


floor  cos  d 


(¥)^ 


W 


2Cg  cos  0 


B  -  iK  COS0B  - 


^9 

floor  cos d  fkreM\ 
Cg  VTJ 


(37) 

(38) 


where  we  have  assumed  that  the  dissipative  effects  act  uniformly  and  equally  on  the  incident 
and  reflected  waves.  The  new  dissipative  term  may  be  absorbed  in  the  previous  derivation 
by  defining  a  new  complex  detuning  factor 

W 


Cg  COS  6 


=  —2K  COS  d  -  t 


Cg  COS  0 


(39) 


From  this  point,  the  derivations  in  section  2  follow  through  exactly,  with  the  exception  that 
the  definition  of  P  in  (29)  is  replaced  by 


P  = 


floCOS® 


V 

.flo/ 


+ 1 


iiw 


fig  cos*  6 


\2flocos*^/  \  k  ) 


(40) 


For  the  case  of  laminar  bottom  boundary  layer  damping,  an  expression  for  W  is  given 
bv  (Liu,  1986) 

gk^ 


VF  =  (l-i)- 


_/Z 

i^khy2u 


(41) 


wcosh 

For  laminar  damping  alone,  the  effects  on  the  results  given  above  are  quite  small,  owing 
to  the  relatively  short  length  of  the  bar  field.  The  reflection  coefficients  with  damping  can 
be  either  smaller  or  larger  than  the  coefficients  without  damping,  due  to  the  effect  that 
damping  has  on  the  wave  speed  and  wavelength.  It  is  not  anticipated  that  local  damping 
over  the  bar  field  has  any  significant  effects  in  any  of  the  available  laboratory  data  sets. 


6  Discussion 

The  results  of  this  and  previous  studies  have  shown  that  numerical  solutions  of  a  generalized 
mild-slope  equation  may  be  used  to  model  the  reflection  of  waves  from  simple  (nearly 
sinusoidal)  bar  fields.  A  body  of  analytic  approximations  are  also  now  available  for  providing 
estimates  in  the  same  situation. 
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Several  recent  studies  have  indicated  that  bottom  configurations  having  multiple  Fourier 
components  can  lead  to  additional  reflection  peaks  appearing  at  difference  wavenumbers 
resulting  from  the  interaction  of  two  or  more  individual  bottom  components  (Belzons  et 
al,  1991).  This  situation  is  not  handled  properly  by  the  existing  generalized  mild-slope 
equation.  The  problem  of  these  subharmonic  resonances  is  being  pursued  further  by  the 
present  author  in  the  context  of  a  generalization  of  the  mild-slope  equation  to  handle  the 
presence  of  non-propagated  modes. 

Acknowledgement.  This  work  was  supported  by  the  Office  of  Naval  Research,  con¬ 
tract  N00014-90-J-1678. 

Appendix  A.  Extension  of  Liu  and  YL  results  to  the  case  of 
obliquely  incident  waves 

The  problem  of  specifying  the  evolution  equations  for  the  detuned  interaction  between 
surface  waves  and  the  sinusoidal  bar  field  has  been  considered  previously  by  Liu  and  YL,  as 
well  as  by  Kirby  and  Vengayil  (1988).  The  study  of  Liu  was  restricted  to  normal  incidence 
of  linear  waves  in  a  channel  with  both  bottom  and  side-wall  corregations,  while  the  theory 
of  YL  was  ior  the  same  channel  configuration  but  for  weakly  nonlinear,  weakly  dispersive 
waves  obeying  the  Boussinesq  equations.  Here,  we  consider  the  case  of  normal  or  oblique 
incidence  of  linear  dispersive  waves  over  an  infinitely  wide  bar  field,  as  in  section  2,  and 
extend  the  theories  of  Liu  and  YL  accordingly.  Section  3  of  the  main  text  gives  a  comparison 
between  various  features  of  these  extended  results  and  the  results  in  section  2. 

The  YL  approach  is  conceptually  simpler  than  the  Liu  approach,  and  so  we  follow  it 
first.  We  start  by  rewriting  the  governing  equation  as 

where  we  have  used  the  fact  that  the  mean  depth  is  taken  to  be  uniform  and  the  bar  fidd 
varies  as  S(x)  only.  Employing  the  separation  (4)  in  the  main  text  reduces  (Al)  to 

"  u,cJosbHh 

We  now  make  the  separation  into  incident  and  reflecte-i  waves  with  slowly  varying  ampli¬ 
tudes,  as  specified  in  (10)  and  (17).  These  are  substituted  directly  into  (A2)  without  a  prior 
splitting.  An  ^parent  ordering  is  made  which  equates  the  importance  of  slow  derivatives 
of  A  and  B  with  terms  that  are  first  order  in  h.  Terms  of  0(^^)  are  then  dropped.  Then, 
substituting  the  expression  (1)  for  £  and  retaining  terms  closest  to  resonance  leads  to  the 
coupled  equations 

A,  = 


(A3) 


where 


oi  =  (costf)"^  ^2  ^  -  1^ 


and  where  all  other  notations  are  as  in  the  main  text.  The  resulting  model  is  similar  in  form 
to  (18)-(19)  in  the  main  text,  except  for  a  change  in  the  interaction  coefficient.  We  remark 
that  the  principle  reason  for  the  differences  between  the  present  model  and  the  model  in 
the  main  text  arrises  due  to  neglect  of  terms  involving  the  factors  6 Ax  and  6Bx.  If  these 
terms  were  retained  by  effectively  dividing  through  by  the  quantity  p  as  in  the  main  text, 
the  resulting  model  would  be  equivalent. 

For  the  case  of  normal  incidence  and  in  the  limit  of  long  wave  theory,  in  interaction 
coefficient  becomes 

(A6) 


as  in  YL.  At  resonance,  this  reduces  to  XD/Sh  as  in  Kirby  and  Vengayil  (1988),  but  it 
disagrees  with  their  theory  (again,  based  on  a  splitting  procedure)  away  from  resonance. 
Turning  to  the  work  of  Liu,  we  rewrite  (A2)  as 

ip4>x)s  +  l^CCa<i>  =  0  (A7) 

We  then  employ  the  Liouville  transformation 

=  (A8) 


and  obtain 


C*x-2P“‘PxxC  +  /V*C’QC  =  0 


Retaining  terms  to  0(j),  this  becomes 


(AlO) 


We  then  employ  the  split  into  incident  and  reflected  waves  according  to 


C  =  +  V’  {i)e  *'*. 


(All) 


Direct  substitution  and  the  neglect  of  second-order  terms  in  the  slow  derivatives  leads  to 
the  set  of  coupled  equations 


0+  =  _-2a20"e’^* 


(A12) 

(A13) 
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where 


I 


I 


02  =  COS  6 


(A14) 


For  the  case  of  normal  incidence,  the  interaction  coefficient  in  (A12-A13)  reduces  to 


flo  ^^0 


(A15) 


as  in  Liu. 

The  main  drawback  in  the  derivation  of  (A12-A13)  lies  in  the  fact  that  the  amplitudes 
and  are  not  properly  decoupled.  For  example,  compare  the  transformed  incident 
wave  to  an  untransformed  form: 


I  =  p^/^Ae’'*  (A16) 

TaJdng  x  derivatives  of  both  sides  of  (A16)  leads  to  the  expression 

+  ^p~^^^pxA.  (A17) 

We  see  that  the  derivative  of  V’"*'  contains  a  term  which  would  resonate  with  the  reflected, 
B  wave,  indicating  that  the  incident-reflected  wave  separation  is  incomplete  in  this  method. 

I 

» 

I 
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Abstract 

This  paper  describes  the  development  of  a  model  fcH-  the  shoaling  and  refraction  of 
an  incident  directional  spectrum  over  beach  topography.  The  present  model  is  limited 
to  topography  which  varies  only  in  the  on-offrhore  direction,  but  no  restriction  is  made 
on  allowed  angles  of  incidence  with  respect  to  the  shore-normal  direction.  Several 
approximations  for  restricted  directional  bandwidth  are  discussed,  and  several  numerical 
examples  are  presented  which  illustrate  the  growth  of  signihcant  harmonics  resultmg 
from  the  interaction  of  waves  travelling  at  different  directions.  The  model  is  verified  in 
comparison  to  laboratory  data  for  Mach  reflection  of  cnoidal  waves  from  a  vertical  plane 
wall.  It  is  shown  that  the  model  provides  a  more  accurate  representation  of  the  evolving 
wave  field  than  does  an  earlier  parabolic  approximation,  using  the  same  laboratory  data. 


1  Introduction 

As  ocean  surface  waves  propagate  towards  shore,  they  pass  through  a  shoaling  zone  prior 
to  breaking  in  which  nonlinear  interactions  become  strong  and  can  significantly  modify  the 
wave  train.  This  zone  is  characterized  by  weak  frequency  dispersion,  since  waves  become 
relatively  long  compared  to  local  water  depth.  In  addition,  the  weak  dependence  on  wave¬ 
length  in  the  wave  phase  speed  leads  to  the  occurence  of  strong  nonlinear  interactions  at 
second-order  in  wave  height.  The  combined  effects  of  weak  dispersion  and  nonlinearity  may 
be  modelled  by  the  Boussinesq  equations  (discussed  below  in  section  2),  which  serve  as  a 
reasonably  accurate  general  purpose  model  for  the  domain  in  question  as  long  as  waves  do 
not  become  too  high. 

Since  the  area  being  studied  here  can  be  quite  extensive  in  comparison  to  the  character¬ 
istic  wavelength  of  the  waves  being  modelled,  the  choice  of  an  efficient  numerical  solution 
technique  is  crucial.  The  Boussinesq  equations  may  be  successfully  solved  by  time-stepping 
techniques  (Wu  and  Wu,  1983;  Rygg,  1988  for  example),  but  these  methods  are  too  ineffi¬ 
cient  for  the  spatial  and  temporal  spans  being  considered  in  applications  to  an  open  coastal 
zone.  Instead,  various  prior  studies  have  found  that  it  is  efficient  to  decompose  the  time- 
dependent  wave  train  into  a  stack  of  frequency  components  by  Fourier  decomposition,  and 
then  solve  for  the  spatial  evolution  of  each  frequency  mode.  This  is  exactly  the  same  kind 
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of  manipulation  that  is  done  in  Foorier  analysis  of  data,  and  is  valid  as  long  as  the  wave 
trmn  satisfies  the  requirements  for  the  transform  to  exist.  In  fidd  applications,  the  Fourier 
decomposition  into  a  finite  set  of  discrete  frequencies  under  the  assumption  of  periodicity 
is  exactly  the  same  as  applying  an  FFT  to  data.  Assumptions  about  stationarity  of  the 
process  are  implicit  and  should  be  kept  in  mind  when  determining  if  the  present  approach 
is  useful. 

The  solution  for  the  spatial  behavior  of  each  frequency  component  is  itself  invdved.  A 
number  of  studies  (Freilich  and  Guza,  1984;  Elgar  and  Guza,  1985  and  1986  and  subsequent 
contributions)  have  considered  only  the  shoreward  propagation  in  one  space  dimension,  and 
have  shown  that  the  restilting  one-dimensional,  first-order  coupled  equations  are  ca¬ 
pable  of  modelling  the  evolution  of  the  spectrum  and  bispectrum  of  choaling  waves,  as 
verified  by  comparison  with  field  measurements.  Liu,  Yoon  and  Kirby  (1985;  referred  to 
here  as  LYK)  alternately  proposed  a  parabolic  model  for  the  shoreward  propagation  of  each 
frequency  component  over  two-dimensional  topogr^hy,  and  verified  the  method  in  com¬ 
parison  to  a  laboratory  experiment  on  the  focussing  of  periodic  long  waves  by  topography. 
The  parabolic  model  approach  has  never  been  used  to  study  a  field  application,  although 
such  work  is  presently  underway  (Freilich  et  al,  1990b). 

In  this  study,  we  derive  a  model  for  the  evolution  of  each  frequency  mode  in  two  di¬ 
mensions,  based  on  the  angular  spectrum  approach.  This  approach  has  been  applied  to 
the  study  of  intermediate  depth  waves  by  Dalrymple  and  Kirby  (1988)  and  Dalrymple  et  al 
(1989).  For  the  case  of  linear  theory,  the  angular  spectrum  is  usually  posed  as  the  contin¬ 
uous  Fourier  transform  of  the  wave  field  in  the  longshore  direction.  Here,  in  keeping  with 
the  mental  framework  associated  with  the  discrete  transforms  being  employed  in  time,  we 
impose  longshore  periodicity  as  well  and  obtain  a  discrete  spectrum  in  longshore  wavenum¬ 
ber.  This  frequency  -  longshore  wavenumber  spectrum  is  subsequently  referred  to  as  the 
discrete  anjpilar  spectrum,  as  referred  to  in  the  title  of  this  paper. 

In  section  2,  we  describe  the  reduction  of  the  Boussinesq  equation  to  a  general  model 
for  the  on-offshore  behavior  of  each  wave  train  resulting  from  Fourier  decomposition  of  time 
and  longshore  direction.  In  this  derivation,  bathymetry  is  assumed  to  be  uniform  in  the 
longshore  direction.  In  section  3,  this  general  form\ilation  is  restricted  to  the  case  of  waves 
propagating  towards  shore,  and  the  angular  spectrum  is  introduced  as  the  slowly  varying 
amplitude  of  each  component  wave.  The  evcdution  equations  are  described.  An  equation 
describing  the  gmn  or  loss  of  energy  in  each  shoaling  wave  mode  is  derived,  and  an  angular 
bispectrum  is  introduced  which  could  be  used  to  study  the  strength  of  coupling  across  a 
range  of  wave  directions.  Several  reductions  of  the  full  model  are  considered  which  would 
reflect  a  narrowing  of  allowed  range  of  wave  directions  in  the  modelled  wave  field. 
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In  section  4,  we  introduce  some  simple  examples  which  illustrate  the  strength  of  interac¬ 
tion  that  may  occur  between  waves  travelling  in  different  directions,  relative  to  interactions 
occuring  in  waves  travelling  in  the  same  direction.  In  particular,  the  characteristic  hexag¬ 
onal  stem  and  saddle  structure  associated  with  solutions  for  short  crested  waves  (S^ur 
and  Finkel,  1985)  is  seen  to  be  related  to  the  presence  of  a  strong  second-harmonic  compo¬ 
nent  resulting  from  the  interaction  of  the  two  obliquely  interacting  primary  waves.  (This 
observation  would  also  follow  from  an  examination  of  the  spectrum  of  a  Genus  2  KP  equar 
tion  solution  (Segur  and  Finkel),  but  this  interpretation  has  not  been  discussed  extensively 
before.) 

The  model  developed  here  is  applicable  both  to  simple  periodic  waves  generated  in  the 
laboratory  environment  and  to  irregular,  “random  ”  waves  in  the  field  environment.  In  this 
paper,  we  concentrate  on  a  model  verification  conducted  from  the  first  viewpoint.  In  section 
5,  the  model  is  tested  against  laboratory  data  obtained  by  Hammack,  Scheffner  and  Segur 
(1990)  who  investigated  the  development  of  a  Mach  stem  arising  during  the  glancing- angle 
reflection  of  a  cnoidal  wave  by  a  vertical  wall.  The  model  developed  here  is  tested  against 
the  laboratory  data  along  with  an  existing  model  based  on  the  parabolic  approximation 
described  in  LYK.  (The  case  of  Mach  reflection  of  a  cnoidal  wave  has  been  studied  using 
the  parabolic  approximation  by  Yoon  and  Liu  (1989);  however,  no  comparison  to  data  was 
provided  in  that  study.)  The  spectral  model  is  shown  to  ^ve  a  more  accurate  representation 
of  the  data  than  the  parabolic  model  over  the  entire  range  of  angles  of  incidence  considered. 


2  The  basic  model 


k 
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VVs  first  establish  the  form  of  a  model  for  waves  in  a  laterally  unrestricted  domain.  A 
Cartesian  coordinate  system  is  adopted  which  has  x  pointed  in  the  onshore  direction  and 
y  pointing  alongshore.  Depth  is  assumed  to  vary  as  h(x)  only.  We  take  as  a  starting  point 
the  variable  depth  Boussinesq  equations  as  given  by  Peregrine  (1967): 


ijt  +  V  •  (hu)  +  (c)V  •  (»7u)  =  0 
u,  +  (Ou .  Vu  +  gVr,  =  (m"){|v(V  •  (hu,))  -  y  V(V  ■  u,)} 


(1) 

(2) 


Here,  rf  is  the  surface  displacement  and  u  is  the  horizontal  wave-induced  velocity.  The 
equations  are  kept  in  dimensional  form;  the  scaling  parameters  e  for  nonlinearity  and  for 
weak  dispersion  are  present  only  schematically  and  will  be  subsequently  dropped.  We  will 
assume  that  either  bottom  slope  or  the  amplitude  of  bottom  features  (as  scaled  by  water 
depth)  are  also  small  and  hence  the  model  will  be  developed  to  leading  order  in  nonlinearity, 
dispersion,  and  domain  inhomogeneity.  This  leads  to  immediate  neglect  of  bottom  slope 
effects  in  the  dispersive  terms  of  (2).  Using  the  linear  portion  of  (1),  we  may  then  write  (2) 
in  the  reduced  form 


u,  -I-  u  •  Vu  -i-  gVrj  +  =  0. 


(3) 


We  now  make  the  following  assumptions.  First,  the  model  will  be  applied  to  time- 
periodic  wave  trains,  where  periodicity  is  in  the  sense  of  either  a  regular  wave  train,  or  of 
a  discrete  FFT  over  a  finite  length  of  sampled  data.  This  assumption  has  been  routinely 
applied  in  the  spectral  sense  in  the  one- directional  shoaling  model  of  Freilich  and  Guza 
(1984)  (see  also  Elgar  and  Guza:  1985,  1986).  Its  interpretation  in  the  regular  wave  case  is 
straightforward,  with  the  wave  being  separated  into  its  harmonic  components  (Rogers  and 
Mei,  1978;  LYK). 

Secondly,  the  wave  field  will  be  assumed  to  be  periodic  in  the  transverse  y  direction. 
This  corresponds  again  to  a  fixed  longshore  wavelength  in  the  regular  wave  case,  or  to 
periodicity  over  a  long  spatial  interval  in  the  spectral  sense. 

The  governing  equations  are  first  split  into  coupled  elliptic  models  for  separate  harmonic 
components.  Following  LYK,  we  write  the  surface  displacement  and  velocity  as 


V-  2^  e  +  c.c. 

(4) 

n=0  ^ 

n=sO  ^ 

(5) 
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Substitution  of  (4)  and  (5)  in  (1)  and  (3)  and  subsequent  elimination  of  the  vdocity  leads 
to  the  following  modd  equation  for  the  rjn  in  the  horizontal  plane: 

n  +  V .  (Gn^nn)  +  [»./.<.]«  =  0;  n  =  1, JV  (6) 

Here,  [n./.t]„  denotes  the  nonlinear  interactions  with  other  discrete  frequency  components 
which  are  sorted  by  means  of  the  rules  for  triad  interactions  applied  to  the  time  dependence. 
The  term  is  given  in  Appendix  A.  The  mode  n  =  0  corresponding  to  the  steady,  wave- 
induced  setdown  is  neglected  since  it  is  at  most  second  order  in  the  largest  wave  amplitudes 
present  (see  LYK).  Also, 


Gnix)  =  gh{x)  ~  |nVh^(x).  (7) 

We  now  apply  a  spectral  transform  to  the  y  dependence  of  the  wavehdd,  assuming 
propagation  is  to  be  considered  in  the  on-ofishore  (^x)  direction.  We  consider  here  the  case 
of  an  unbounded  lateral  domain  and  a  wavefield  which  is  periodic  over  the  basic  interval  L. 
We  then  represent  qn(*,y)  as 

-)»(*.»)=  E  (8) 

fn=— M 


where 

Ao  =  y .  (9) 

Substituting  (8)  in  (6)  and  neglecting  x  -  derivatives  of  small  terms  in  Gn  then  leads  to  a 
set  of  coupled  second-order  ODE’s  for  the  q",  ^ven  by 


=  0; 


(10) 


where  [n./.t.]^  now  represents  triad  interactions  satisfying  resonance  conditions  in  i  and  y. 
The  nonlinear  term  is  again  given  in  Appendix  A.  Here,  k  is  the  wavenumber  determined 
by  the  lowest  order  dispersion  relation 

=  gk^h.  (11) 


Also, 

W)’  =  -  m’Al  (W) 

For  fixed  n,  ib,  Aq,  large  values  of  m  will  make  ima^nary,  which  corresponds  to  modes 
which  are  exponential  rather  than  osdllatory  in  x  in  the  linear  approximation.  In  the  linear 
case,  the  presence  of  these  modes  in  the  initial  conditions  would  be  interpreted  in  the  same 
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light  as  the  presence  of  evanescent  modes  in  the  general  wavemaker  problem  (see  Dalrymple 
and  Kirby,  1988).  However,  the  interpretation  in  the  case  of  possible  nonlinear  forcing  of 
the  offshore  portion  of  trapped  modes  in  the  nearshore  region  is  non-trivial  and  will  need 
to  be  considered  carefully  in  applications  where  the  inclusion  of  this  effect  is  desired.  In 
addition,  nonlinearity  could  force  the  propagation  of  modes  that  would  not  be  present  in  a 
linearized  wave  field,  and  which  could  affect  a  detailed  representation  of  an  individual  wave. 
At  present,  the  range  of  M  at  each  value  of  n  may  be  restricted  to  Af„  <  nk/Xo  in  order  to 
eliminate  forcing  of  these  modes  arbitrarily. 
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3  Shoaling  waves 


I 


The  model  developed  in  the  previons  section  allows  for  the  onshore  and  offshore  propagation 
of  the  directional  spectrum  components.  For  the  remainder  of  the  present  paper,  attention 
is  restricted  to  waves  propagating  onshore,  or  in  the  positive  sense  with  respect  to  the  x 
coordinate. 


3.1  The  Model  Equations 

Based  on  the  linear,  nondispersive  portion  of  the  model  (10),  we  assume  that  the  incident 
wave  may  be  written  in  the  form 

•;?(*)  =  /trcxy*/”- (13) 


where  it  is  assumed  that  the  x  dependence  of  A,  k  and  7  is  on  a  slow  scale  of  0(e),  and 
where 


n 


(14) 


(where  the  positive  root  is  taken).  The  amplitudes  A  represent  the  discrete  angular  spectrum 
being  considered  here,  and  are  allowed  to  vary  owing  to  refraction,  shoaling,  dispersion  and 
nonlinear  interaction.  (It  would  be  possible  to  absorb  shoalin?  and  refraction  effects  by  the 
use  of  the  usual  linear  refraction  formulae;  this  step  is  not  taken  here.)  Substitution  of  (13) 
in  (10)  leads  to  the  spectral  model  for  incident  waves,  given  by 


+ 


ink 

8h 


+  2 


(**7^): 


1 . 


2kir^  - 


( /=1  p=/»i  i=l  p=Pi  ) 


n  = 


Tn  — 


0; 

(15) 


Here,  (•)*  denotes  the  complex  conjugate.  The  limits  of  summation  P\  —  P4  are  pven  by 

Px  =  max(-M/,-Af„_i  +  m) 

P2  =  mia(Mi,  M„_/  +  m) 

P3  =  max(-Mi,-Mn+i  -  m) 

F4  =  min(Af/,Afn+J  -  ”*)  (16) 


The  interaction  coefficients  I  and  J  are  given  by 


cr- 


(/^  -Kn  - 


](17) 

(18) 
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The  phase  arguments  6  and  T  represent  the  basic  mismatch  in  the  z  direction  of  the  triads 
chosen  based  on  perfect  matching  in  y  and  t.  Generally,  the  only  components  which  expe¬ 
rience  complete  resonance  in  the  long  wave  limit  must  have  parallel  propagation  directions; 
all  obliquely  interacting  components  are  somewhat  detuned.  The  phase  arguments  are  given 
by 


ey  =  ltTf  +  (n-l)t7r-?-"*7r  (19) 

T”/  =  (M) 


The  spectral  model  (15)  is  a  set  of  coupled  first  order  ODE’s  which  are  solvable  by 
standard  techniques.  Results  presented  in  sections  4  and  5  were  obtained  using  a  standard 
4th-order  Runge-Kutta  scheme  with  fixed  step  size  and  no  error  checking.  Presently,  exten¬ 
sions  of  the  model  using  error  control  and  adaptive  step  size  are  being  tested  and  wiR  likdy 
be  the  vehicle  for  further  field  testing  of  the  model  beyond  the  scope  of  the  present  study. 

The  model  developed  here  allows  for  coupling  and  energy  transfer  between  spectral 
modes  with  different  directions  and  frequencies  as  the  wave  field  evolves  in  the  shoreward 
direction.  For  the  case  when  modes  with  real  wavenumbers  7^  are  retained,  an  energy 
equation  for  each  mode  may  be  constructed  by  multiplying  (15)  by  ^4^*  and  adding  the 
resulting  equation  to  its  complex  conjugate,  to  get 


nk  f"”!  ^ 

-JjELC 

U=lp=F, 


7;rMril  +  wi* 


l=l  p=P3 


=  0 

(21) 


The  quantity  B  introduced  here  is  given  by 

Kf  =  (K) 


and  may  be  thought  of  as  being  a  discrete  representation  of  a  directional  bispectrum.  Anal¬ 
ysis  of  the  magnitude  of  this  term  should  provide  a  measure  of  the  strength  of  nonlinear 
coupling  across  a  range  of  directions  in  a  field  spectrum.  This  aspect  of  the  analysis  is  not 
pursued  further  here.  In  (21),  the  first  appearance  of  the  directional  bispectrum  corresponds 
to  coupling  between  the  present  mode  n  and  modes  of  lower  frequency;  the  present  wave 
takes  the  role  of  a  sum  wave  in  a  higher-harmonic  excitation.  The  second  appearance  of 
bispectral  coupling  corresponds  to  cases  where  the  present  mode  has  either  the  lowest  or 
middle  frequency  in  a  triad. 
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3.2  Approximations  for  small  angles  of  incidence 

Parabolic  modeling  schemes  have  drawn  a  reasonable  amount  of  interest  as  a  way  of  cmn- 
puting  wave  fields  in  the  coastal  environment.  Since  these  models  presuppose  a  narrow 
range  of  wave  propagation  directions  relative  to  the  full  180°  band  allowed  here,  it  is  of 
some  interest  to  see  how  the  present  model  could  be  reduced  in  scope  to  correspond  to  the 
assumptions  made  in  developing  the  parabolic  approximation.  (Note  that  the  parabolic 
model  computations  presented  in  section  5  were  made  using  the  standaird  model  in  physical 
space,  and  are  not  an  application  of  the  results  derived  here). 

We  consider  first  the  effect  of  assuming  that  the  spread  of  wave  directions  about  a  central 
angle  is  small.  In  particular,  we  choose  the  central  angle  as  9  =  0°,  coincident  with  the  z 
axis.  The  effect  of  having  waves  propagating  at  close  to  normal  incidence  is  to  lengthen 
the  longshore  interval  needed  in  order  to  assume  spatial  periodicity.  Thus,  Ao/Jb  in  (14'20) 
would  be  effectively  small.  Further  restricting  attention  to  small  an^e  deviations  from  the 
normal  direction  limits  consideration  to  relatively  small  values  of  m.  Consequently,  the  7 
may  be  approximated  according  to 

which  is  consistent  with  the  usual  smaU-angle  form  of  the  parabolic  approximation.  The 
resulting  modifications  to  the  phase  arguments  and  interaction  coefficients  are  ^ven  by  (to 
the  same  level  of  approximation) 


(24) 

(25) 


where  the  second  form  in  equation  (24)  results  from  the  need  to  avoid  a  zero  arising  at 
finite  angles  in  the  first  approximation.  This  model  then  corresponds  to  the  case  where 
angular  variations  between  interacting  wave  components  are  small,  but  where  no  direct 
correspondence  is  made  between  the  smallness  of  this  dependence  and  the  other  small 
parameters  in  the  problem.  A  further  reduction  in  complexity  would  occur  if  we  were 
to  make  the  full  set  of  restrictions  appropriate  to  the  parabolic  approximation  of  the  KP 
equation  (KadomtseV'Petviashvili,  1970),  which  is  available  in  variable  depth  form  from 
LYK.  Carrying  the  ordering  restrictions  to  their  fullest  extreme,  we  would  conclude  that, 
since  the  nonlinear  terms  are  already  of  small  size,  only  0(1)  information  should  be  retained 


9 


in  their  coefficients.  This  further  assumption  reduces  the  coefficient  I  to 

C/  =  3  (M) 

for  all  components  of  the  wave  train. 

To  date,  the  various  forms  of  small-angle  approximations  have  not  been  fully  tested 
for  their  ability  to  model  the  general  (unrestricted  angle)  evolution  of  a  wave  train.  The 
existence  of  various  approximations  of  the  full  governing  equations  leads  to  the  question  of 
whether  the  approximate  interaction  coefficients  may  be  applied  in  a  directional  wave  calcu¬ 
lation  with  finite  directional  bandwidth.  For  example,  if  the  presence  of  detuning  effectively 
nullifies  the  nonlinear  interaction  in  a  particular  triad  before  the  point  is  reached  where  the 
interaction  coefficient  deviates  significantly  from  one  of  its  small- angle  forms,  then  the  small- 
angle  forms  may  be  retained  in  the  full  governing  equations  and  sums  over  transverse  mode 
numbers  may  be  limited  to  the  bandwidth  contributing  to  resonance.  Support  for  this  hy¬ 
pothesis  is  fragmented,  but  comes  from  both  numerical  and  experimental  results  indicating 
the  large  angular  range  for  which  the  KP  equation  produces  reasonably  valid  results.  For 
example,  Kirby  et  al  (1988)  have  demonstrated  that  the  KP  model  accurately  describes 
reflection  of  solitary  waves  by  a  wall,  in  comparison  with  data  obtained  by  Melville(1980) 
and  with  numerical  solutions  of  the  full  Boussinesq  theory  by  Funakoshi  (1981).  However, 
these  computations  are  limited  to  the  near  field  of  the  start  of  the  reflecting  wall  and  do  not 
adequately  describe  far  field  behavior.  The  reflection  pattern  is  reproduced  for  angles  of 
incidence  well  beyond  the  critical  angle  for  producing  a  Mach  stem,  as  defined  theoretically 
by  Miles  (1977),  at  least  for  short  distances  in  a  wave  channel.  Similarly,  Hammack  et  al 
(1989)  have  demonstrated  that  the  theory  of  Segur  and  Finkel  (1985)  for  KP  solutions  of 
biperiodic  (short-crested)  waves  may  be  applied  accurately  to  wave  trains  having  essentially 
no  apparent  scale  distinction  in  the  underlying  x  —  y  parallelogram. 

The  determination  of  the  validity  of  the  KP  approximation  relative  to  the  full  Boussi¬ 
nesq  equation  will  require  careful  studies  with  models  of  the  present  type  and  with  time- 
dependent  studies  for  transient  phenomena.  The  results  in  section  4  and  5  are  a  start 
towards  such  a  study,  and  indicate  that  the  range  of  validity  of  the  smaU  angle  approx¬ 
imations  may  be  quite  limited,  when  actual  computational  evolution  of  the  wave  field  is 
con»dered. 
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4  Analysis  for  a  simplified  wave  field 


Before  moving  to  model  -  data  comparison,  we  consider  a  simple  example  which  illustrates 
the  strength  of  possible  triad  interactions  occuring  between  waves  travelling  at  different 
directions,  relative  to  interactions  occuring  between  waves  travelling  in  the  same  direction. 
For  this  example,  we  consider  the  case  of  the  oblique  interaction  of  two  synchronous  n  =  1 
modes  oriented  at  ±0  with  respect  to  x  (m  =  ±1,  representing  a  linear  short-crested  wave 
train),  which  interact  to  form  2n  =  2  modes  oriented  with  z  (m  =  0,  the  sum  interaction) 
and  at  with  respect  to  z  (m  =  ±2,  the  self-interaction  of  each  linear  mode).  Water 
depth  is  assumed  to  be  constant.  We  further  assume  that  the  oblique  modes  are  equal  in 
amplitude  (Aj  =  =  -dj^)  and  obtmn  a  reduced  set  of  three  equations.  For  this 

reduced  case  involving  only  one  wave-angle  0,  it  is  convenient  to  express  Aq  in  terms  of  k 
and  6.  We  take  Aq  =  ksind  for  n  =  l,m  =  ±1.  We  further  nondlmensionalize  all  id's  by 
i4i(0)  and  z  by  and  obtain  the  equations 


7^1^  - 

-i- ^€(  -1- 3A}M^)  =  0 

(27) 

(28) 

’rAl^-^fi'^Al  +  j€{A\f  =  0 

(29) 

where 

7 

=  cosd 

(30) 

J 

=  1-1-  cos  d(6  —  4  cos  d) 

(31) 

I 

=  3  —  5  sin’ d  4-  2  sin^  d 

(32) 

e 

=  2(1 -cosd) 

(33) 

Equations  (27-29)  admit  a  corresponding  conservation  law 

2C06d(|y4i| 

’  +  Mai*)  +  yMjl’  =  constant 

(34) 

For  the  case  of  a  small  angle  approximation  consistent  with  the  information  in  equations 
(24)  and  (25),  we  obtain  the  approximations  coed  -►  1  -  Agin^d,  J  3  sin’tf,  I 
3/(1  -i-  (5/3)sin*d)  and  0  -♦  sin*d.  For  the  parabolic  model  used  in  section  5,  I  and  J 
further  reduce  to  I  =  J  =  3.  The  resulting  evolution  is  expected  to  be  dependent  on 
=  {kh)‘^,€  =  i4i(0)/h  and  9.  We  are  interested  in  testing  the  apprmcimate  models  as 
0  increases  for  fixed  values  of  and  c,  and  in  studying  the  resulting  interactions.  Figure 
1  gives  plots  of  I{0),  J{0)  and  6(d)  and  their  small  an^e  approximations  over  the  range 
0"  <  d  <  90*. 


11 


4.1  Steady  solutions 


The  three  coupled  evolution  equations  may  be  8(^ved  simply  for  the  case  of  a  steady  wave 
of  permanent  form.  Letting 

A\  =  aexp(^) 

=  /3€xp{|t/i2i)  (35) 

M  =  PexpCj  — ) 

leads  to  elimination  of  the  dispersive  terms  and  ^ves  the  revised  equations 


70,  -1-  -1-  =  0 

4  4 

(36) 

0,  ^  ~o  =  0. 

(37) 

7Px  +  ~Q^e"*^»*  =  0. 

4 

(38) 

where 

(39) 

> 

II 

(40) 

We  then  seek  steady  solutions  of  the  form 

Q  =  oe‘^‘*;  0  =  6e‘*»';  p  =  ce*^* 

(41) 

where  u,6,c  are  real  amplitudes  and  the  6's  are  phases.  Defining  t  =  i4i(0)/h,  we  then  take 
0=1.  Substitution  of  (41)  in  (36-38)  and  use  of  the  properties  of  the  amplitudes  leads  to 
the  conditions 


62  —  2Si  -f  Ai  =  0, 

^3  —  2Si  -f-  A2  =  0. 


-  <Jb  3<c  -  c/  3c 

47  47  2b  47c 

The  solution  is  then  given  by 

Z(b 

27[/f-^2A'6]’ 

where  A'  =  A2  —  Ai,  and  a  cubic  equation  for  5  given  by 

/n  r  I  /.Tt  I  47A1A  m.2  .  rn^  r/ a  a /mi. 


(2yA')6^  +  (c/y  + 


+  ^)6’  l27/(Ar  -  A016  -  =  0. 


(42) 

(43) 

(44) 

(45) 

(46) 
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In  the  parameter  ranges  of  interest  (0  <  (/**»«)  <  0-5,  0“  <  <  90®),  the  roots  of  (46)  are 
typically  real  and  distinct.  However,  only  one  root  leads  to  a  positive  value  for  c  and  hence 
a  co-directional,  in-phase  harmonic  for  each  individual  primary  wave.  This  root  is  taken  to 
be  the  physically  plausible  solution. 

As  $  -*  90®,  the  solutions  of  (46)  become  difficult  to  evaluate  numerically.  0  =  90® 
represents  primary  waves  travelling  in  opposite  directions;  these  waves  do  not  interact  at 
leading  order  (Rogers  and  Mei,  1978)  and  hence  b  drops  to  0  in  this  limit.  Therefore, 
solutions  for  c  are  obtained  directly  from  a  truncated,  two  equation  model,  which  gives 


where  the  positive  root  for  c  is  chosen.  This  result  is  utilized  if  either  of  the  coefficients  of 
the  zeroth  and  first  order  terms  in  (46)  fall  below  10~^  in  absolute  value.  This  limit  does 
not  have  a  great  deal  of  practical  importance  in  the  present  context. 

Figure  2  presents  solutions  for  the  harmonic  amplitudes  b  and  c  for  =  0.36  and 
a  range  of  f  and  0  values.  The  value  of  is  chosen  to  correspond  approximately  to 
the  tests  in  section  5.  The  most  striking  feature  of  the  result  is  that  the  amplitude  b  of 
the  sum-wave  component  at  normal  incidence  exceeds  the  amplitude  c  of  the  inline  second 
harmonics  of  each  oblique  primary  wave  for  sufficiently  small  angles  of  incidence.  This  result 
is  not  immediately  intuitive,  since  any  triad  based  on  oblique  components  would  necessarily 
be  more  detuned  than  the  corresponding  triad  based  on  second-harmonic  resonance.  See 
Figure  3  for  the  geometry  of  several  hypothetical  detuned  triads  whose  component  lengths 
are  determined  from  the  lowest  order  relation  (11).  The  results  in  Figure  2  indicate  that  the 
result  of  the  detuned  triad  interaction  actually  dominates  the  higher  harmonic  contents  for 
angles  of  intersection  between  the  two  primary  waves  as  large  as  50®  —  60®.  Of  course,  the 
conclusion  that  detuned  interactions  are  dominating  over  more  closely  tuned  interactions 
is  an  artifact,  in  this  case,  of  the  assumption  that  the  component  wavenumbers  are  given 
by  the  linear  non-dispersive  theory.  In  a  permanent-form  wave,  nonlinearity  distorts  all 
component  wavenumbers  in  such  a  way  that  components  satisfying  resonance  conditions 
in  frequency  and  longshore  wavenumber  also  satisfy  resonance  conditions  in  on-offshore 
wavenumber.  If  this  were  not  true,  the  components  would  have  different  wavelengths  in 
the  T  direction  and  permanent  form  would  not  be  retained.  The  short  crested  wave  field 
computed  here  could  be  regarded  as  being  a  truncated  form  of  the  Genus  2  KP  solution  of 
Segur  and  Finkel  (1985)  in  the  event  that  the  intersection  angle  of  the  two  primary  waves 
approaches  zero  (and  KP  scaling  restrictions  are  satisfied).  Conversely,  reta'ning  higher 
hsT'^nonlcs  and  determining  solutions  of  the  full  problem  would  provide  the  large  angle 
extension  of  the  KP  sdution. 


13 


Figure  4  shows  a  plot  of  two  periods  of  the  waveheld  for  the  case  =  .36,  c  =  .1,7  = 
0.88,  which  corresponds  reasonably  closely  to  the  test  condition  CR38  discussed  in  section 
5.  For  this  case,  6  =  0.15522  and  c  =  0.18204;  the  in-line  and  center  angle  second  harmonic 
components  are  nearly  balanced.  The  resulting  wave  field  exhibits  the  hexagonal  form 
characteristic  of  short-crested  shallow  water  waves,  as  discussed  in  Segur  and  Finkel  (1985) 
and  Hammack  et  al  (1989).  Wavefields  dominated  by  the  center  angle  component  b  give 
the  appearance  of  being  long  crested  with  phase  jumps  confined  to  a  more  narrow  saddle 
region;  see  Figure  8  below.  In  contrast,  wavefields  dominated  by  inline  harmonics  would 
instead  appear  more  diamond  shaped,  as  is  the  case  in  linear  theory;  see  Figure  9  below. 

In  order  to  get  an  idea  of  how  much  the  solution  of  a  small-angle  approximation  would 
deviate  from  the  predictions  of  the  full  theory,  we  have  recomputed  the  solution  for  h  and  c 
for  the  parameters  used  to  obtain  the  solution  in  Figure  2.  The  solution  is  obtained  using  the 
small  angle  values  for  7  and  G  and  the  parabolized  values  of  7  =  7  =  3,  which  represents  the 
most  extreme  reduction  of  the  theory  in  terms  of  directional  resolution,  barring  reduction 
to  one-dimension.  Results  are  shown  in  Figure  5.  The  figure  is  similar  to  Figure  2  except 
that  the  plots  of  6  and  c  seem  to  be  stretched  out  along  the  6  axis;  the  waves  behave  as  if  the 
ngle  between  the  primary  components  is  less  than  it  actually  is.  This  is  due  to  the  fact  that 
the  small  angle  estimate  of  the  value  of  7  drops  off  less  quickly  than  does  the  true  value, 
and  to  the  fact  that  the  use  of  /  =  7  =  3  causes  the  nonlinear  interactions  to  behave  more 
like  they  would  at  vanishingly  small  angle  of  incidence  than  would  be  expected.  Note  that 
tests  in  .‘iection  5  cover  a  corresponding  range  of  9's  up  to  about  45®;  deviations  between 
the  full  solution  and  the  parabolic  form  are  significant  in  this  range. 

4.2  Unsteady  evolution 

The  results  for  steady  periodic  waves  indicate  that  nonlinear  interactions  can  distort  the 
underlying  components  of  the  wave  train  to  the  extent  that  components  that  would  be 
expected  to  be  only  weakly  tuned  based  on  linear  parameters  are  actually  dominant  parts  of 
the  wave.  In  this  section,  we  consider  instead  the  unsteady  evolution  of  similar  wave  trains, 
in  which  second  harmonics  are  allowed  to  grow  from  zero  amplitude.  We  are  interested  in 
the  competition  between  the  growing  second-harmonic  amplitudes  and  the  dependence  of 
the  evolution  on  the  intersection  an(^e  between  the  primary  components.  The  numerical 
experiments  described  here  are  similar  in  spirit  to  the  studies  of  second-harmonic  generation 
in  one  space  dimension  performed  by  Mei  and  Unluata  (1972)  and  Boczar-Karakiewicz 
(1972),  with  the  exception  that  the  two-dimensionality  of  the  problem  introduces  a  set  of 
competing  second  harmonics.  Depending  on  the  angle  of  incidence  of  each  primary  with 
respect  to  the  z  axis,  one  of  the  two  possible  second-harmonic  components  may  become 
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dominant  in  the  hannonic  generation  process.  The  dominance  corresponds  closely  to  which 
of  the  components  6  or  c  takes  on  the  larger  valne  in  the  steady  solution;  larger  values  of 
b  indicate  that  the  sum- wave  interaction  will  become  dominant,  larger  values  of  c  indicate 
that  the  in-line  components  will  be  dominant,  and  values  of  b  and  c  which  are  close  in 
magnitude  indicate  that  the  two  harmonics  will  evolve  in  a  complex  fashion  with  neither 
being  dominant.  An  iUustration  of  these  effects  is  given  in  Figure  6,  which  shows  the  pattern 
of  second  harmonic  evolution  for  the  case  of  =  0.36,  e  =  0.1  over  an  x  distance  of  100 
water  depths.  The  pattern  shifts  from  one  where  the  sum-wave  interaction  dominates  the 
wave  field  evolution  at  small  angles,  to  one  where  the  in-line  components  dominate  and  the 
sum  wave  is  only  of  minor  importance  at  large  angles. 

The  complexity  of  the  wave  field  evolution  at  intermediate  angles  is  of  interest,  and 
sensitive  laboratory  experiments  on  this  topic  could  provide  a  way  of  discriminating  between 
the  detailed  accuracy  of  various  future  versions  of  schemes  like  the  one  developed  here.  The 
relevance  of  these  results  to  natural  events  is  in  question,  though.  Recent  computational 
tests  using  a  one-dimensional  model  (Elgar  et  al,  1990)  have  shown  that  an  evolving  triad 
(or  degenerate  triad  involving  a  primary  component  and  its  second  harmonic)  is  unstable  to 
broad-banded  noise  and  would  thus  not  be  expected  to  be  manifested  in  the  natural  wave 
field  which  has  arrived  after  travelling  some  distance.  This  result  does  not  invalidate  the 
use  of  second-harmonic  generation  tests  in  the  laboratory  as  a  means  of  studying  model 
accuracy. 

Finally,  use  of  the  full  model  (15)  with  a  large  number  N  of  components  retained  would 
allow  for  the  generation  of  a  more  complex  pattern  of  higher  harmonics;  the  resulting  wave 
fields  are  not  qualitatively  different  from  the  ones  shown  here  and  still  exhibit  recurrence 
or  quasi-recurrence. 
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5  Comparison  with  Data  and  a  Parabolic  Model  Approxi¬ 
mation 

In  order  to  verify  the  basic  computational  model  provided  by  (15)  and  to  test  whetha  the 
present  angular  spectrum  provides  a  more  accurate  representation  of  the  wave  field  relative 
to  earlier  parabolic  models  (LYK;  Yoon  and  Liu,  1989),  we  have  compared  model  predictions 
to  laboratory  data  obtained  by  Hammack  et  al  (1990)  for  the  case  of  glancing,  or  Mach, 
reflection  of  a  cnoidal  wave  by  a  vertical  wall.  The  experimental  tests  were  conducted  using 
the  directional  wave  maker  at  the  Coastal  Engineering  Research  Center,  Vicksburg,  MS.  A 
prior  use  of  this  facility  to  study  the  properties  of  intersecting  cncadal  waves  is  described  in 
Hammack  et  al  (1989),  referred  to  here  as  HSS. 

5.1  Layout  of  Experimental  Facility 

For  the  tests  considered  here,  the  wave  basin  was  operated  with  a  water  depth  of  20cm  in 
a  constant  depth  region  extending  12.55m  in  front  of  the  wavemaker,  after  which  a  beach 
with  1:30  slope  provided  an  efficient  wave  absorber  giving  little  reflection.  The  basin  floor 
was  leveled  to  a  tolerance  of  0.01  ft  in  order  to  remove  some  of  the  spatial  irregularity  of  the 
waves  reported  in  HSS  resulting  from  local  refraction  effects.  For  the  Mach  stem  tests,  two 
parallel  false  walls  were  installed  perpendicular  to  the  wavemaker  axis  in  order  to  provide 
a  closed  channel.  The  channel  walls  were  situated  13.26m  apart,  which  fixes  the  width  of 
the  numerical  domain  to  be  considered. 

Instrumentation  and  data  aquisition  are  described  in  HSS,  and  readers  are  referred  there 
for  greater  detail.  In  the  present  tests,  an  array  of  18  wave  gages  were  installed  in  the  basin 
as  shown  in  Figure  7.  In  this  figure,  the  directional  wavemaker  occupies  the  y-axis,  and 
the  channel  sidewalls  lie  along  the  x-axis  and  the  line  y  =  -13.26m.  The  positions  of 
the  gages  are  given  in  Table  1.  The  linear  gage  array  13-9-8-7-6-4-3-2-1 
provides  a  transect  perpendicular  to  the  reflecting  sidewall  which  allows  a  determination  cA 
the  width  and  structure  of  the  Mach  stem  and  additional  crests  in  the  reflection  pattern. 
This  array  is  used  to  provide  most  of  the  information  described  below.  An  additional  array 
18-17-16-15-14-13-5  provides  measurements  of  the  evolution  of  the  reflected 
stem  wave  along  the  wall.  For  each  gage,  data  consists  of  a  time  series  of  1250  points  with 
a  sampling  rate  of  2bHz. 

5.2  Specification  of  the  Incident  Wave 

The  generation  of  oblique  cnoidal  waves  using  the  directional  wavemaker  has  been  described 
in  HSS.  In  the  present  laboratory  tests,  waves  were  initially  specified  as  having  a  wavelength 
of  2m  and  a  crest  elevation  4cm  above  mean  water  level.  The  algorithms  given  by  Goring  and 


Raichlen  (1980)  were  used  to  generate  a  time  series  of  paddle  displacement  corresponding  to 
one-dimensional  generation.  Oblique  waves  were  then  generated  by  phase  lag^ng  adjacent 
paddles.  The  relation  between  paddle  phase  shift  angle  0  and  directed  wave  an^e  a  is  given 
by 

0  =  arcsin(QtL/360Vr)  (48) 


where  L  is  the  wave  length  and  W  =  45cto  is  the  individual  paddle  width.  Tests  were 
conducted  for  six  paddle  phase  lags,  and  are  denoted  CRxx0204,  where  xx  denotes  paddle 
phase  lag.  Note  that  (48)  corrects  a  typographic  error  appearing  in  HSS.  Table  1  gives  a 
list  of  a  and  0  values  for  the  six  tests. 

Prior  to  running  the  Mach  reflection  experiments,  the  nominally  4cm  high  wave  was 
generated  in  the  normally  incident  direction  (traveling  parallel  to  sidewalls)  in  order  to 
study  its  characteristics.  It  was  found  that  the  wave  actually  had  a  crest  elevation  close 
to  3.3cm  above  mean  water  level.  This  value  was  used  to  specify  incident  waves  in  the 
numerical  computation,  along  with  a  wave  period  of  1.478s  as  specified  by  KdV  cnoidal 
wave  theory.  We  note  that  it  is  impossible  to  tell  whether  the  value  of  3.3cm  was  invariant 


under  changes  of  angle  of  incidence  in  the  laboratory  experiment.  Variation  of  this  quantity 
would  add  an  untraceable  source  of  error  in  model-data  comparison. 

It  has  been  shown  (Kirby,  1990)  that  spectral  models  of  the  type  described  here  (with 
lowest  order  non-dispersive  results  used  to  develop  dispersive  and  nonlinear  terms)  have 
one-dimensional  solutions  that  correspond  exactly  to  solutions  of  the  modified  KdV-type 


equation 


In  order  to  account  for  this  factor,  input  data  for  the  computational  models  w'as  specified 
as  follows.  First,  KdV  theory  was  used  to  determine  a  waveheight  and  wave  period  cor¬ 
responding  to  the  wave  described  above.  Then,  a  Fourier  series  solution  of  the  modified 
KdV  equation  (49)  was  performed  using  that  waveheight  and  period  to  specify  the  solution. 
This  solution  led  to  a  predicted  wavelength  of  1.9977m,  which  is  only  slightly  less  than  the 
KdV  prediction  of  2m.  (Recall,  however,  that  the  2m  value  was  for  a  crest  elevation  of  4cm 
rather  than  3.3cm).  This  wavelength  was  then  used  in  (48)  along  with  the  phase  lag  angle 
a  in  order  to  determine  the  directed  angle  0  for  each  test.  A  table  of  resulting  angles  is 
given  in  Table  2. 

The  computed  wave  angle  and  the  Fourier  coefficients  derived  from  the  solution  of  (49) 
were  used  to  compute  the  surface  displacement  for  oblique  cnoidal  waves  along  the  boundary 
X  =  Om  corresponding  to  the  wavemaker.  For  parabolic  model  calculations,  information 
in  the  range  -  13.26m  <  y  <  Om  was  used  to  start  the  computation.  The  model  was  run 
with  reflective  sidewalls  at  y  =  0,- 13.26m.  For  the  spectral  modd,  a  periodic  interval 
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was  constructed  by  using  a  mirror  image  about  y  =  Om;  the  computational  domain  thus 
corresponds  to  the  region  —13.26m  <  y  <  13.26m.  The  computed  waveform  was  then 
Fourier  transformed  over  this  interval  and  the  resulting  frequency-wavenumber  spectrum 
was  fed  into  the  spectral  model.  (It  is  noted  that  the  problem  as  stated  could  be  handled 
directly  by  means  of  a  cosine  transform  over  the  true  model  domain;  this  wais  not  done 
because  of  the  reprogramming  of  the  basic  model  that  would  have  been  required.) 

For  the  examples  shown  here,  the  parabolic  model  was  run  with  a  grid  spacing  Ax  = 
Ay  =  0.0625m,  and  N  =  9  frequency  components  were  used.  Tests  were  performed  for 
two  cases  for  both  half  the  grid  spaicing  and  twice  the  number  of  harmonics  to  insure  that 
convergence  was  adequate  for  the  parameters  finally  used.  In  the  spectral  model,  we  used 
N  =  9  and  Af  =  64.  The  large  value  of  M  insures  that  all  freely  propagating  modes  of 
the  solution  aoe  retained  for  the  highest  harmonic  considered.  The  forward  grid  space  step 
was  also  Az  =  0.0625m.  Again,  these  parameters  were  found  to  give  sufficient  convergence 
of  solutions  when  compared  to  runs  with  smaller  grid  steps  and  more  retained  frequency 
components. 

5.3  Results 

Results  of  model  calculations  are  presented  here  in  two  forms;  gray  level  contour  plots  of  in¬ 
stantaneous  surface  over  the  model  basin,  and  time  series  computed  from  model  output  and 
compared  to  experimental  data.  The  gray  level  plots  are  actually  of  the  quantity  -drjidx, 
and  the  pictures  thus  mimic  the  visual  image  that  would  be  obtained  in  an  overhead  pho¬ 
tograph  resulting  from  lighting  at  a  low  angle  from  the  direction  of  the  wavemaker.  (This 
is  similar  to  the  photographic  arrangement  in  HSS).  In  all  cases,  the  gray  level  plots  from 
parabolic  model  computations  are  similar  to  those  from  spectral  model  computations,  and 
only  spectral  model  results  are  shown.  Figure  8  presents  results  for  the  test  CR150204, 
which  clearly  shows  the  evolution  of  a  wide  Mach  stem  wave  along  the  reflecting  boundary. 
In  contrast,  Figure  9  shows  the  other  extreme  example  of  test  CR580204,  where  the  angle  of 
incidence  is  about  45°  and  the  reflection  pattern  is  regular  (i.e.,  nearly  linear  superposition.) 

In  order  to  compare  time  series  from  experiment  and  model  calculations,  the  following 
procedure  is  used.  First,  after  allowing  some  time  for  nearly  periodic  motion  to  be  estab¬ 
lished  in  the  experiment,  a  reference  time  to  is  established  at  the  position  of  a  wave  crest  at 
gage  13  (z  =  11m,  y  =  -0.07m).  This  start  time  is  used  for  all  other  gsiges  as  well,  in  order 
to  maintain  synchronization.  Then,  a  start  time  is  also  established  for  the  model-predicted 
time  series  by  identifying  a  crest  at  model  gage  13.  Synchronization  between  model  and 
data  time  series  is  thus  based  on  correlating  the  series  at  a  gage  location,  rather  than  (more 
correctly)  correlating  at  the  wavemaker.  This  step  is  necessary  since  the  wavemaker  contnd 
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is  not  available  and  since  absolute  time  in  the  model-constructed  time  series  is  arbitrary. 

After  determining  the  synchronous  start  times  for  model  and  data  at  gage  13,  time  series 
were  plotted  for  the  transect  along  the  wall  (  gages  18-17-16-15- 14-13-5)  and  perpendicular 
to  the  wall  (gages  13-9-8-7-6-4-3-2-1).  Aside  from  the  presence  of  experimental  errors,  a 
correct  model  result  would  be  indicated  by  complete  agreement  between  model  and  data 
time  series  at  each  gage.  As  examples,  two  periods  of  time  series  for  the  perpendicular 
transect  are  shown  in  Figure  10  for  the  spectral  model  predictions  for  CR150204.  The 
plot  shows  a  stem  developed  near  the  wall  in  the  area  spanned  by  gages  13-9-8-7,  beyond 
which  there  is  a  clear  phase  lead  at  each  subsequent  gage,  indicating  the  approach  of  the 
incident  wave  at  a  small  angle  of  incidence.  In  this  plot,  model  predictions  are  indicated 
by  solid  lines  and  experimental  data  by  dashed  lines.  At  the  opposite  extreme.  Figure  11 
shows  spectral  model  predictions  and  data  for  the  case  CR580204.  This  figure  indicates  the 
structure  of  a  short  crested  wave  field,  with  one  complete  diamond  over  the  range  of  gages 
13-9-8-7-6-4-3,  whereas  the  signal  in  gages  3-2-1  indicating  a  strong  progressive  phase  lag, 
corresponding  to  a  wave  travelling  away  from  the  wall  at  close  to  45°.  The  structure  of 
the  plots  in  Figures  10  and  11  may  be  further  clarified  by  comparing  them  to  the  contour 
plots  in  Figures  8  and  9,  respectively. 

In  order  to  quantify  the  comparison  between  model  predictions  and  measurements,  a 
standard  correlation  coefficient  p  was  calculated  for  each  time  series  on  the  perpendicular 
transect.  For  i= corresponding  to  the  discrete  time  series  over  two  wave  periods  and 
j  =  1,  ...,9  corresponding  to  the  9  gages,  we  first  define  a  correlation  for  each  gage  by 


Pj  = 


C\2j 


(50) 


where  subscripts  1  and  2  refer  to  data  and  model,  respectively,  and  where 


1  / 

=  7  5^»hj(»)»?2j(*) 

^  i=i 

=  \  7  5I’?»j(*>i,>(«);  /=i,2 


(51) 

(52) 


are  measures  of  covariance  of  model  and  data  signals  and  of  the  individual  standard  devia¬ 
tions.  Table  3  gives  values  of  pj  for  the  nine  gages  for  both  spectral  and  parabolic  models, 
as  well  as  a  composite  average  p  given  by 


P  = 


(53) 


which  becomes  a  measure  of  accuracy  at  each  angle  of  incidence.  A  plot  of  ^  versus  incident 
angle  is  given  in  Figure  12.  The  plot  indicates  that  the  correlation  for  the  spectral  model 
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is  high  and  does  not  fall  off  with  increasing  an^e  until  the  last  test,  CUSS.  This  would 
be  expected,  since  the  model  does  not  have  any  restriction  on  angle  of  incidence.  The  loss 
of  correlatira  in  CR58  relative  to  the  other  tests  is  due  to  the  strong  apparent  mismatch 
in  phase  on  the  last  three  gages  3-2-1  in  Figure  11,  which  biases  the  overall  average  low 
(see  results  for  this  test  in  Table  3).  In  contrast,  the  correlation  between  parabolic  model 
and  data  is  weaker  and  falls  off  progressively  with  increasing  angle  except  for  an  anomalous 
agreement  on  the  last  test,  CR58.  We  remark  here  that  the  standard  correlation  coefficient 
p,  which  penalizes  mismatches  in  wave  phase  heavily,  can  also  give  high  values  of  correlation 
when  information  is  in  phase,  even  though  there  is  strong  disagreement  between  details  of 
the  compared  time  series.  To  illustrate  this  effect,  the  parabolic  model  time  series  are 
compared  to  data  in  Figure  13  for  CR58.  The  mismatch  in  crest  structures  is  striking,  but 
the  waves  are  phased  relatively  correctly  and  so  the  correlation  coefficient  becomes  high. 

In  order  to  provide  an  independent  check  on  these  results,  an  rms  error  measure  was 
constructed  for  the  comparison.  This  measure  e  is  ^ven  by 


(54) 


for  each  gage,  where  the  normalization  of  the  mesm  square  error  is  with  respect  to  the  true 
(data)  standard  deviation.  A  composite  value  e  for  each  angle  or  test  is  constructed  from 
the  ratio  of  the  rms  error  for  all  points  divided  by  the  total  standard  deviation  for  all  pcunts. 
Table  4  gives  the  computed  error  estimates  for  all  tests,  and  Figure  14  shows  the  composite 
error  as  a  function  of  incidence  angle  for  the  two  models.  In  this  case,  the  trend  towards 
increasing  error  with  increasing  angle  is  clear  in  the  parabolic  model  results.  The  spectral 
model  also  shows  an  increasing  trend  in  this  case,  which  would  not  necessarily  be  expected. 

One  possible  reason  for  an  increase  in  error  in  the  spectral  model  as  wave  angle  becomes 
large  rests  in  the  fact  that  the  reflected  wave  crests  at  the  wall  and  in  the  short-crested 
wave  pattern  away  from  the  wall  result  from  the  interaction  between  waves  that  are  colliding 
head  on  to  the  same  extent  as  they  are  interacting  colinearly.  Tests  and  analysis  of  head  on 
collision  and  vertical  wall  reflection  of  solitary  waves  have  indicated  that  the  leading  order 
theory  (as  in  the  Boussinesq  model  employed  here)  is  not  capable  of  predicting  the  height 
of  the  runup  or  maximum  elevation,  or  the  phase  lag  associated  with  the  opposite-going 
interaction.  See  Su  and  Mirie  (1980)  for  an  example  of  this  type  of  analysis.  Effects  of 
this  nature  may  be  present  in  the  experiments  being  considered  here.  Figure  15  displays 
the  comparison  of  the  time  series  on  the  transect  18-17-10-15-14-13-5  along  the  wall  for 
test  CR580204,  using  spectral  model  results.  Synchronization  is  again  with  respect  to  gage 
13.  It  is  clear  that  the  height  of  the  wave  along  the  wall  is  nnderpredicted  by  the  model, 
as  may  be  expected  for  the  case  <d  low-order  Boussinesq  theory  with  waves  interacting  at 
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large  angles.  This  underestimation  of  the  crest  height  could  in  turn  have  a  detrimental 
effect  on  the  prediction  of  phase  speed;  the  picture  indicates  that  the  model  waves  are 
relatively  slower  than  the  experimental  waves.  The  discrepancies  between  model  and  data 
on  the  transect  along  the  wall  are  even  more  drastic  for  the  parabdic  model,  as  would  be 
expected. 

In  the  comparisons  made  above,  no  attempt  was  made  to  find  a  relative  allignment 
of  model  and  data  results  that  would  minimize  error  in  each  case,  or  to  define  an  error 
measure  that  is  “best”  in  any  sense.  The  results  of  the  comparison  generally  indicate  that 
the  spectral  model  is  better  capable  of  modelling  the  mach  reflection  than  is  the  parabolic 
model. 
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6  Discussion 

We  have  described  the  development  cd  a  solution  technique  for  the  Boussinesq  model  of  long 
waves,  based  on  a  discrete  representation  of  the  angular  spectrum  for  waves  incident  towards 
shore  from  the  ocean.  The  model  is  similar  in  intent  to  the  parabolic  model  developed  earlier 
by  LYK,  with  the  exception  that  the  present  model  does  not  impose  a  restriction  on  the 
range  of  directions  that  can  be  accurately  modelled.  This  advantage  is  counterbalanced 
(at  this  stage  of  development)  by  the  need  to  impose  periodic  boundary  conditions  on  the 
modelled  problem.  This  restriction  must  be  alleviated  before  the  various  representations 
of  angular  spectrum  models  (present  model;  Dalrymple  and  Kirby,  1988;  Dalrymple  et  al, 
1989  )  become  generally  applicable  as  coastal  wave  models. 

A  simple  application  of  the  angular  spectrum  model  in  section  4  showed  that  nonlinear 
interactions  between  modes  travelling  in  different  directions  can  interact  strongly  to  produce 
a  significant  mode  travelling  at  an  intermediate  direction.  In  the  example  here,  it  was 
shown  that  the  amplitude  of  second  harmonic  waves  resulting  from  the  sum  interaction  of 
two  intersecting  primary  waves  may  actually  exceed  the  amplitude  of  the  in-line  second 
harmonics  resulting  from  self-interaction  in  each  primary  wave.  This  result  would  seem  to 
be  counter-intuitive,  since  the  interaction  between  modes  travelling  at  different  directions 
would  seem  to  be  weaker  than  the  interaction  between  modes  with  similar  frequencies 
travelling  in  the  same  direction.  In  the  overall  interplay  of  the  nonlinear  interactions,  this 
simple  picture  does  not  win  out,  either  for  the  case  of  steady  solutions  or  solutions  evolving 
by  means  of  harmonic  generation.  A  recent  paper  (Freilich  et  al,  1990a)  has  shown  evidence 
that  distinct  peaks  resulting  from  higher-harmonic  generation  by  sum  interactions  of  waves 
travelling  in  different  directions  are  observable  in  measured  field  data.  Initial  indications 
are  that  the  parabolic  modelling  scheme  (as  in  LYK)  is  capable  of  modelling  this  interaction 
and  reproducing  the  measured  directional  spectra  after  shoaling  (Fteilich  et  al,  1990b).  A 
comparison  of  the  present  and  the  parabolic  modelling  schemes  and  their  application  to 
field  data  are  underway  and  will  be  reported  separately. 

Finally,  comparison  of  the  present  model,  a  parabolic  model  and  laboratory  data  in¬ 
dicates  that  differences  do  occur  between  predictions  of  a  small-angle  approximation  and 
the  fully  directional  expansion  employed  here,  and  that  the  deviation  is  apparent  even  at 
small  angles  of  incidence.  This  result  indicates  that  the  further  development  of  the  angular 
spectrum  model  is  worth  pursuing  as  a  means  of  providing  accurate  prediction  of  coastal 
wave  fields.  The  first  extension  of  the  present  model  to  include  weak  longshore  topographic 
variation  and  on-offshore  reflection  is  presently  underway  and  will  be  described  shortly. 
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Appendix  A:  Nonlinear  terms 

The  nonlinear  term  [n./.f.]„  appearing  in  (6)  U  ^ven  by 
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Table  1:  Wave  gage  positions  in  Mach  reflection  experiments 
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Table  2:  Paddle  phase  angles  and  directed  wave  angles;  Mach  reflection  experiments 
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Table  3:  Correlation  coefficients  p  for  spectra]  and  parabolic  model  runs 
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Table  4:  RMS  error  coefficients  e  for  spectral  and  parabolic  model  runs 
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1  Introduction 


The  propagation  of  water  waves  over  an  irregular  bathymetry  has  served  as  an  important 
area  of  research  for  mathematicians,  oceanographers  and  ocean  engineers  for  a  long  time 
due  to  the  the  importance  of  the  problem  and  the  mathematical  difficulties  encountered  in 
its  solution.  For  example,  the  design  of  shore  facilities  (harbors,  breakwaters  and  piers) 
and  offshore  structures  requires  a  good  knowledge  of  the  wave  heights  and  directions  to  be 
encountered  at  the  site  of  interest. 

Numerous  models  have  been  proposed  in  the  past  to  solve  various  aspects  of  wave  prop¬ 
agation.  The  first  models  were  ray  tracing  models,  which  determine  the  path  foUowed  by 
the  waves  as  they  traverse  irregular  bathymetry  (e.g.,  Munk  and  Arthur,  1952;  Noda,  1974). 
These  models,  when  used  either  from  offshore  towards  the  shore  or  in  reverse,  provide  rea¬ 
sonable  estimates  for  wave  heights  when  refraction  (due  to  variations  in  water  depth  normal 
to  the  direction  of  propagation)  and  shoaling  (due  to  changes  in  water  depth  in  the  wave 
direction)  are  the  dominant  phenomena  affecting  the  waves. 

Computer  models  of  nearshore  circulation  have  required  knowledge  of  the  wave  field  at 
discrete  grid  points  over  an  offshore  region,  which  led  to  the  development  of  models  based 
on  the  irrotationality  of  the  wave  number.  Examples  of  these  models  are  Perlin  and  Dean 
(1983)  and  Dalrymple  (1988).  Despite  the  inclusion  of  numerous  other  effects  (such  as  wave- 
current  interaction  and  bottom  friction  by  Dalrymple),  none  of  these  models  incorporates 
the  process  of  diffraction.  Ebersole  (1985)  extended  these  models  to  include  diffraction. 

Diffraction,  which  is  the  turning  of  the  wave  rays  due  to  gradients  in  wave  amplitude, 
as  would  occur  as  waves  pass  a  surface-piercing  obstacle  (e.g.,  breakwaters),  was  largely 
ignored  for  coastal  situations  due  to  the  lack  of  a  suitable  model  which  would  incorporate 
refraction  and  diffraction  simultaneously.  In  the  vicinity  of  breakwaters,  the  optics  solution 
of  Sommerfeld  (1898)  was  adapted  by  Penney  and  Price  (1952)  for  waves;  however,  the 
procedure  by  which  refraction  models  would  be  patched  locally  into  the  diffraction  solution 
remained  an  art  for  many  years. 

Of  fundamental  importance  to  the  combined  refraction /diffraction  problem  was  the  devel¬ 
opment  of  the  mild-slope  equation  by  Berkhoff  (1972).  This  equation,  which  is  the  vertically- 
averaged  equation  for  wave  motion,  reduces  the  3-D  problem  of  solving  the  Laplace  equation 
to  a  2-D  problem  for  cases  where  the  bottom  does  not  vary  greatly.  This  second  order 
partial  differential  equation  permits  refraction  and  shoaling,  as  well  as  diffraction,  to  occur 
simultaneously.  Numerous  finite  element  models  were  developed  to  solve  this  equation  (e.g., 
Berkhoff,  1972,  Houston,  1981,  and  Bettess  and  Zienkiewicz,  1977). 

Liu  and  Mei  (1976)  first  developed  a  parabolic  approximation  to  the  problem  of  wave 
fields  in  the  vicinity  of  shore-parallel  breakwaters,  in  order  to  study  the  wave-induced  circu¬ 
lation  caused  by  the  presence  of  the  structure. 

Radder  (1979)  developed  the  first  parabolic  representation  of  the  mild-slope  equation, 
leading  to  greater  computational  simplicity  and  the  ability  to  neglect  the  downwave  boundary 
condition,  which  in  practice  is  very  difficult  to  specify  a  prion.  The  parabdic  modelling  has 
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been  explored  more  folly  by  Kirby  and  Dalrymple  (1983)  amd  Lin  and  Tsay  (1984),  who 
showed  how  to  develop  a  weakly-nonlinear  version  of  the  mild-slope  equation,  allowing  for 
amplitude  dispersion  (bigger  waves  travel  faster  than  smaller  ones). 

One  of  the  drawbacks  of  the  parabolic  modelling  has  been  the  restriction  that  the  waves 
travel  almost  in  a  prescribed  direction.  The  small-angle  parabolic  models  for  example  require 
that  the  waves  travel  within  ±30®  of  the  x  axis.  Wider  angle  models  have  been  developed  by 
Booij  (1981)  and  Kirby  (1986a,b);  however,  each  successive  approximation  only  opens  the 
rsmge  of  allowed  angles  by  a  finite  amount,  and  the  limit  of  full  ±90®  is  never  approached. 
An  alternative  path  taken  by  Lozano  and  Liu  (1980)  was  to  base  the  parabolic  model  on 
the  underlying  refracted  wave  field,  such  that  the  large  angles  due  to  refraction  were  taken 
care  of  properly  prior  to  the  diffraction  calculations.  However,  di&action  effects  still  only 
occurred  in  a  narrow  range  of  angles  around  the  principle  ray  direction. 

Angular  spectrum  modelling  has  in  principle  no  limitation  on  wave  angle;  hence,  its 
development  for  water  waves  by  Dalrymple  and  Kirby  (1988),  Dalrymple,  Suh,  Kirby  and 
Chae  (1989)  and  Suh,  Dalrymple  and  Kirby  (1990).  The  angular  spectrum  approach  has 
been  used  in  other  fields,  such  as  radio  astronomy  (Booker  and  Clemmow,  1950)  and  electro¬ 
magnetic  fields  (Clemmow,  1966).  Some  references  for  the  development  of  angular  spectra 
are  Goodman  (1968;  optics)  and  Stamnes  (1986;  optics  and  water  waves). 

The  basic  ideas  behind  the  angular  spectrum  cam  be  obtained  by  examining  a  single  wave 
train  on  the  surface  of  the  ocean,  which  cam  be  described  at  a  point  (x,y)  in  the  horizontal 
plane  as 

r}{x,y,t)  =  (1.1) 

where  the  real  part  of  rj  is  the  displacement  of  the  water  surface  about  its  mean  position,  a 
is  the  wave  amplitude,  k  and  u  are  the  wave  number  (defined  as  2x  over  the  wave  length 
L)  and  the  amgular  frequency  (2x  over  the  wave  period  T),  respectively.  The  angle  0  is  the 
direction  that  the  wave  direction  makes  with  the  x  axis,  which  points  onshore.  At  z=0,  the 
wave  form  can  be  separated  into  two  parts, 

q(0,y,0  =  (1.2) 

Clearly  there  is  a  oscillatory  variation  of  q  in  the  y  direction,  which  will  be  denoted  by  ^y). 
Defining  A  as  A;  sin  0  we  have 

^y)  =  oe‘^«'  (1.3) 


For  a  more  realistic  sea  state,  the  y  variation  of  ij  will  be  more  complicated,  say,  /(y), 
resulting  from  the  superimposition  of  many  wave  trains  with  different  directions,  but  with 
the  same  frequency  u.  To  determine  the  contributions  of  each  of  the  many  wave  trains  in 
/(y),  we  can  decompose  it  through  the  use  of  a  Fourier  transform  in  the  y  direction.  The 
Fourier  transform  and  its  inverse  for  any  function  /(y)  defined  on  an  infinitely  wide  domain 
are  given  by 


/(A)  =  r 

/(»)  ==  hr  (*•<) 
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where  the  Foaiier  transform  parameter  is  A.  The  angalar  spectrum  is  ^A),  which  consists  of 
the  amplitudes  of  the  wave  trains  travelling  in  directions,  0  =  sin~*(A/ib).  Thus,  the  angular 
spectrum  is  a  (continuous)  collection  of  wave  trains,  each  travelling  in  a  different  directi<», 
determined  by  the  Fourier  parameter,  A.  The  free  surface  displacement  is  now  expressed  as 

>?(*,  ^  (1-5) 

where  k  cos  $  is  rewritten  as  y/k^  —  A^  using  the  definition  of  A.  (We  note  that  the  contri¬ 
butions  from  the  integral  from  |A|  >  k  are  not  propagating  waves,  as  they  decay  in  the  x 
direction.) 

2  The  Angular  Spectrum  Model 

The  sea  surface  for  many  years  has  been  described  by  a  superposition  of  individual  wave  trains 
travelling  in  different  directions,  leading  to  a  directional  spectrum.  The  angular  spectrum  is 
very  similar  except  that  the  directions  are  prescribed  by  Fourier  analysis  for  each  frequency, 
as  we  now  show. 

The  directional  frequency  spectrum  may  be  written 

q(x,y,t)  =  jT  (2.1) 

where  the  real  part  of  q  is  the  water  surface  elevation  as  a  function  of  time  at  position 
{x,y).  F(u;,0)  is  the  amplitude  spectrum  for  the  waves.  If  we  now  assume  that  we  can 
separate  ^(0;,  0)  into  a  separate  frequency  and  direction  components,  F  =  5(u>)Zl(9)  then 
the  above  expression  can  be  rewritten  (replacing  Fsin9  with  the  parameter  A  and  kcos0 
with  y/k^  -  A^)  as 

77(r,  y,  0  =  S(u>)  e'-'du;  (2.2) 

For  the  present  time  we  will  restrict  ourselves  to  a  single  frequency,  with  an  amplitude, 
So  =  S{uo)duf,  and  define  the  angular  spectrum  representation  of  the  waves  at  that  frequency 
as 

qu,(*,  y,  0  =  /_*  (2 J) 

The  angular  spectrum  is  defined  as  the  complex  amplitude  of  the  waves, 

_  2irSoD(A) 

The  wave  trains  making  up  the  angular  spectrum  are  progressive  (that  is,  y/k^  -  A^  is  real); 
however,  for  convenience  and  for  analogy  to  the  Fourier  transform,  we  change  the  limits  cf 
integrati<m  to  ±00. 

»JL,(*,y,0  =  ^  >4(A)e^*’-*’*e’^dA 


(2.4) 


The  new  additional  Mnm  trains  that  are  added  by  extending  the  range  of  integration  are  all 
evanescent,  decaying  in  the  x  direction  as  A  is  greater  than  k  and  the  x  dependency  becomes 
exp{-VA^  —  fc*  i).  In  this  form,  the  wave  spectrum  now  looks  like  the  Fourier  transform 
(1.5). 

For  a  domain  which  repeats  periodically  in  the  transverse  direction,  we  have  the  following 
Fourier  transform  pair, 

/«(*,0  =  (2.5) 

/(*,y,0  =  (2.6) 

— oo 

where  now  \  —  icjh  and  26  is  the  width  of  the  domain.  The  equivalent  angular  spectrum  to 
Eq.  2.4  for  a  given  frequency  is 

»?(*,  y.  0  =  £  (2.7) 


where  u’ave  trains  of  amplitude.  An,  propagate  in  discrete  directions,  measured  by  the  angle 
On  to  the  X  axis. 


On  =  tan 


for  nX  <  k  and  n  =  1, 2, 3, . . . 


(2.8) 


Again,  evanescent  modes  occur  when  nA  >  k.  Figure  1  shows  a  schematic  of  the  angular 
spectrum,  showing  a  number  of  wave  directions.  Sometimes  these  individual  wave  trains  will 
be  referred  to  as  Fourier  modes  or  simply  modes. 

For  domains  bounded  laterally  by  impermeable  barriers,  solutions  are  sought  in  the  form 


OO 


E  /  cos  nXy 

n=0 


(2.9) 


where  again  A  =  x/6,  in  order  that 


=  0,  at  y  =  ±6, 


or,  for  solutions  which  are  non-symmetric  about  y  =  0, 


(2.10) 


i»aO  nsO 

where  the  first  series  represents  the  even  solution  (symmetric  about  y  =  0)  and  the  second 
series  is  the  antisymmetric  part  of  the  solution.  The  parameter  7„  =  (n  + 
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Figure  1;  The  Angular  Spectrum 


2.1  Constant  Depth 


To  fix  ideas,  consider  the  onshore  direction  as  the  x  direction  and  the  transverse  direction  is 
the  y  direction;  the  coordinate  z  will  point  upward  from  the  mean  water  level.  The  velocity 
potential  for  the  wave  motion  is  ♦(x,y,2,t),  from  which  the  velocities  and  water  surface 
elevation  can  be  determined  (see,  e.g.,  Dean  -and  Dalrymple,  1984).  For  linear  plane  waves 
on  constant  depth,  the  boundary  value  problem  is  separable  and  the  reduced  potential  must 
satisfy  the  Helmholz  equation. 


S*0**’-" 


(2.11) 


where  $(z,  y,  z,  t)  =  0(x,  y)  cosh  Jt(/i  +  z)  and  the  wave  number  k  and  the  angular  wave 
frequency  u  are  related  through  the  dispersion  relationship. 


=  gk  tanh  kh 

Fourier  transforming  the  differential  equation  in  the  y  direction  yields 

(jP^ 


dx* 


+  (*’  -  A*)  ^  =  0 


which  is  an  ordinary  differential  equation  for  ^x).  The  solution  is 

^x)  =  Ae*'^*’-^** 


(2.12) 


(2.13) 


(2.14) 


where  A  is  a  function  of  A  (but  constant  in  x)  which  must  be  determined.  The  Fourier 
inverse  in  an  infinitely  wide  domain  is 


(2.15) 
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At  X  =  0, 


^o,y)=  r  A(xy^ydx  (2.16) 

•/— oo 

Comparing  with  Equations  (1.4),  we  see  the  boundary  condition  for  A(A), 

A(A)=  r  4>(0,y)e-'^'>dy  (2.17) 

•/—OP 

It  should  be  emphasized  that  the  solution  we  have  obtained  (2.15,  2.17)  for  this  linear  wave 
field  is  y  the  initial  condition  at  x  =  0.  as  there  is  no  coupling  between  the  various 

modes  com^.  ising  ihe  angular  spectrum.  If  the  problem  of  interest  is  wave  propagation 
behind  a  breakwater,  with  all  the  attendent  diffraction,  the  angular  spectrum  model  shows 
that  the  wave  behavior  is  solely  determined  by  the  initial  condition,  t.e.,  the  diffraction  the 
waves  experience  behind  a  breakwater  is  determined  by  the  initial  condition  at  the  breakwater 
as  each  of  the  wave  trains  comprising  the  angular  spectrum  propagate  independently  of  each 
other.  The  diffraction  pattern  observed  is  created  by  the  radial  spreading  of  the  wave  trains 
comprising  the  angular  spectrum. 

For  a  periodic  domain  of  lateral  extent  26,  the  Fourier  transform  leads  to 

^  +  (fc’-(«A)^)^=0  (2.18) 

»*ith  the  solution  _ 

^x,nA)  =  (2.19) 

The  inverse  is 

♦(x,y,z,t)  =  v-<^)  cosh  k{h  +  z)  (2.20) 

-OO 

The  wave  number  vector  k  =  {kcos6,ksin0)  =  (y/k^  -  (nA)^,nA). 

Again,  the  behavior  of  the  wave  trains  as  they  propagate  is  fixed  by  the  initial  conditions 
in  this  problem,  which  are  ascertained  by  a  Fourier  transform  of  the  initial  conditions, 
^n(0,y)  =  A„.  For  example,  for  a  plane  wave  train  propagating  at  angle,  7  to  the  x  axis, 
the  initial  value  of  ^0,y)  =  where  Xo  =  k  an  t.  The  Fourier  transform  of  this  initial 

condition,  again  for  a  periodic  domain,  is 

An  =  2Z16sin(Ao  -  nA)6, 

unless  Aq  =  mA  for  some  integer,  m,  in  which  case  all  the  An  are  zero,  except  for  Am,  which 
is  2Db.  For  computational  purposes,  it  is  convenient  to  chose  incident  plane  wave  trains 
which  correspond  to  a  Fourier  mode. 

The  conservation  of  energy  equation  can  be  used  as  a  check  for  numerical  calculations  in 
periodic  domains,  as  the  energy  fluxes  through  the  lateral  boundaries  cancel  and  therefore 


the  flux  of  energy  across  any  plane  parallel  to  the  y  axis  is  a  constant.  If  the  energy  flux  is 
defined  as 

then  the  flux  past  any  location  in  terms  of  the  Fourier  coefficients,  An,  can  be  shown  to  be 

CC  r 

pu - ^6y^|.4n|^yA;^  — (nAP  =  constant  (2.22) 

^  n 


2.1.1  Wavemakers 

The  solution  for  <6(x,y)  in  constant  depth  water  (2.15,2.17)  can  be  used  to  determine  the  wave 
field  inside  directional  wave  basins,  which  are  typically  rectangular  basins  with  segmented 
wavemakers  along  one  wall  (Dalrymple  and  Greenberg,  1985).  For  an  example,  we  will 
assume  the  basin  is  infinitely  wide  (y  direction)  with  the  x  axis  located  at  the  center  of  the 
wavemaker  pointing  into  the  basin. 

Equation  (2.16)  indicates  that  the  wave  field  is  known,  once  the  velocity  potential  is 
specified  along  the  wavemaker;  however,  the  usual  linearized  wavemaker  boundary  condition 
specifies  the  velocity  in  the  x  direction  at  z=0  (see  Dean  and  Dalrymple,  1984,  for  example). 
Therefore  we  need  to  treat  this  problem  slightly  differently. 

We  will  take  the  horizontal  velocity  created  by  the  wavemaker  of  length  2a  (associated 
with  the  progressive  wave  mode)  to  be  adequately  described  by 


<^<Ko.y)  f  |y|  <  o 

dx  \  0,  |y| >  o 

(2.23) 

where  Aq  =  ifesin7  and  7  measures  the  desired  wave  direction.  The  Fourier  transform  of  this 

condition  is  _ 

d^0,A)  2sin(Ao-A)a 

dx  ~  (Ao-A)  ^ 

From  the  solution  for  <t>{x,y),  we  have 

=  T-  -  A2  A(A)  dA 

of  2x  V— 00 

(2.25) 

and,  after  transforming, 

=  tV*»  -  AM(A) 

(2.26) 

Equating  these  two  expressions  for  the  transformed  velocity,  we  find 

,,,,  2i6in(Ao-A)o 

^  (Ao  -  A)v'*’ -  A3 

(2.27) 
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which  is  also  obtainable  through  a  Greens  function  approach.  Dalrymple  and  Greenberg 
(1985)  also  treated  the  evanescent  modes  which  are  important  in  the  vicinity  of  the  wave- 
maker.  Stamnes  (1986)  shows  how  to  obtain  the  evanescent  modes  in  the  context  of  an 
angular  spectrum. 

Dalrymple  and  Kirby  (1988)  point  out  that  the  wave  held  behind  a  thin  island  of  width 
2a  can  be  easily  found  by  taking  the  plane  wave  solution  and  subtracting  the  wavemaker 
solution  from  it,  according  to  Babinet’s  Principle. 

Dalrymple  (1989)  has  used  this  constant  depth  solution  and  the  non-symmetric  Fourier 
expansion  (2.10)  to  determine  the  wave  fields  generated  by  directional  wavemakers  in  basins 
of  finite  width,  taking  into  account  side  wall  reflection.  This  is  the  so-called  ‘designer  waves' 
solution.  (His  solution  is  also  valid  for  basins  with  straight  and  parallel  bottom  contours.) 

To  illustrate  the  ‘designer  wave’  concept,  see  Figure  2,  which  show  the  instantaneous 
water  surface  for  a  wave  train  with  30**  angle  of  incidence  which  is  generated  at  the  wavemaker, 
located  at  the  bottom  of  the  figure.  Due  to  the  presence  of  the  reflecting  wall  at  the  left  of 
the  figure,  a  short-crested  wave  pattern  is  created.  The  presence  of  the  side  wall  at  the  right 
of  the  figure  creates  a  diffraction  zone.  There  is  only  a  limited  portion  of  the  wave  basin 
near  the  wavemakers,  where  the  deared  wave  train  exists.  At  the  far  end  of  the  tank  (at 
the  top  of  the  figure),  the  wave  field  is  clearly  no  longer  a  uniform  wave  train.  Now,  if  this 
long-crested  wave  train  is  desired  to  occur  at  the  far  end  of  the  tank,  then  the  power  of  a 
directional  wavemaker  can  be  used.  By  generating  a  nonuniform  wave  field,  the  presence  of 
the  side  walls  can  be  incorporated  in  the  wavemaker  signals.  For  example,  to  generate  a  wave 
train  with  30°angle  of  incidence  at  the  far  end  of  the  basin,  the  wavemaker  must  generated 
short-crested  and  diffracted  waves.  For  this  case,  if  we  now  imagine  the  wavemaker  at  the 
top  of  the  same  figure,  we  see  that  by  the  time  the  wave  field  propagates  to  the  far  end  of 
the  tank  (now  the  bottom  of  Figure  2),  the  wave  becomes  long-crested.  This  technique, 
including  shoaling  waves,  is  used  in  the  Ocean  Engineering  Laboratory  at  the  University  of 
Delaware. 


2.1.2  Waves  Behind  Narrow  Gaps 

Another  interesting  application  of  this  constant  depth  solution  is  the  propagation  of  a  direc¬ 
tional  sea  state  through  a  small  gaq>.  This  problem  can  be  treated  in  several  ways.  First  the 
influence  of  the  waves  on  the  seaward  side  of  the  breakwater  could  be  ne^ected,  in  which 
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Figure  2;  Instantaneous  Wave  Field  Generated  by  a  Directional  Wavemaker  in  a  Basin  with 
Reflecting  Sidewalls,  from  Dalrymple  (1989) 


case  the  wavemaker  solution  could  be  used  directly  with  the  appropriate  specification  of  the 
velocity.  Alternatively  the  offshore  wave  field  can  be  specified  and  matched  at  the  gap  to 
the  wavefield  on  the  sheltered  side  of  the  breakwater.  This  is  the  case  we  will  treat. 

Consider  either  a  natural  gap  through  which  waves  must  propagate  into  an  embayment  or 
a  man-made  gap,  as  between  two  breakwater  segments.  It  has  been  noticed  that,  regardless 
of  how  directional  the  sea  state  seaward  of  the  gap  is,  the  waves  inshore  of  the  gap  are 
long-crested  and  these  crests  are  circular.  To  show  that  this  is  true,  we  postulate  an  angular 
spectrum  for  the  waves  offshore  of  the  (infinitessimally  narrow)  gap  for  a  given  frequency  as 

U^,y)  =  ^  (2.30) 

where  the  A{X)  are  imposed  by  the  nature  of  the  offshore  sea  state.  The  velocity  at  the  gap 

^  A(A)e'^  dX  Siy)  (2.31) 

where  the  S(y)  is  the  Dirac  delta  ftmction,  which  represents  the  fact  that  the  velocity  is  zero 
except  at  the  narrow  gap.  The  Fourier  transform: 


dMo^y)  _  J_ 

dx  2it 


(2.32) 


Inside  the  (infinitely  thin)  breakwater,  the  wave  field  for  this  frequency  is  given  by  another 
angular  spectrum, 

4>i{x,y)  =  r  B{X)  dX  (2.33) 

J-^OO 

The  transformed  horizontal  velocity  (z  direction)  at  the  gap  is 


(2.34) 

(2.35) 


We  can  make  several  assumptions  here  for  A(A).  If  we  assume  that  the  wave  trains  are 
coming  from  all  directions  with  the  same  height  (a),  then  A(A)  =  o;  alternatively  if  the 
waves  come  only  from  the  direction  associated  with  Ao,  then  A{X)  =  A(Xo)S{Xo  -  A).  For 
both  cases,  the  integral  in  the  above  expression  can  be  evaluated,  such  that 


B{X)  = 


ak^ 

Wk^  -  X^ 


(2.36) 


or 


B{X)  = 


-iAjXo) 

-  A5 


(2.37) 
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Substituting  these  expression  into  the  angular  spectrum  inside  the  breakwater,  (2.33),  and 
integrating,  we  have 

^i(*.y)  =  QfHo(kr)  (2.38) 

where,  for  the  first  case,  a  =  ak^/2  and,  for  the  second  case,  a  =  -*>l(Ao)/2.  The  Hankel 
function  of  the  first  kind  and  order  zero,  which  is  a  function  of  the  wavenumber  and  the 
radial  distance  r  =  y/x^  +  y^,  arrives  through  the  following  identity: 


1  d\ 

2x  y_oo  y/k^  -  A2 


=  I  Ho(kr) 


(2.39) 


The  nature  of  our  solution  for  one  frequency  is  that  the  wave  form  inside  the  bay,  which 
is  only  dependent  on  the  radial  distance  from  the  gap,  is  described  by  circular  wave  crests 
(via  the  Hankel  function).  For  the  linear  superposition  of  many  frequencies,  the  same  result 
obtains-all  frequencies  are  described  by  circular  wave  crests. 

For  waves  through  a  breakwater  gap  much  smaller  than  a  wavelength,  the  result  is  not 
new,  Penney  and  Price  (1952). 

Dalrymple  and  Martin  (1990)  examine  the  wave  field  inshore  of  a  line  of  breakwaters  that 
are  separated  by  gaps  of  the  same  length  (similar  to  an  optical  grating).  The  influence  of 
the  offshore  waves  is  included  in  their  analysis.  They  find  that  the  wave  field  inshore  of  the 
breakwater  can  be  very  complicated  as  new  wave  modes  are  generated  when  the  gap  spacing 
is  less  than  a  wave  length,  due  to  superposition  of  the  dif&action  patterns  behind  each  gsq>. 
The  presence  of  multiple  wave  trains  of  the  same  frequency  can  lead  to  the  formation  of  rip 
currents  behind  such  structures. 


2.1.3  Wave  Focussing 

Stamnes  et  al.  (1983)  carried  out  a  laboratory  and  theoretical  study  of  the  focussing  of 
waves  behind  a  submerged  shoal,  designed  to  act  as  a  Fresnel  lens.  The  field  experiment 
consisted  of  generating  a  circular  wave  with  a  point  wavemaker;  these  waves  were  focussed 
by  the  lens  to  another  (focal)  point.  Wave  height  measurements  were  made  with  densely 
spaced  transects  in  the  vicinity  of  the  focal  point.  The  comparisons  to  the  data  were  made 
using  the  (linear)  angular  spectrum  approach,  along  with  a  nonlinear  parabolic  model.  The 
comparison  was  carried  out  by  assuming  the  wave  field  over  the  lens  was  representable  by  a 
arc  of  a  circular  wave  converging  on  the  focal  point. 


3  Variable  Bathymetry 

Most  coastal  refpons  of  interest  are  not  characterized  by  uniform  depth,  but  instead  have 
spatial  nonuniformity  in  bathymetry.  At  the  simplest,  this  spatial  variation  is  characterized 
by  a  trend  of  decreasing  depth  in  the  shoreward  direction,  with  a  superposed,  irregular 
depth  variation  in  both  the  on-offshore  and  longshore  direction.  For  the  case  oi  intermediate 

I 


I 


12 


water  depth,  it  is  convenient  to  model  small  amplitude  waves  using  the  mild-slope  equation 
(Berkhoff,  1972;  Smith  and  Sprinks,  1975), 

V  •  CCgV<l>  +  k^CCg(t>  =  0  (3.1) 

which  is  applied  separately  to  each  frequency  component  in  the  wave  train.  Here,  V  = 
(^t  +  ^j)  is  the  horizontal  gradient  operator,  C  is  the  wave  phase  celerity  and  Cg  is  the 
group  velocity  of  the  waves,  given  by 

C{x,y)  =  u,/k{x,y)  (3.2) 

Again,  the  complete  velocity  potential  is  #(i,p,z,t)  =  ^i,y)coshfc(/»  +  z)c~“^.  The  coeffi¬ 
cients,  C  and  Cf,  in  (3.1)  are  determined  based  on  the  local  value  of  the  water  depth  h(x,y) 
at  each  point  in  the  domain  of  interest  and  the  dispersion  relationship  (2.12). 

3.1  Straight  and  Parallel  Contours 

The  Fourier  transform  of  (3.1)  for  an  infinitely  wide  domain  where  the  depth  only  varies  in 
the  X  direction  is 

^ 

This  ordinary  differential  equation  for  ^  has  variable  coefficients  as  the  wave  number,  wave 
phase,  and  group  velocities  vary  with  the  depth,  h{x). 

An  assumed  form  of  the  solution  for  ^  is 

^i)  =  <^(x)e'/'/^^  (3.5) 


where  d  is  assumed  to  vary  slowly  in  x  as  the  exponential  term  carries  most  of  the  phase  in¬ 
formation  and  the  integral  is  necessary  to  get  the  phase  change  with  x  correctly.  Substituting 
into  equation  (3.4)  leaves  an  equation  for 


> 

'  dx  ax 

(3.6) 

where  we  have  neglected  two  small  terms. 

CC^ 

<  1 

(3.7) 

dCCgd^ 
dz  dx 

<  1 

(3.8) 

as  the  bottom  varies  slowly  in  x  as  does 


I 


The  resulting  equation  is  a  first  order  differential  equation  for  the  wave  modes  which  can 
be  solved  by  separation,  resulting  in 


,  .  ^CC.^Jk^  -  x\ 

4>  =  Ml-  — 

sjCCgy/k'^  -  A2 


AoA'.(A)Jirr(A) 


(3.9) 


or 

^  (3.10) 

where  the  terms  with  a  zero  subscript  axe  evaluated  at  the  origin.  This  equation  shows 
that  the  (transformed)  wave  form  is  ^ven  by  an  initial  amplitude,  A,  evaluated  at  x  =;  0 
and  then  the  wave  changes  with  x  in  accordance  with  standard  shoaling  and  refraction 
coefficients,  K»Kr-  Furthermore  there  is  no  interaction  between  the  different  Fourier  modes; 
any  diffraction  observed  for  x  >  0  is  due  to  the  spreading  of  the  angular  spectrum,  as  each 
wave  train  propagates  in  a  different  direction. 


The  inverse  Fourier  transform  yields 

§  =  r  (3.11) 

ZK  Joo 

for  an  infinitely  wide  domain.  This  is  a  generalization  of  the  result  of  Mei  et  ai  (1968),  who 
used  a  multiple  scales  approach,  to  now  include  all  possible  wave  trains  in  a  directional  sea. 


Dalrymple  and  Kirby  (1988)  studied  the  case  of  waves  propagating  through  gaps  in  a 
row  of  periodically  spaced  offshore  breakwaters,  assuming  that  the  potential  through  the 
gaps  could  be  ^ven  as  the  wave  potential  in  the  absence  of  the  breakwaters  (known  as  the 
Kirchoff  approximation  in  optics).  Shoreward  of  the  breakwater  the  assumed  planar  bottom 
sloped  upwards  to  the  shorefine.  Their  results,  see  Figure  3,  show  diffraction  patterns  behind 
each  of  the  breakwater  gaps  as  would  be  expected,  including  the  refraction  of  the  waves  due 
to  the  sloping  bottom.  There  is  also  a  region  of  short-crested  waves  due  to  the  influence  of 
neighboring  gaps. 


3.2  Realistic  Bathymetry 

For  more  realistic  bathymetry,  we  must  permit  a  variation  in  the  bathymetry  in  the  y  di¬ 
rection  as  well.  It  is  convenient  to  modify  the  mOd-slope  equation  to  treat  this  problem. 
Introducing  p(x,y)  =  CCg  and  ^  we  obtain  a  Helmholz  equation  for  the  modified 

potential, 

(3.12) 

where 

=  (3.13) 

Fdlowing  Dalrymple  et  a/.,  1989,  a  laterally  averted  wave  number  is  introduced  for  a  domain 
of  width,2fi. 
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This  pennits  rewriting  the  Helmholtz  equation  as 

=  0  (3.14) 

where  y)  represents  the  lateral  deviation  of  it(*)  from  Are(*,  y): 

Fi/2  =P-jk2. 

For  bathymetries  which  have  very  little  variation  in  the  y  direction,  the  is  small. 

Fourier  transforming  (3.14)  yields  an  equation  for  the  Fourier  modes, 

0  +  (F  -  A^)^  -  Ff’(i/2o)  =  0  (3.15) 

where  the  Fourier  transform  of  the  product,  is  shown  symbdically  with  the  F  operator. 


This  second  order  differential  equation  can  be  separated  into  two  equations  (as  is  done  for 
developing  parabolic  equations),  one  equation  governing  waves  propagating  in  the  positive 
X  direction  and  another  for  the  opposite  direction,  which,  after  assuming  that  the  nega¬ 
tively  propagating  wave  motion  is  small,  due  to  only  small  amounts  of  reflection,  yields  one 


equation. 


dx 


*2-  A2 


2\/F- A2 


t 

2\/W^ 


(3.16) 


The  second  term  on  the  right  hand  side  is  the  shoaling  and  reixaction  term,  which  occurs  for 
planar  bathymetry,  and  the  last  term  represents  the  effect  of  the  irregular  bathymetry.  This 
bathymetric  effect  represents  a  coupling  between  all  Fourier  wave  modes  and  the  bottom  due 
to  refraction.  Modes  which  do  not  exist  can  be  created  by  this  interaction. 


The  example  used  by  Dalrymple,  Sub,  Kirby  and  Cbae  (1989)  is  that  of  a  submerged 
circular  shoal  located  in  a  region  of  constant  depth,  as  studied  in  the  laboratory  by  Ito 
and  Taiumoto  (1972).  The  evolution  of  the  modes  over  and  behind  the  shoal  are  shown  in 
Figure  4.  The  incident  wave  train  is  specified  to  have  only  one  direction  (mode).  As  the  wave 
train  encounters  the  shoal,  the  waves  refract  and  focus  due  to  the  depth  changes.  Through  the 
bottom  term,  new  Fourier  modes  are  spawned  and  grow  over  the  island.  Behind  the  island, 
in  the  constant  depth  region,  the  bottom  coupling  term  is  no  longer  active  and  the  angular 
spectrum  is  unchanged  with  further  propagation  distance,  yet  the  wave  field  experiences  the 
formation  of  a  strong  focal  re^on,  where  diffraction  is  important.  The  apparent  diffraction 
once  again  is  explained  by  the  radial  propagation  of  the  individual  members  of  the  angular 
spectrum,  generated  by  refraction  over  the  shoal.  The  envelope  of  the  wave  field  in  the 
vicinity  of  the  shoal  is  shown  in  Figure  5. 


3.3  Alternative  Formulation  of  the  Bottom  Coupling  Term 

Due  to  the  efficiency  of  the  FFT  algorithm,  the  Fourier  transform  of  the  last  term  in  (3.16) 
and  the  inverse  transform  needed  to  obtain  ^  from  ^  are  usually  evaluated  over  a  periodic 
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Figure  4:  Evolution  of  the  Angular  Spectrum  over  a  Circular  Shoal,  from  Dalrymple  et  al., 
1989.  (a)  Seaward  of  Shoal,  (b)  Over  the  Top  of  the  Shoal,  (c)  Just  Landward  of  the  Shoal, 
(d)  Far  Landward  of  the  Shoal  (note:  no  change  from  (c)  to  (d)). 


Figure  5:  Envelope  of  the  Wave  Amplitude  Behind  a  Circular  Shoal  for  45*of  Incidence,  fran 
Dalrymple  et  al.,  1989 


interval  in  a  finite  width  domain.  Using  this  fact  abont  the  computational  procedure  allows 
the  effect  of  the  coupling  term  to  be  directly  interpreted  in  terms  of  mode  coupling  between 
transverse  surface  wave  modes  and  transverse  bottom  modes.  In  order  to  show  this,  let 


f{y)  =  (3.17) 

denote  the  real-valued  argument  of  F  on  the  line  —b<y<b.  The  Fourier  transform  pair 
for  /  is  then  defined  by  the  transform  pair,  (2.5,  2.6). 


/(»)=  E ^  =  T 

f|=  — oo 

(3.18) 

(3.19) 

Likewise,  the  discrete  transforms  4>n^  physical  variables 

are  defined  by 

(3.20) 

(3.21) 

The  discrete  Fourier  transform  of  the  governing  equation  (3.15)  gives 

^  +  (F  _  («A)2)<^  -  F/n  =  0  (3.22) 

By  direct  substitution,  we  obtain 

/(y)  =  (3-23) 

I  "» 

Then 

I  in 

=  EE‘'i'^"‘«(^ (3-24) 

/  m 

where  d(s)  =  1  for  s  =  0  and  zero  otherwise.  We  thus  require  the  condition  /  -I-  m  —  n  =  0 
in  order  to  obtain  a  contribution  to  the  surface  wave  mode  by  the  bottom.  For  an  arbitrary 
choice  of  I  and  a  subsequently  fixed  value  of  m  =  n  - 1,  we  obtain  the  expressicm 

(3.25) 

(s— oo 
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Figure  6:  Amplitude  Spectrum  of  the  Bottom  Modes,  i/^(z,mA),  from  Dalrymple  and  Suh, 
1988 

When  this  term  is  substituted  back  into  (3.22),  it  is  clear  that  the  effect  of  Fn  is  to  scatter 
energy  into  (or  out  of)  mode  n  through  resoawt  interaction  between  a  surface  wave  at  mode 
n  -  I  and  a  bottom  mode  L  The  interaction  is  a  linear  process,  since  the  amplitudes  vf 
are  fixed.  The  formulation  may  be  interpreted  as  a  generalization  of  the  Bragg  scattering 
mechanism  described  by  Mei  (1985)  for  forward  propagation  over  a  sinusoidal  bed  (actually  - 
not  explicitly  covered  here)  and  by  Naciri  and  Mei  (1988)  for  waves  propagating  over  a  bi- 
periodic  bottom.  In  this  extension,  the  entire  problem  of  wave  deformation  by  a  non-uniform 
bottom  may  be  viewed  as  the  result  of  a  complicated  multiple  scattering  problem  involving 
the  entire  set  of  Fourier  modes  resolved  by  the  system. 

In  practice,  it  is  possible  either  to  use  the  FFT  and  inverse  FFT  algorithms  to  evaluate 
the  last  term  in  (3.22),  or  to  evaluate  the  sum  in  (3.25)  after  evaluating  the  FFT  of  The 
operation  count  of  the  second  option  is  smaller,  although  all  exponentials  involved  in  the 
computation  of  the  FFT’s  may  be  evaluated  once  and  then  stored,  making  the  subsequent 
FFT  calls  quite  efficient. 

Dalrymple  and  Suh  (1988)  show  this  interaction  between  the  bottom  and  surface  modes, 
using  an  idealized  bathymetry  consisting  of  sinusoidal  corrugations  perpendicular  to  the  t 
axis,  so  that  the  major  component  of  the  bottom  was  vf.  The  surface  wave  train  was  incidoit 
at  an  angle  corresponding  to  ^4  at  zsO.  According  to  (3.25),  the  next  modes  to  be  forced 
by  the  bottom  are  ^  and  de-  These  modes  then  create  additional  modes.  These  results  are 
shown  on  Figures  6-8,  which  show  the  evolution  of  the  angular  spectrum  with  distance,  «, 
and  the  instantaneous  water  surface.  Note  that  wave  rays,  obtuned  from  refraction  theory 
and  depicted  with  sc^d  lines,  are  also  shown. 
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Figure  7:  Evolution  of  the  Angular  Spectrum  at  Different  x  Locations,  from  Dalrymple  and 
Suh,  1988 

0 

3.4  Nonlinear  Effects 

The  solutions  discussed  so  far  have  been  for  Hnear  theory,  which  permits  the  summation 
over  many  wave  directions.  The  water  wave  problem  however  is  nonlinear,  due  to  the  free 
surface  boundary  conditions.  There  have  been  two  approaches  to  solving  a  more  nonlinear 
problem.  The  first  is  to  replace  the  linear  dispersion  relationship  (2.12)  used  to  find  the  wave 
number  given  the  wave  frequency  with  a  nonlinear  dispersion  relationship,  which  includes 
amplitude  dispersion.  Candidates  for  this  disperaon  relationship  could  be  the  Stokes  third 
order  relationship  or  the  modified  dispersion  relationship  of  Kirby  and  Dalrymple  (1986). 
This  last  relationship  was  used  directly  in  the  equations  in  Section  3.2,  with  very  good 
agreement  with  laboratory  data  by  Dalrymple  et  al.  (1989). 

The  second  approach  is  to  rigorously  satisfy  the  nonlinear  governing  equations,  which 
are 


vH  = 

0  {-h  <  z  <  ij) 

(3.26) 

0  (z  =  1?) 

(3.27) 

♦l  +  “|V*|*  +  STJ  = 

a 

N 

II 

(3.28) 

♦.  = 

(z  =  -h) 

(3.29) 

where  V  and  V4  are  the  three>dimensional  and  the  horizontal  gradient  operatcm. 
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Figure  8:  Instantaneous  Water  Surface  Elevations  Over  Periodic  Ridges,  from  Dalrymple 
and  Sub,  1988 
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The  methodology  is  to  expand  $  and  17  in  terms  of  a  small  parameter,  e,  which  is  the 
Stokes  steepness  parameter, 


00  00 

^  =  11  »?  =  (3.30) 

n=l  n=l 

Also  the  method  of  multiple  scales  is  used,  which  has  been  shown  to  be  very  useful  for  wave 
evolution  problems.  Suh  et  al.  (1990)  chose  the  following  slow  variables 

xi  =  =  et,i2  =  (3.31) 

A  mild  slope  bottom  was  assumed,  such  that  hx  ~  .  This  leaves  the  bottom  effectively 

characterized  by  straight  and  parallel  contours  up  until  third  order.  Next,  the  Taylor  series 
is  used  to  expand  the  nonlinear  free  surface  boundary  conditions  about  z  =  0  and  the 
bottom  boundary  condition  about  z  =  —  A,  which  provides  for  series  representations  of  these 
nonlinear  conditions  in  linear  form,  ordered  by  the  Stokes  parameter.  Finally,  grouping  the 
coefficients  of  each  order  of  the  Stokes  parameter  gives  a  boundary  value  problem  for  each 
order,  v: 


<1 

II 

{-h<Z<  0) 

(3.32) 

+  ^n.,  =  Gn 

(z  =  0) 

(3.33) 

+  9  =  fffi 

(z  =  0) 

(3.34) 

<l>n,  =  Bn 

(z  =  -h) 

(3.35) 

where  the  forcing  terms,  F„,Gn,ffnil3n,  are  determined  from  solutions  of  lower  order.  See 
Suh  et  ai  (1990)  for  the  complete  expressions  for  the  right  hand  sides.  The  first  order 
solution,  as  expected,  is  that  of  Dalrymple  and  Kirby  (1988),  for  waves  on  straight  and 
parallel  contours.  At  the  second  order,  the  interaction  of  components  of  the  first  order 
directional  spectrum  leads  to  the  appearance  of  sum  and  difference  frequencies,  including  the 
usual  Stokes  harmonic.  Additionally,  forcing  from  the  bottom  occurs  through  the  coupling 
term  in  Dalrymple  et  ai  (1989).  At  the  third  order,  the  forcing  of  the  wave  field  is  due  to 
third-order  terms  proportional  to  the  first  harmonics  and  a  cubic  resonance  which  results 
from  the  interaction  of  the  primary  waves  and  the  second-order  sum  and  difference  waves  or 
among  three  primary  waves.  New  bottom  coupling  terms  arise  as  well. 

Suh  et  al.  (1990)  compared  the  third-order  model  to  data  (Bomj,  Berkhoff  and  Radder, 
1982,  Ito  and  Tanimoto,  1972),  showing  good  agreement,  far  better  than  that  obtained 
by  linear  theory.  Another  comparison  was  made  to  the  nonlinear  model  of  Dalrymple  et  al. 
(1989)  on  the  BBR  data  set,  also  showing  good  s^eement  and  indicating  that,  at  least  for  the 
example  shown,  the  two  methods  of  including  nonlinearity  are  almost  equivalent.  In  fact, 
the  Dalrymple  et  al.  (1989)  model  has  a  computational  advantage  for  most  applications; 
however  it  does  not  include  the  wave-wave  interactions  as  does  Suh  et  al.  (1990).  Both 
models  have  the  disadvantage  that  for  very  large  angles,  say  greater  than  50*,  the  model 
results  be^n  to  show  discrepancies. 
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4  Shallow  Water  Waves 


As  waves  propagate  towards  shore,  they  enter  a  re^on  where  the  wavelength  becomes  long 
relative  to  the  water  depth,  and  the  product  kh  in  the  dispersion  relation  (2.12)  becomes 
small.  In  this  case,  the  dispersion  relation  approaches  the  limiting  form 

=  gk^h  (4.1) 

and  wave  speed  C  =  u/k  becomes  only  weakly  dependent  on  frequency.  In  this  limit,  all 
waves  are  travelling  at  nearly  the  same  speed,  and  strong  nonlinear  coupling  exists  between 
waves  of  different  frequency  and  direction.  Prior  to  wave  breaking,  numerous  studies  have  in¬ 
dicated  that  the  Boussinesq  equations  (in  which  all  frequencies  are  treated  together)  provide 
an  adequate  model  for  the  wave  field.  To  date,  this  problem  has  only  been  treated  (from  the 
angular  spectrum  point  of  view)  for  the  case  of  topography  varying  in  one  direction;  the  case 
of  two-dimensional  topography  has  not  been  described  as  of  yet.  Kirby  (1990)  has  developed 
a  model  for  waves  on  an  infinitely  long  beach  and  has  compared  model  results  to  data  and 
parabolic  model  results.  This  model  will  be  presented  here. 

We  first  establish  the  form  of  a  model  for  waves  in  a  laterally  unrestricted  domain.  Again, 
a  Cartesian  coordinate  system  is  adopted  which  has  x  pointed  in  the  onshore  direction  and 
y  pointing  alongshore.  Depth  is  assumed  to  vary  as  h{x)  only.  We  take  as  a  starting  point 
the  variable  depth  Boussinesq  equations  as  ^ven  by  Peregrine  (1967): 


fh  +  V  •  (hu)  -f  cV  •  (qu)  =  0  (4.2) 

u,  +  oi  •  Vu  +  pVq  =  /1^{-V(V  •  (hut))  -  y  V(V  •  u*)}  (4.3) 

Here,  q  is  the  surface  displacement  and  u  is  the  horizontal  wave-induced  velocity  vector.  The 
equations  are  kept  in  dimensional  form;  the  scaling  parameters  e  for  nonlinearity  and  for 
weak  dispersion  are  present  only  schematically  and  will  be  subsequently  dropped.  We  will 
assume  that  either  bottom  slope  or  the  amplitude  of  bottom  features  (as  scaled  by  water 
depth)  are  also  small  and  hence  the  model  will  be  developed  to  leading  order  in  nonlinearity, 
dispersion,  and  domain  inhomogeneity.  This  leads  to  immediate  neglect  of  bottom  slope 
effects  in  the  dispersive  terms  of  (4.3).  Using  the  linear  portion  of  (4.2),  we  may  then  write 
(4.3)  in  the  reduced  form 


Ut  -I-  u  •  Vu  -I-  pVq  -I-  -  Vq,i  =  0. 

o 


(4.4) 


We  now  make  the  following  two  assumptions.  First,  the  model  will  be  applied  to  time- 
periodic  wave  trains,  where  periodicity  is  in  the  sense  of  either  a  regular  wave  train,  or  of  a 
discrete  FFT  over  a  finite  length  of  sampled  data.  This  assumption  has  been  routinely  ^plied 
in  the  spectral  sense  in  the  one-directional  shoaling  model  of  Freilich  and  Gnza  (1984)  (see 
also  Elgar  and  Guza,  1985).  Its  interpretation  in  the  regular  wave  case  is  strai^tforw^, 
with  the  wave  being  separated  into  its  harmonic  components  (Rogers  and  Mei,  1978;  Lin, 
Yoon  and  Kirby,  1985,  hereafter  referred  to  as  LYK).  Secondly,  the  wave  field  will  be  assumed 


to  be  periodic  in  the  transverse  y  direction.  This  corresponds  again  to  a  fixed  longshore 
wavelength  in  the  regolar  wave  case,  or  to  periodicity  over  a  long  spatial  interval  in  the 
spectral  sense. 

The  governing  equations  are  first  split  into  coupled  elliptic  models  for  separate  harmonic 
components.  Following  LYK,  the  surface  displacement  and  velocity  are  written  as 

_  \  V) inwl  .  -  .  /  A  r\ 

17  =  2^  — 5 — e  +  C.c.  (4.5) 

n=0  ^ 

+  (4.6) 

n=0  ^ 

Substitution  of  (4.5)  and  (4.6)  in  (4.2)  and  (4.4)  and  subsequent  elimination  of  the  velocity 
leads  to  the  foUowing  model  equation  for  the  in  the  horizontal  plane: 

+  V  •  (Gn^ih)  +  [n./.t.)„  =  0;  n  =  1, N  (4.7) 

Here,  [n./.t]„  denotes  the  nonlinear  interactions  with  other  discrete  frequency  components 
which  are  sorted  by  means  of  the  rules  for  triad  interactions  applied  to  the  time  dependence; 
complete  expressions  may  be  found  in  Kirby  (1990).  The  mode  n  =  0  corresponding  to  the 
steady,  wave-induced  setdown  is  neglected  since  it  is  at  most  second  order  in  the  largest  wave 
amplitudes  present  (see  LYK).  Also, 


<?«(i)  =  gh(x)  -  |n*w*/i2(i),  (4.8) 

We  now  Fourier  transform  the  wave  field  in  the  y  direction,  assuming  propagation  is  to 

be  considered  in  the  on-offshore  (±x)  direction.  We  consider  here  the  case  of  an  unbounded 
lateral  domain  and  a  wave  held  which  is  periodic  over  the  basic  interval  —b<y<b.  We 
then  represent  qn(*»y)  as 

E  (4.9) 

m=-M 

where 

A  =  ^  (4.10) 

as  before.  Substituting  (4.9)  in  (4.7)  and  neglecting  x  -  derivatives  of  small  terms  in  Gn  then 
leads  to  a  set  of  coupled  second-order  ODE’s  for  the  17”,  given  by 


^’C„  +  +  W)’C  +  imVAH'fcV  +  =  0; 


n  = 


(4.11) 


where  [n.f.t.]”  now  represents  triad  interactions  satisfying  resonance  conditions  in  i  and  y. 
Here,  k  is  the  wavenumber  determined  by  the  lowest  order  dispersion  relation  (4.1).  Also, 

(7r)’  =  nV-m*A»  (4.12) 
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For  fixed  n,ib,A,  large  values  of  m  will  make  7"  ima^nary,  wliich  corresponds  to  modes 
which  are  exponential  rather  than  oscillatory  in  x  in  the  linear  approximation.  In  the  linear 
case,  the  presence  of  these  modes  in  the  initial  conditions  would  be  interpreted  in  the  same 
light  as  the  presence  of  evanescent  modes  in  the  wavemaker  problem  discussed  in  Section 
2.1.1.  However,  the  interpretation  in  the  case  of  possible  nonlinear  forcing  of  the  offshore 
portion  of  trapped  modes  in  the  nearshore  region  is  non- trivial  and  will  need  to  be  considered 
carefully  in  applications  where  the  inclusion  of  this  effect  is  desired.  In  addition,  nonlinearity 
could  force  the  propagation  of  modes  that  would  not  be  present  in  a  linearized  wave  field, 
and  which  could  affect  a  detailed  representation  of  an  individual  wave.  At  present,  the  range 
of  M  at  each  value  of  n  may  be  restricted  to  Mn  <  nk/X  in  order  to  eliminate  forcing  of 
these  modes  arbitrarily. 


4.1  Shoaling  waves 


The  model  developed  in  the  previous  section  allows  for  the  onshore  and  offshore  propagation 
of  the  directional  spectrum  components.  Here,  attention  is  restricted  to  waves  propagating 
onshore,  or  in  the  positive  sense  with  respect  to  the  x  coordinate. 

Based  on  the  linear,  nondispersive  portion  of  the  model  (4.11),  we  assume  that  the 
incident  wave  may  be  written  in  the  form 

7"(x)  =  (4.13) 


where  it  is  assumed  that  the  x  dependence  of  A,  k  and  7  is  on  a  slow  scale  of  0(e),  and 
where 

=  =  ^  (4.14) 

(where  the  positive  root  is  taken).  The  amplitudes  A  represent  the  discrete  angular  spectrum 
being  considered  here,  and  are  allowed  to  vary  owing  to  refraction,  shoaling,  dispersion  and 
nonlinear  interaction.  (It  would  be  possible  to  absorb  shoaling  and  refraction  effects  by  the 
use  of  the  usual  linear  refraction  formulae;  this  step  is  not  taken  here.)  Substitution  of  (4.13) 
in  (4.11)  leads  to  the  spectral  model  for  incident  waves,  given  by 


j-m  jm  ,  (^^7n  )»  jm  _  jtn 


,._L  f  »»~1  ft 


fr.4r.r+ 


2H 


1=1  J>=/>3  J 

Ti  ***9  9  m  ^^91  f  •  ••  y  • 


=  0; 

(4.15) 


Here,  (•)*  denotes  the  complex  conjugate.  The  limits  of  summation  Pi  —  P4  are  ^ven  by 


Pi  =  max(-Afj,-Mn-i  +  m) 

Pj  =  min(M|,Afn-(  +  ”») 

P3  s  -  m) 

Pi  =  min(Jlf|,Jlfn.|.{  -  m)  (4.16) 
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The  interaction  coefficients  /  and  J  are  given  by 


I 


» 


I 


» 


cr 

7"».P 

•'n,/ 


1  +  Rf-rr-?  +  •  I*  +  <f 

(4.18) 


fm,-p 

*n,-l  • 


The  phase  arguments  0  and  T  represent  the  base  mismatch  in  the  x  direction  of  the  triads 
chosen  based  on  perfect  matching  in  y  and  t.  Generally,  the  only  components  which  expe¬ 
rience  complete  resonance  in  the  long  wave  limit  must  have  parallel  propagation  directions; 
all  obliquely  interacting  components  are  somewhat  detuned.  The  phase  arguments  are  given 
by 


©nf  =  /Hr+(n-/)fc7r-7-«Hr  (4.19) 

T™;”  =  e“L7  (4.20) 


The  spectral  model  (4.15)  is  a  set  of  coupled  first  order  ODE’s  which  are  solvable  by 
standard  techniques.  Results  were  obtained  using  a  4th-order  Runge-Kutta  scheme. 


4.2  Comparison  with  Laboratory  Data 

In  order  to  verify  the  basic  computational  model  provided  by  (4.15),  Kirby  (1990)  compared 
model  predictions  to  the  laboratory  data  obtained  by  Hammack  tt  al  (1990)  for  the  case 
of  glancing,  or  Mach,  reflection  of  a  cnoidal  wave  by  a  vertical  wall.  Additional  compar¬ 
isons  were  made  with  parabolic  models.  The  experimental  tests  were  conducted  using  the 
directional  wave  maker  at  the  Coastal  Engineering  Research  Center,  Vicksburg,  MS.  A  prior 
use  of  this  facility  to  study  the  properties  of  intersecting  cnoidal  waves  is  described  in  Ham- 
mack  et  al.  (1989),  referred  to  here  as  HSS,  who  also  discuss  the  instrumentation  and  data 
acquisition  used. 

For  the  tests  considered  here,  the  wave  basin  was  operated  with  a  water  depth  of  20cm 
in  a  constant  depth  region  extending  12.55m  in  front  of  the  wavemaker,  after  which  a  beach 
with  1:30  slope  provided  an  efficient  wave  absorber  giving  little  reflection.  For  the  Mach 
stem  tests,  two  parallel  false  walls  were  installed  perpendicular  to  the  wavemaker  axis  in 
order  to  provide  a  dosed  channel.  The  channel  walls  were  situated  13.26m  apart,  which  fixes 
the  width  of  the  numerical  domain  to  be  considered. 

The  generation  of  oblique  cnoidal  waves  using  the  directional  wavemaker  has  been  de¬ 
scribed  in  HSS.  In  the  present  laboratory  tests,  waves  were  initially  spedfled  as  having  a 
wavelength  of  2m  smd  a  crest  elevation  4cm  above  mean  water  levd. 

A  sample  of  the  modd  calculations  are  presented  here  as  gray  levd  contour  plots  of 
instantaneous  surface  over  the  modd  basin.  The  gray  level  plots  are  actually  of  the  quantity 
-drjidx,  and  the  pictures  thus  mimic  the  visual  image  that  would  be  obtained  in  an  overhead 
photograph  resulting  from  lighting  at  a  low  angle  from  the  direction  of  the  wavemaker.  (This 
is  similar  to  the  photographic  arrangement  in  HSS).  Figure  9  presents  results  for  the  test 
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Figure  9:  Reflection  of  Cnoidal  Wave  from  Barrier,  Corresponding  to  TEST  CR150204 
(Kirby.  1990) 


» 


CR15020-1.  which  clearly  shows  the  evolution  of  a  wide  Mach  stem  wave  along  the  reflecting 
boundary.  In  contrast,  Figure  10  shows  the  other  extreme  example  of  test  CR580204.  where 
the  angle  of  incidence  is  about  45°  and  the  reflection  pattern  is  regular  (i.e.,  nearly  Linear 
superposition.)  The  reader  is  referred  to  Kirby  (1990)  for  a  detailed  comparison  betw'een 
model  results  and  the  laboratory  data. 


5  Conclusions 

The  angular  spectrum  method  for  both  intermediate  and  shallow  water  depths  provides  a 
useful  tool  for  the  propagation  of  water  waves.  It  has  the  advantage  of  including  refraction, 
diffraction,  and  shoaling,  but  permitting  larger  angles  of  propagation  than  the  parabolic 
models  permit. 

The  angular  spectrum  provides  an  interesting  interpretation  of  diffraction.  For  the  case 
of  diffraction  through  a  gap,  treated  by  the  Kirchoff  approximation,  the  initial  condition 
contains  aU  the  information  for  diffraction.  The  circular  spreading  of  waves  behind  a  gap 
is  simply  the  radial  spreading  associated  with  the  different  propagation  directions  of  the 
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Figure  10:  Reflection  of  Cnoidal  Wave  from  Barrier,  Coi responding  to  TEST  CR580204 
(Kirby,  1990) 
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Fourier  modes.  There  is  no  coupling  between  the  modes  when  there  is  no  lateral  variation 
in  bathymetry.  Farther  for  cases  where  diffraction  occurs  within  the  computational  domain, 
say,  behind  a  shoal  surrounded  by  constant  depth  water,  the  coupling  between  the  bottom 
modes  and  the  surface  wave  modes  over  the  shoal  forces  new  modes  to  grow  which,  behind 
the  shoal,  lead  to  the  characteristic  focus  and  diffraction  regions.  However,  for  linear  wave 
models,  there  is  no  coupling  between  the  modes  after  the  waves  pass  over  the  (refractive) 
shoal  region. 

In  shallow  water,  the  mode  coupling  is  very  strong  and  plays  a  major  role  in  the  compu¬ 
tations.  However,  comparisons  to  data  show  that  the  angular  spectrum  model  is  very  good; 
in  fact,  better  than  parabolic  representations  of  the  Boussinesq  equations. 

Future  work  in  angular  spectrum  modelling  will  be  directed  to  the  development  of  an 
elliptic  angular  spectrum  model  to  provide  for  reflection  upwave  and  the  application  to  a  full 
directional  spectrum  by  superposing  many  frequencies. 
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